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ABSTRACT 


Four candidate navigation systems for the 
Space Shuttle Orbiter approach and landing phase are 
evaluated in detail. These include three conventional 
navaid systems and a single -station one-way doppler 
system. In each case a Kalman filter is assumed to 
be mechanized in the onboard computer, blending the 
navaid data with IMU and altimeter data. Filter state 
dimensions ranging from 6 to 24 are involved in the 
candidate systems. Comprehensive truth models with 
state dimensions ranging from 63 to 82 are formulated 
and used to generate detailed error budgets and sen- 
sitivity curves illustrating the effect of variations in 
the size of individual error sources on touchdown 
accuracy. The projected overall performance of each 
system is shown in the form of time histories of posi- 
tion and velocity error components. The detailed re- 
sults are summarized., compared and interpreted, the 
filter designs are reviewed and suggestions are made 
concerning possible software improvements. 
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1 . 


INTRODUCTION 


The reusable Space Shuttle Orbiter will be required to navigate 
accurately in multiple mission phases, based on data received from a 
variety of sources. Considerable effort is underway within NASA and 
among NASA contractors, aimed at producing reasonable solutions to the 
varices multi -sensor navigation problems involved. A best overall solu 
tion is lik ely to involve a commonalty of navigation hardware elements 
and software modules. This document describes the results of a 
program of research which is one element in that overall effort. 


1. 1 OBJECTIVES AND SCOPE 

The general purpose of this study is to aid in the evaluation and 
design of multi -sensor navigation schemes proposed for the Orbiter. The 
scope of the effort described herein is limited to the post-entry, energy 
management and approach and landing mission phases . Figure 1. 1-1 illus- 
trates one candidate system based on conventional navigation aids, in- 
cluding two DME (Distance Measuring Equipment) stations and ILS (In- 
strument Landing System) glide slope and localizer antennae . Some key 
elements of the system not shown in the figure are the onboard IMUs 
(Inertial Measurement Units), altimeters and a computer. The latter is 
programmed to mix together (filter) the IMU data and the externally- 
derived data 0 A completely automatic, all-weather landing capability 
is required. Since no air-breathing engines will be carried on orbital 
flights, tiier e will be no chance to "go around" and try again following 
a missed approach. 
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Specific objectives of this effort are: 

• Define a nd develop mathematical models describ- 
all potentially significant sources of error for 
each candidate navigation system eonsideredo 

• Develop a detailed, quantitative understanding 
of the contributions of individual error sources 
to overall system performance. 

9 Present results in a form which will help NASA 
choose or specify hardware elements (IMXJ, 
navaids) or their characteristics. 

• Evaluate alternative software approaches to the 
navigation filter design problem. 

Table 1. 1-1 summarizes the main features of four candidate systems 
which have been evaluated in detail. The first two, Systems A and B , 
consist of identical sets of hardware elements, but very different filters. 
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TABLE 1. 1-1 
CANDIDATE SYSTEMS 



IMli 

External Aids 

Filter 


System A 

KT-70 

2 DMEs 
Baro Altimeter 
ILS (3 deg) 
Radar Altimeter 

24-State 

Multi-Phase Filter 

No 

System B 

KT-70 

2 DMEs 
Baro Altimeter 
ILS (3 deg) 
Radar Altimeter 

6 -State 

Square-Root Filter 

No 

System C 

KT-70 

One-Way Doppler 
Baro Altimeter 
Radar Altimeter 

22-State 

Multi-Phase Filter 

No 

System D 

KT-70 

1 DME 
1 VOR 

ILS (Localizer Only) 
Baro Altimeter 
Radar Altimeter 

23 -State 

Multi-Phase Filter 

Yes 


The key hardware elements are the DME and ILS navaids pictured in 
Fig. 1. 1-1. System C s ubstitutes a one-way doppler beacon for the con- 
ventional navaids used in Systems A and B. System D is similar to 
System A except that a single VOR/DME (VHF Omnidirectional Range/ 
Distance Measuring Equipment) station is used instead of 2 DME stations 
and the glide slope portion of the ILS system is omitted. The principle 
computational results of this study are the error budget tables given in 
Section 5, one for each of these four candidate systems. 

The projected system navigation accuracies are compared in 
Section 5 with the specification values summarized in Table 1. 1-2. 

These values are taken from Schiesser (Ref. 1), who states that they 
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are preliminary numbers and are a function of vehicle dynamics, structure, 
runway size and braking. A basic groundrule, which defines allowable 
navigation error to be some fraction of allowable total error, is also in- 
volved. Although the values tabulated are likely to change, they provide 
a useful standard of comparison. 


TABLE 1.1-2 

LANDING NAVIGATION SPECIFICATIONS 



Allowable RMS Navigation 
Errors at Touchdown 

Vertical 

Downrange 

Crossrange 

Position (ft) 

3.0 

80.0 

mam 

Velocity (ft/sec) 

0.2 

3.0 

2.0 j 


1. 2 RELATIONSHIP TO OTHER EFFORTS 

An important foundation of this study is a family of computer 
programs (and the associated experience in using them) generated by 
TASC in connection with two Air Force sponsored programs. The first 
was a CIRIS (Completely Integrated Reference Instrumentation System) 
study in which TASC performed an independent evaluation (Ref. 2) of two 
filter designs. This called for the generation of system performance 
projections and detailed error budgets for each candidate system over 
a common trajectory and measurement geometry with respect to a. net- 
work of precision ranging transponders. The second was a CLASS 
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(Close Air Support System) 'weapon delivery system evaluation (Ref. 3), 
In both studies a high-dimensional mathematical truth model was 
formulated and the linear covariance technique outlined in Section 1.3 
was used to obtain the desired results. In the latter case (Ref. 3) it 
was possible to make detailed comparisons of computer -generated co- 
variance results with actual flight test results. Good agreement was 
was noted, leading to high confidence in the truth model and method- 
ology employed. 

Major inputs to this study have come from three types of efforts 
conducted by others. 

• Filter Design Studies 

e IMU Error Propagation Studies 

• Error Model Developments 

A number of Space Shuttle Landing navigation filter studies have been 
conducted or supported by NASA/JSC. Two which have directly impacted 
this effort are the continuing activities of TRW (Refs. 4, 5 and 10) and 
MIT Draper Laboratory (Refs. 6 and 18). Both of these have roots in 
previous work on the Apollo program. The filters evaluated as part of 
Systems A, C and D of this study are specializations of the 28 -state 
"multi-phase” filter described by Lear in Ref. 4. The 6 -state filter of 
System B is based on the square -root filter algorithms employed in the 
MIT Draper Laboratory studies. 

The starting point of the candidate system evaluations reported 
nerein is the end-of-bJackout point, which occurs at an altitude of 
130, 000 ft along a nominal reentry trajectory. Considerable use has 
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been made of data generated by Clark and Mitchell (Ref* 7)* This data 
represents initial navigation state errors and their correlations with 
basic IMU error sources, based on a shuttle "one-rev" mission with pure 
inertial navigation prior to entry. That effort made use of previously 
developed programs, notably the GEAP (General Error Analysis Program) 
described in Ref. 19. 

Error models representing the performance of IMUs and the 
various navigation aids considered have been the subject of numerous 
reports. A number of these are cited as references in the "truth model 
data base" tables of Section 3. A particularly useful NASA document, 
for our purposes, has been Ref. 9, which pulls together much data from 
a broad variety of sources and presents it in readily usable form. That 
document has in turn drawn much from Ref* 20 prepared by the Lockheed 
Electronics Company. 

The impact which this study is expected to have on other efforts 
involves both software and hardware decisions to be made by NASA or its 
major contractors. The software decisions involve methods for mixing 
inertial data with externally -derived data -- choice of filter states, filter 
update algorithms, etc. The hardware decisions involve the selection of 
an IMU and the devices needed to obtain the external data, as well as the 
specification of particular hardware characteristics, required survey 
accuracy, etc. The detailed, quantitative results and their interpretation 
provide valuable insights which are a useful input to the decision-making 
process. 
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1.3 APPROACH 

Figure 1. 3-1 is a generalized representation of all candidate 
navigation systems evaluated in this study. The hardware elements are 
an inertial measurement unit (IMU), a computer and a collection of de- 
vices (some airborne, some on the ground) needed to make external 
measurements. The computer function is divided into two parts, as shown. 
First, there is a standard set of algorithms needed to mechanize an iner- 
tial navigation system (INS). Second, there is a set of filter algorithms 
needed to blend in the external data. The filter is of the "feedback” 
variety, in which updated state estimates are used, after each measure- 
ment incorporation, to immediately correct appropriate variables of the 
INS rou tin e and the measurement routine. (The latter corrections take 
place only if the filter states include estimates of correlated measurement 

errors. ) 


EXTERNAL AIDS 


R-ao« 



STATE 

ESTIMATES 


Figure 1.3-1 Aided Inertial Navigation with a Feedback Filter 
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The overall methodology used in each candidate system eval- 
uation is summarized in Fig, 1.3-2. The upper half of the diagram 
represents the recursive solution of the filter covariance propagate and 
update equations. These are solved once for a given trajectory and 
measurement schedule. Certain elements of the filter dynamics and 
measurement matrices are functions of the Shuttle position and velocity 
vectors and the relative geometry between the Shuttle and the ground- 
navaid antenna locations. The outputs are the time histories of the filter - 
indicated performance and the Kalman filter gain matrices „ The latter 
are stored on tape and called a "gain file”. The lower half of the diagram 
represents the recursive solution of the linear system covariance equa- 
tions, given in Section 2„4 0 These are solved repeatedly to produce an 
error budget, the same gain file being used each time. The trajectory 
dependent matrices, describing the real-world error model are of much 
higher dimension than the corresponding filter matrices. In individual 
error budget runs specific elements of input matrices, corresponding to 
specific error sources, are set to non-zero values, with all other ele- 
ments set to zero. The output time history of the system error covari- 
ance matrix is then a statistical measure of the effect of that particular 
error source or small group of error sources, and generates one row of 
the error budget table. When the entire table is filled in, an overall sys- 
tem performance projection can be calculated from the detailed error- 
s ounce -by -err or -source breakdown. 

It should be mentioned that while the methodology outlined above 
is conceptually very simple, considerable complication could arise in the 
detailed computer programming involved, for two reasons. First, there 
is often a requirement to work with multiple coordinate frames and to 
transform various vectors and matrices from one frame to another in the 
course of accepting inputs, making calculations and expressing outputs. 
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The detailed mathematics given in Appendix A and Appendix B of this 
report are a good example of this type of complication. Second, the desire 
to hold down the cost of computer running time has led to numerous com- 
plexities and program options which take advantage of special properties 
of the system truth models. It is not necessary, for example, to work 
with 80 x80 matrices when the truth model state vector is 80 -dimensional. 
Because TASC had developed a considerable amount of software to do 
these calculations prior to the start of this study, the computer -running- 
time cost of the error budgets generated in this report is at least an order 
of magnitude less than would be the case if tne most straight-forward ap- 
proach to the methodology outlined above were followed. 
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Figure 1.3-2 Realistic System Performance Projections 

1.4 ORGANIZATION OF THE REPORT 

Section 2 describes the trajectory and measurement schedules 
used in evaluating the four candidate systems. Systems A, B and C were 
evaluated using identical trajectories. The same ground track was used 
in the System D study, but a different descent rate time history was used 
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during final approach., This modification reflects recent changes in the 
Orbiter aerodynamic configuration* Some basic terminology and the 
linear covariance equations used in evaluating candidates systems are 
also presented in this section* 

Section 3 presents the truth models used in evaluating the four 
systems. In each case the states and error sources are listed first; then 
the mathematical structure of the model is defined. Finally, the data base 
listing rms magnitudes, correlation times and data sources is given. 

Section 4 defines the filter states and algorithms used to com- 
pute filter -indicated performance and to generate the filter gain sequences, 
which were filed on magnetic tape. 

Section 5 presents the major results of the study. For each 
candidate system a detailed error budget table is given, listing the con- 
tributions of individual error sources, or small groups of error sources, 
to navigation errors at touchdown. Overall system performance projections, 
expressed as rms position and velocity component errors, are plotted 
from end -of -blackout to touchdown. Sensitivity curves showing the ef- 
fects of variations in the size of major error sources are presented* 

In Section 6 the results of the candidate system evaluations are 
summarized, compared and interpreted, and possibilities for improve- 
ment through software modifications are discussed* 

Section 7 provides a brief summary of major conclusions and 
recommendations for further study. The appendices (Volume n) contain 
mathematical details involving coordinate frames and transformations 
and detailed computer print-outs generated in the course of producing 
the system error budgets* 
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trajectory, measurement schedules 

AND METHODOLOGY 


This section describes the approach and landing trajectory, 
measurement geometry and measurement schedules used in evaluating 
the four candidate systems. A modified final-approach descent -rate 
history was used for System D, reflecting recent changes in the Orbit er 
aerodynamic configuration,, Section ‘LA defines the basic terminology 
and states the Unear covariance equations used in evaluating candidate 

systems . 


2„1 NOMINAL TRAJECTORY 

At the outset of the study JSC furnished a magnetic tape with 
position, velocity and other data representing atypical shuttle reentry 
and landing trajectory, starting at an altitude of 400,000 ft and ending 
at touchdown. A new tape was created covering the portion below 130, 000 
ft and including just that data needed for our purposes. Figure 2. 1-1 
shows the ground trace of the trajectory portion of interest and indi- 
cates several key events, such as the time of the first DMIE measure- 
ment, T dme , the time of the first baro altimeter measurement, T^, 
etc. One -minute time marks, indicating time elapsed after the first 
post -blackout measurement are included in the figure. Table 2.1-1 
defines the other symbols used in Fig. 2.1-1 and provides further de- 
tails, such as altitudes and relative velocity magnitudes at various times 

of interest. 
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Figure 2. 1-1 Trajectory Ground Trace 


2. 2 MEASUREMENT SCHEDULES 

Both the System A and System B evaluations are based on the 
assumption that filter updates occur every 2 seconds for the first 750 
seconds and every 0. 5 seconds for the last 46 seconds (after Tg^). At 
each update time a sequence of scalar measurements are incorporated 
according to the following order. 

1. DME Station 1 

2. DME Station 2 

3. Baro Altimeter ft s T BA } * 

4. ILS Localizer (t £ T LCL ) 


*q he baro altimeter measurements are inhibited in the interval 

950 < < 1020 ft/sec (near Marti 1). 

REJLi 
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TABLE 2ol“l 
KEY EVENT DATA 



System 

Time 

Altitude 

Relative 

Event 

Irajeciury 

Tape 

(sec) 

Elapsed 

(sec) 

(ft) 

Velocity 

(ft/sec) 

t dme or t owd’ First DME or 






OWD Measurement 

ABCD 

1636 

0 

130,000 

4,990 

T„ . , First Baro Altimeter 

JdA 

Measurement 






ABCD 

1806 

120 

100,000 

2,700 

T eM , Energy Management 






Initiation 

ABCD 

2092 

386 

50,000 

889 

T t ni , First Localizer 

JjLrL 

A'B 

2292 

606 

18,100 

696 

Measurement 

D 

2292 

606 

24, 000 

687 

T™,, Switch from 2 sec Update 

C 

2406 

720 

4,690 

495 

SW’ 

AB 

2436 

750 

1,870 

453 

Rate to 0.5 sec Update Rate 

D 

2436 

750 

4,958 

480 

Tj^, First Radar Altimeter 

ABC 

2430 

744 

2,450 

459 

Measurement 

D 

2450 

763 

2,530 

465 

T GLg , First Glide Slope 






Measurement 

AB 

2449.5 

762.5 

741 

440 

t stop’ stop Glide slope 






Measurement 

AB 

2474 

788 

85 

310 

T td , Touchdown 

ABCD 

2482 

796 

0 

302 
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5. Radar Altimeter (t £ T^) 

6. ILS Glide Slope (T GLg ^ t £ T STQp ) 

Therefore the measurement sequences are: 

Measurement Number Time 


1,2 

* t dme s 1 < t ba ) 

1,2,3 

< t ba s * " t lcl>* 

(except near Mach 1) 

1,2, 3,4 

< t lcl s 4 ' t ra> 

1,2,3, 4,5 

< T RA S 4<T GLS> 

1,2, 3, 4, 5, 6 

( T GLS s 4 s t stop* 

1,2, 3, 4, 5 

(t stop < 4 < t td^ 


The update rate changes after radar altimeter measurements have begun, 
but before the glide slope measurements begin. This is well before over- 
flight of the first DME station occurs. 

The System C measurement schedule includes filter updates 
every 2 seconds for the first 720 seconds and every 0. 5 seconds for the 
last 76 seconds (after T„„ 7 ). At each update time a sequence of scalar 
measurements are incorporated according to the following order: 


1 . 


2 . 


3. 


One-Way Doppler 

Baro Altimeter (t £ (except near Mach 1) 

Radar Altimeter (t a T^} 


Therefore the measurement sequences are: 
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Measurement Number Time 


1 

1 , 2 
1, 2, 3 


( T OWD 5 1 K t ba> 

(t . £ t <T d .) (except as noted above) 
BA RA 

( t ra * 1 < t td^ 


The update rate changes before radar altimeter measurements begin and 
well before the Shuttle overflies the one-way doppler station. 

The System D measurement schedule is very similar to that of 
Systems A and B. Updates occur every 2 seconds for the first 750 seconds 
and every 0. 5 seconds for the last 46 seconds. At each update time a 
sequence of scalar measurements is incorporated according to the follow- 
ing order. 

1. DME 

2. VOR* 

3. Baro Altimeter (t * T bA ) (except near Mach 1) 

4. ILS Localizer (t £ T LCL ) 

5. Radar Altimeter (t ^ T^) 

Therefore the measurement sequences are: 

Measurement Number Time 


1, 

2 

( T DME <t ba>* 

1, 

2, 3 

(T , £ t < T t ) (except as noted above) 
v BA LCL 

1, 

CO 

cL 

< T LCL £t<T RA> 

1, 

2, 3, 4, 5 

( T EA St<T TD> 


*The VOR measurements are inhibited in two regions when the elevation 
angle from the station to the vehicle is greater than 60 degrees — at the 
beginning of the spiral and during final approach,, 
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Figure 2.2-1 indicates the ground locations of the external aid 
devices. The values tabulated are incorporated in both the filter co- 
variance programs and the truth model subroutines. 


2.3 TRAJECTORY MODIFICATION FOR SYSTEM D STUDY 

The trajectory used to evaluate System D has exactly the same 
ground trace as that used in the earlier studies and pictured in Fig. 2. 1-1. 
A modified altitude history covering the final five minutes prior to touch- 
down has been employed in the System D study. The difference between 
the original trajectory and the System D trajectory is sketched in 
Fig. 2.3-1. A computer -generated listing describing the descent over 
the final minute was received from NASA. This corresponds to the per- 
formance of the new Or biter configuration and appears as the solid -line 
portion of the System D curve in the figure. The dashed-line portion was 
constructed to match the original trajectory at 5 minutes before touch- 
down and the new trajectory at 1 minute before touchdown, using the fol- 
lowing algorithm: 


Ah = k 2 r 2 -k 3 r 3 (2.3-1) 

where Ah is the difference between the new and the old altitude in 
Fig . 2.3 —1 , and 

t = i - tg (2.3-2) 

where t e is the time at five minutes prior to touchdown, and the param- 
5 

eters k 2 and k g are chosen such that 

^ = 0 1 at r = 0 
Ah = 0) 
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DME1 and VOR 


DME 2 


LOCALIZER ANTENNA 


GLIDE 5 LOPE ANTENNA 


OWD BEACON 

R Lp =6: 
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Figure 2.2-1 


R -7792 a 


N 

i 

m lci 


(2 GLS 
2 OWD 


m dme i 


UP 

EAST 

(ft) 

NORTH 

(ft) 

R L r 

200 

-9840 

r lp 

-9840 

0 

r lp 

0 

10,000 

r lp 

500 

-1300 

r lp 

115 

-3280 


373.3116 km 
0,909,773 ft 
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t (sec) 


Figure 2.3-1 System D Trajectory Altitude Profile 


Ah = 5174 ft ) 

\ at t = 240 sec 

Ah = -134 ft/sec) 

That is, coefficients of the polynomial of Eq, (2.3-1) were chosen to make 
the new trajectory match the desired altitude and altitude rate at both ends 
of the faired -in position. 


2o4 COVARIANCE EQUATIONS 

Detailed system error budgets are generated by repeated solu- 
tion of just two equations, the system covariance propagation and updat e_ 
equations, given at the end of this section. 
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Figure 2.4-1 is a conceptual representation of the mathematical 
structure of all truth models used in this study. An m-dimensional sys- 
tem error state vector, Sx a , propagates between measurement times 

* — D 

according to the linear, time-varying differential equation 


6x g (t) = F s (t)6x g (t)+w s (t) 


(2.4-1) 


where F g is the "system error dynamics matrix" and w g is a zero mean, 
gaussian noise pro .ith 


E Wg(t) w g (t - t ) 


Q s fi(t-T) 


(2.4-2) 


At discrete measurement times, t^, a measurement residual vector, 

6z D , is formed by 
— b 


5z q = H q 6x q + v q 
A S k ~ S k _S k 


(2.4-3) 


where Hg is the "measurement matrix" and v g is a zero mean, gaussian 
white sequence with 


E v Q v Q = R„ 
~k "^k S k 


(2.4-4) 


Filter corrections or updates are then given by 


c = K 6z q 
^k S k ^k 


(2.4-5) 


where K g is a gain matrix with m rows. 
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R-8047 



K g = AK p 




K m-n rows\ 
of zeroes ' 


(2„4-6) 


A is an m x n matrix linearly relating the n filter states to the m real- 
world states, and is the Kalman filter gain matrix (of n rows) pro- 
duced in the filter covariance calculation. The latter is based on the 
simplified n-dimensional model of system error propagation adopted by 
the filter designer. This model is defined by the filter matrices F^, Q F , 
H F and R^, which are analogous to the real-world or ''system" matrices, 

p Q n o and R c , discussed above. The filter and system error co- 
S ? S S b 

variance matrices are defined, respectively, as 


and 


P F (t) = E 


6x F (t) 5x F (t)' 


(2.4-7) 


P g (t) = E 


6xg(t) 6xg(t) J 


(2.4-8) 
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The error budget results given in Section 5 summarize system 
covariance calculations of the elements of P g (t). The major inputs to these 
calculations are the initial value, P g (0), and time-varying elements of F g , 
Q g , H g , R g and K g . Calculation of the non-zero elements of K g requires 
a solution for K_ {see Eq. (2.4-6), which implies a filter covariance calcu- 
lation. The major inputs to the latter are the initial value, P F (0), and the 
time-varying elements of Fp, Qp, Hp and Rp. 

The filter covariance equations used to generate the gain ma- 
trices, K„, can take one of a number of forms, such as the Joseph form, 
the square -root form, etc. The choice is made by the filter designer. 
Section 4 describes the specific forms used in this study. 

The system or truth model covariance equations are as follows: 
For propagation between measurements 



(2.4-9) 


where 



the transition matrix relating conditions just 
after the update at time to conditions just 
before the update at time =tj^ + At, (Fg is 
assumed constant over the interval, ) 


QSl 


At 

\ exp 


'0 


F Sk (At- r) 


Qg exp Fg. (At - r) dT 


the "discrete noise matrix" representing the 
integrated effect of system process noise Qg 
over the time interval from t^ to W* 
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Truncated matrix exponential series' are used in evaluating the trans- 
ition and discrete noise matrices defined above. For update of the 
system covariance matrix at the k^ 1 measurement time 



r i r i T 

p s ■ 


+ K s r s. K s_ 


k k k 


(2.4-10) 



The above is the Joseph form sometimes used in updating the filter co- 
variance matrix. This form is appropriate for any linear filter update 
and is therefore useful in sub -optimal filter evaluations. 

Equations (2.4-9) and (2.4-10) are time -varying, linear equa- 
tions in the system error covariance, Pg. This fact allows easy com- 
putation of overall system performance by root -sum -squaring separate 
contributions and easy development of the sensitivity curves. 


y 
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3. 


TRUTH MODELS 


This section describes the "truth models" used in evaluating 
Systems A, B, C and D.* Each truth model description includes a list 
of error sources, a detailed mathematical structure and a data base. 
The overall mathematical structure is a set of linear differential and 
algebraic equations describing system error propagation, as in Eqs. 
(2.4-1), (2.4-3) and (2.4-5). Most of the detailed structure is given 
in terms of sub-matrices of F g and H g , the system error dynamics 
matrix and measurement matrices, defined in Section 2 .4. The data 
base is the set of numerical values, such as rms values and correlation 
times, used to represent real-world error source statistics. 


3. 1 SYSTEM A TRUTH MODEL 

The hardware elements of System A include a KT-70 IMU, a 
computer and the equipment needed for six external measurements: range 
to two DME stations, baro altitude, radar altitude and ILS localizer and 
glide slope measurements. The ground station locations and the specific 
measurement sequence are given in Section 2. The computer navigation 
program includes a 24-state Kalman filter, called herein the JSC/TRW 
filter. The filter states and algorithms are defined in detail in Section 4. 1 
and Ref. 4. 

*By "truth model", we mean a mathematical model of all potentially 
significant sources of error and the way they affect system per- 
formance in the real world. 
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3. 1. 1 States and Error Sources 

The truth model states and other error sources used in evaluat- 
ing System A are listed in Table 3. 1-1, which divides them into three 
major categories: 

® The 24 estimated states and uncorrelated 
measurement and process noises. 

e Non-estimated states related to the inertial 
system. 

® Non-estimated states related to the external 
aids 

Error budget results corresponding to the first category above 
are generated using a truth model structure (Fg, Hg), which corresponds 
almost exactly to the filter error model structure (Fp, Hp). Results 
corresponding to the other two categories are generated using a higher - 
dimensional truth model, which contains the basic 24-state structure 
plus other states representing time -correlated error sources not modeled 
explicitly in the filter design. These additional error sources are divided 
into 16 groups as indicated in Table 3. 1-1. (The same group numbers 
are used in the Baseline Error Budget table in Section 5„1 and in the de- 
tailed results tabulated in Appendix C. ) 


3.1.2 Truth Model Equations 

The System A truth model requires a 24 x 24 filter dynamics 
matrix, Fp, and an 82 x 82 system dynamics matrix, Fg. The filter 
mechanization treats external data as a sequence of scalar measurements, 
rather than as components of vector measurements. Up to six separate 
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TABLE 3.1-1 


SYSTEM A TRUTH MODEL STATES AND ERROR SOURCES 


Error Source Names 

Number of 
States 

Number of 
Error Sources 

I. ESTIMATED STATES AND UNCORRELATED 



NOISES (Group 1) 



Position Errors 


3 

Velocity Errors 

3 

3 

Platform Misalignments 

3 

3 

Gyro Drifts (First-Order Markovs) 

3 

3 

Accelerometer Scale Factor Errors 



(First-Order Markovs) 

3 

3 

Accelerometer 'Biases' (First-Order 



Markovs 

3 

3 

Correlated Measurement Errors 



(First-Order Markovs) 



2 DMEs 

2 

2 

Baro and Radar Altimeters 

2 

2 

ILS Localizer and Glide Slope 

2 

2 

Uncorrelated Measurement Noise 



2 DMEs 

- 

2 

Baro and Radar Altimeters 

- 

2 

ILS Localizer and Glide Slope 

- 

2 

INS Quantization Noise 

- 

3 

II. NON-ESTIMATED, IMU -RELATED STATES 



Group 2. Accelerometer True Biases 

3 

3 

Group 3. Accelerometer Constant Scale 



Factor Errors 

3 

3 

Group 4. Accelerometer Misalignments 

6 

6 

Group 5. Accelerometer Nonlinearities 

3 

3 
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TABLE 3.1-1 (Continued) 


SYSTEM A TRUTH MODEL STATES AND ERROR SOURCES 


Error Source Names 


Number of Number of 
States Error Sources 


n. NON -ESTIMATED, IMU-RE DATED STATES 


(Continued) 

Group 6. Gravitational Deflections and 
Anomalies 

Group 7. Gyro True Bias Drifts 
Group 8. Gyro Mass Unbalances 
Group 9. Gyro Anisoelasticity 

HI. NON-ESTIMATED, EXTERNAL-AID 
RELATED STATES 

Group 10. DME Bias Errors 

Group 11. DME Scale Factor Errors 

Group 12. Baro Altimeter Errors 

Bias 

Scale Factor 
First-Order Markov 
Static Defect 

Group 13. ILS Bias Errors 

Localizer 
Glide Slope 

Group 14. ILS Second-Order Markov 
Errors 

Localizer 
Glide Slope 

Group 15. Radar Altimeter Errors 

Scale Factor Error 
Bias Error 

Group 16. DME Survey Errors 
Group 17. ILS Survey Errors 
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scalars may be incorporated at one update time. Therefore, there are 

six distinct filter measurement matrices, H„ through -- each is a 

■^1 ~ 6 

24 dimensional row vector. There are six corresponding system measure- 
ment matrices, H„ through H c -- each of these is an 82 dimensional 

C>i 

row vector. 


Figure 3. 1-1 presents the overall structure of the F g and H g 
matrices for System A. The upper-left partition of F g is the 24 x 24 
matrix F^ whose elements define the dynamic interaction between the 
Group 1 error states. This sub-matrix of F g corresponds* to the filter 
matrix, F^. The horizontal row of sub-matrices, F^ 2 through F 1 g , 
defines the effects of the non-estimated, IMU-related error sources 
(Groups 2-9) on the velocity errors (states 4-6) and the platform misalign- 
ment (states 7-9). The sub-matrices along the main diagonal, F 2 2 
through F 1? 17 , define the dynamics of all non -estimated, correlated er- 
ror sources,, For a group of random -constant error sources this sub- 
matrix is zero. The six system measurement matrices (row vectors) 
are outlined in the lower half of Fig. 3. 1-1. In each case the first 24 
elements are the same as those of the corresponding filter measurement 
matrix. Other non -zero sub -matrices, H. ., define the effects of the 

lj J 

correlated error sources (Group j) in question on the measurement error 
(6z^) in question. 

The sub -matrices F. . and H. . are defined in detail below, 

^ > J 3 

using the group -number designations given in Table 3. 1-1. 


F^ ^ and F^, are identical except for the presence of a 3 x 3 coordinate 
transformation. See T^p in Eq. (3. 1-1). 
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Figure 3. 1-1 System A Truth Model Structure 
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1=3x3 identity matrix 

r = position vector (in I frame; see Appendix A) 
= specific force vector in the X frame 




(3. 1-4) 



specific force vector in the P frame 


3x3 matrix which transforms vector in the P 
(platform) frame to the I{navigation error analysis) 
frame; see Appendix A. See also, the discussion of 
the A matrix, at the end of Section 3. 1.2. 


Group 2 •— Accelerometer True Biases, (3 sensors r 1 state each) 



"t 

1 I/P 

r 

0 

0 

5 2, 2 

6x3 



J3 x 3 


(3. 1-5) 


Group 3 — Accelerometer Constant Scale Factor Errors, 
(3 sensors, 1 state each) 


1,3 


1 

i 

GO 

I 

„ = 

L 0 . 

’ 3,3 

6x3 L 


J3 x 3 


(3. 1-6) 
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where 



(3. 1-7) 


Group 4 — Accelerometer Misalignments , (3 sensors, 2 states 

each) 



I 

p-i 

Eh 

i 

• TT = 

o 

« 

’ *4,4 


U 

6x6 

- 


(3. 1-8) 


6x6 


where 



0 


0 0 



0 


0 


0 0 



(3. 1-9) 


Gr oup 5 — Accelerometer Nonlinearities , (3 sensors, 1 state 

each) 





-16 x 3 



0 


J 3 x 3 


(3. 1-10) 
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where 


F 5 “ 


(3.1-11) 


Group 6 - Gravity Deflections and Anomaly , (up, east and north 
components) 


6x3 


1/t 


0 -l/r 0 (3. 1-12) 


0 -l/r. 


where 


t = d /V i 
u u ret 

T e ‘ V V rel 

T n ' V V rel 


and d u , d Q , d^ are correlation distances 


Group 7 - Gyro Bias Drifts, (3 sensors, 1 state each) 


= [i] 


; f, 


6x3 


= [ ° ] 3 X 


(3. 1-13) 
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Group 8 - Gyro Mass Unbalances, (3 sensors, 2 states each) 


0 



L 0 0 0 0 

1 2 

0 ° f 2 f 3 0 0 


L 


0 


0 


0 


0 



6 




J6 x 6 


(3. 1-14) 


The elements of F ± 8 and F ± g , below, are determined by the gyro input 
and spin axis directions shown in Fig. 3. 1-2. 


VERTICAL 
AT LAUNCH 





B-eou 


DOWNRANGE 

AT launch 
( 50.141 Degrcet 
Eon ol North) 


Figure 3. 1-2 Orientation of Gyro Axes 
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Group 9 — Gyro Ani so elasticity, (3 sensors, 1 state each) 




0 


¥2 

0 

0 

0 

f f 

V3 

0 

I 

0 

0 

f f 
2 3 



o] ( 3 . 1 - 15 ) 

-6x6 


Group 10 - DME Bias Errors, (2 stations, 1 state each) 


F 10, 10 



2x2 


(3. 1-16) 



(3. 1-17) 



(3. 1-18) 


Group 11 — DME Scale Factor Errors, (2 stations, 1 state each) 


11,11 


/T 


0 1 


sf 


0 - 1 / 


sf 


-12 x 2 


H l,ll ” P 1 
H 2,ll = P 2 


(3.1-19) 


(3. 1-20) 
(3.1-21) 
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where p^ and pg are the distances from the vehicle to DME stations 1 and 
2, respectively 

Group 12 — Baro Altimeter Errors, (4 states) 


12,12 


0 

0 

0 

0 


0 0 
0 0 
o -l/r a 
0 0 


0 

0 

0 

0 


bias 

scale factor 
markov 
static defect 


(3. 1-22) 


where 


H 3, 12 


1 


h 


1 



(3. 1-23) 


h = altitude 
^rel = relative velocity 

Group 13 — ILS Bias Errors, (2 states: localizer and glide slope) 


f„ 0 = r o 


(3. 1-24) 

13, 13 [ 

-2x2 

H 4, 13 ■ 

1 

(3. 1-25) 

H 6,13 = 

1 

(3. 1-26) 
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Group 14 — ILS Second Order Markov Errors , (4 states: 2 lo- 
calizer, 2 glide slope) 


'14,14 



0 

0 - 1/t lcl 



2 

~ 1/r GLS " 1/r GLS 

0 

1 

0 " 1/r GLS 


->4 x 4 


H 


4, 14 


= [1 0 0 0] 


(3. 1-27) 
(3.1-28) 


H, 


.g 14 = [0 0 1 0] (3.1-29) 

Group 15 — Radar Altimeter, (2 states: bias, scale factor) 


15,15 


0 


(3. 1-30) 


-*2x2 


H 5,15 = C 1 h ] t 3 ' 1 - 315 

Group 16 — D ME Su r vey Errors , (2 stations, 3 components each) 

(3. 1-32) 


F 


16,16 


0 


J6 x 6 


H, 


1,16 


■Pi Jh -PIq/pi ' p l n / p l 0 0 0 


(3. 1-33) 
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H 


■2,16 


0 0 0 " P 2 r / P2 ~ P 2 c / P2 ~ P 2 d / P2 


(3. 1-34) 


where 


Pi = 


Pi » P-i > P - 1 
R C D 


Pc 

— a 


p 9 7 Pc 7 Pc 

Z R v 


relative positive vector: station 1 to vehicle 
components of p^ in R frame 

relative position vector: station 2 to vehicle 
components of p^ in R frame 


each) 


Group 17 — ILS Survey Errors, (2 antenna locations, 3 components 


17,17 


(3. 1-35) 


-*6x6 


H 


4, 17 


V% 'V 


p ooo 

s 4 


H 


6,17 


0 0 0 -u /p„ 

0R/ 6 


-u /p -U /[ 

p c/ 6 P d/ 


(3. 1-36) 


(3. 1-37) 


where 


u , u , u 
e R e C e D 


components of a cross-runway unit vector; 
see Figure 3. 1-3. 

projection onto the horizontal plane of p^, the 
relative position vector from localizer antenna 
to the vehicle; see Figure 3. 1-3- 
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U 


P 


R 




components of a unit vector defined in 
Figure 3. 1-4. 


p = length of the relative position vector 
b from glide slope antenna to the vehicle; 
see Figure 3. 1-4. 


A Matrix 


The 82 x 24 matrix relating the filter states to the truth model 
states has the form: 


A = 


A' 

A 24 x 24 
0 


(3. 1-38) 


where 


A' 


69 18 24 


^6x6 


0 

0 

0 

T P/I 

0 

T p/I 

0 

0 

0 

0 

T p/I 

0 

0 

0 

T p/I 

0 

0 

h x 6 . 


(3. 1-39) 


The filter states 7 through 18, corresponding to platform misalignments, gyro 
drifts; accelerometer scale factor errors and accelerometer biases, are 
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defined in the I frame. The corresponding truth model states are de- 
fined in the platform (P) frame because a considerable amount of use- 
ful input data (from Ref. 7) represents the initial <h = 130,000 ft) 
contributions of individual sensors in the particular orientations pic- 
tured in Fig. 3. 1-2. Thus, the transformation Tp^j, from I to P 
coordinates, appears four times in the A matrix. 

The mathematical structure of the 82- state System A truth 
model and its relation to the 24 filter states has been completely defined 
in this section. The filter model is described in Section 4. 1. 


3.1.3 Truth Model Data Base 

In this section numerical values are assigned to the System A 
truth model matrix elements describing error source statistics. Most of 
these values are elements of the following matrices: 

P (0) — the initial system (real world) error covariance matrix 
S 

Q — the system process noise matrix 

® 2 
Re (i=l ... ,6) — the system measurement error variances, a. (in 

i ’ ” this case each R g is a 1 x 1 matrix, or scalar). 

l 

In addition to the above, some elements of F g , left unspecified in the 
previous section, are needed to complete the truth model. These are 
main-diagonal elements, such as the elements of Fg q, which involve 
correlation times of markov processes associated with various error 
source groups. These elements of F g and corresponding elements of 
P g (0)and Qg are normally chosen together to define random processes 
with desired properties. Frequently, a first-order markov process 
(exponentially correlated random process) is modeled with a specified 
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rms value, o., and correlation time, r.. In this case the relevant truth 
7 1 l 

model matrix elements are: 


p s <°> 

2 

= cr. 

l 

(3. 1-40) 

11 

F 

S.. 

-1/T i 

(3. 1-41) 

11 

Q s.. 

2 a ?/r. 

l i 

(3. 1-42) 


i. • 

li 


A constant error source can be considered a special case of the above 

with t. - oo and F„ = Q s = 0. 

1 b ii b ii 

Table 3. 1-2 summarizes the truth model rms values and 
correlation times used in evaluating System A. The data is organized 
according to the three major categories and 17 groups of error sources 
already used in Table 3. 1-1. The final column in the table indicates the 
pricipal references or sources of data. 

The 9x9 initial condition error covariance matrix is taken 
directly from page 203 of Ref. 7 (Table D-EH-c). This matrix re- 
presents error statistics at an altitude of 130,000 ft, following a "one- 
rev" mission in which pure inertial navigation is used from launch to the 
130,000 ft point. The assumed error sources correspond to KT-70 IMU 
performance, and are exactly those given in Table 3. 1-2 for groups 2, 3, 
4, 5, 7, 8 and 9. The assumed prelaunch alignment technique consists 
of gyrocompassing for azimuth alignment and accelerometer tilt leveling. 
This matrix is used in Section 5. 1 to produce one row in the "Alternative 
Error Contributions” table, showing how the initial errors, considered 
independent of root causes, propagate through the approach and landing 
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SYSTEM A TRUTH 


Error Source 


I. ESTIMATED STATES AND UNCOR- 
CORELATED NOISE (Group 1) 

Initial Condition Errors 
Position \ 

Velocity > 

Misalignment ) 

IMU Quantization Errors 


DME Measurement White Noise 

Baro Altimeter Measurement 
White Noise 

Radar Altimeter Measurement 
White Noise 

Radar Altimeter First Order 
Markov Error 

ILS First Order Markov Errors 
Localizer 
Glide Slope 

H. NON- ESTIMATED, IMU-RELATED 
STATES 

Group 2. Accelerometer Biases 

Group 3. Accelerometer Scale 
Factor Errors 

G~oup 4. Accelerometer Misalignments 
Group 5. Accelerometer Nonlinearities 
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TABLE 3.1-2 (Continued) 
SYSTEM A TRUTH MODEL DATA BASE 


Error Source 

Value 

Correlation 

Time 

(see) 

Data 

Source 

H. (Continued) 




Group 6. Gravity 

-4 9 



Anomaly 

<r= 13.3x10. ft/sec 

492,000/V r pt) 


Deflection about East 

cr = 5. 1 sec 

144, 000/V rel i 

Ref. 12 

Deflection about North 

a = 6. 5 sec 

98,50Q/V reI J 


Group 7. Gyro Bias Drifts 

o = 0. 01 deg/hr 

CD 

Ref. 9 

Group 8. Gyro Mass Unbalances 

cr=0. 015 deg/hr/g 

00 

Ref. 9 

Group 9. Gyro Anisoelasticity 

a = 0. 005 deg/hr/g^ 

CO 

Ref. 9 

HL NON-ESTIMATED, EXTERNAL-AID 




RELATED STATES 




Group 10. DME Bias Errors 

a = 295 ft 

CO 

Ref. 9 

Group 11. DME Scale Factor Errors 

cr = 100 ppm 

3600 

Ref. 13 

Group 12. Baro Altimeter Errors 




Bias 

a = 100 ft 

CO } 

Ref, 3 

Markov 

cr — 20 ft 

100) 


Scale Factor 

a = 0. 03 2 

CO | 

Ref 14 

Static Defect 

a = 1. 52 x 10“ 4 ft/V . 

rel 

CO / 


Group 13. ILS Bias Errors 




Localizer 

cr = 0. 5 mrad 

00} 

see text 

Glide Slope 

cr = 0. 5 mrad 

CD / 


Group 14. ILS Second Order 




Markov Errors 




Localizer 

see Eq. (3. 1-44) 

1.5*) 

Ref. 9 

Glide Slope 

see Eq. (3. 1-45) 

1.5*/ 


Group 15. Radar Altimeter Errors 




Bias 

a = 1 ft 

CO} 

see text 

Scale Factor 

o=0. 025 

CD 1 


Group 16. DME Survey Errors 

0 = 1 ft 

CO 

Ref. 15 

Group 17. ILS Survey Errors 

0 = 1 ft 

OD 

Ref. 15 


*These are values for ‘ r LClj and t gls In Eq. (3. 1-27) 
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phase. In the ’’Baseline Error Budget" initial position, velocity and align- 
ment errors are taken into account implicitly by proper treatment of 
accelerometer (Groups 2, 3, 4, 5)and gyro (Groups 7, 8, 9) errors. This 
is done by inserting the contributions to the initial errors of individual 
error so?irees, along with the error sources themselves. These contribu- 
tions are taken from the "1 sigma error tables” of Ref. 7. 

The most convenient way to express the effect of IMU quantiza - 
tion errors is in terms of the forcing function Q g of the velocity estimation 
error variances between updates. The value given in Table 3. 1-2 is based 
on the following typical values taken from Ref. 8, describing a current 
KT-70 application. 


One AV pulse = 1 cm/ sec 
Navigation cycle time = 0. 2 sec. 

Assuming that each AV count has a random error selected from a uniform 
distribution ranging from -0. 5 cm/sec to +0. 5 cm/ sec, the rms error is 

cr^y. = 1//I^ cm/ sec 
= 0. 00945 ft/ sec 

Assuming, further, that successive errors are uncorrelated, the addition 
of five such errors per second causes a navigation error growth rate of 

Qg = 5(0, 00945) 2 (i= 4, 5, 6) 

li -4 9 

= 4. 48 x 10 (ft/ sec) /sec 


as given in the table. 

Table 3. 1-2 does not list values for some of the Group 1 


J 


I 


l 
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(estimated) first-order markov states, which are included in the truth 
model mathematical structure. These are the gyro, accelerometer, 

DME and baro altimeter first-order states. The error budget given in 
Section 5. 1 is based on the assumption the gyro and accelerometer errors 
are random constants, as assumed in Ref. 7 and for Groups 2, 3, 4, 5, 

7, 8 and 9, and that DME and baro altimeter errors are represented by 
a more complicated structure, as in Groups 10, 11, 12 and 16. 


More complicated models are also used to represent ILS and 
radar altimeter errors, as in Groups 13, 14, 15 and 17. However, 
Group 1 states are included in Table 3.1-2 for these categories since the 
filter models are thought to be equally valid— results based on alternative 
models are given in Section 5. 1. The filter model for the radar altimeter 
correlated error takes the rms value to be the following function of altitude 
(units in feet): 



| 

, 0. 05 Alt (Alt > 500) 



°RA = ! 

j 0. 02 Alt (500 s Alt > 100) 

(3. 1-43) 


1 

i 2. 0 (100 * Alt) 


The filter models 

for the ILS correlated errors are based on (ct in millira- 

dians, p in feet) 

[5.6 

(p ■> 27300) 


cr LCL 

= 1. 4 + 4. 2 (f 

TSWr) (273 °° * p > 3500) 

(3. 1-44) 


*1.4 

(3500 > p) 
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2. 5 (p > 27300) 

1.4 + 1.1 frjgg) (27300 sp> 3500) (3.1-45) 

1. 4 (3500 ^ p) 

where p is the slant range from the antenna in question to the vehicle. 

Gravity deflections and anomalies are adequately modeled as 
first-order markov processes (Ref. 12) with fixed correlation distances. 
These are converted to equivalent correlation times by dividing by earth re- 
lative velocity, as shown in Table 3. 1-2, Group 6„ 

DME correlated errors are separated into bias, scale factor 
and survey errors — Groups 10, 11 and 16, respectively. The bias com- 
ponents represent the combined errors in calibrating both the airborne and 
ground equipment. The scale factor components represent the inaccuracies 
in calibrating the index of refraction. Values ranging from 10 parts per 
million to 100 parts per million have been mentioned in the literature, 
varying with the extent to which knowledge of local atmospheric conditions 
and measurement geometry is used in the calibration procedure „ Survey 
errors represent inaccuracies in knowledge of the ground transponder 
locations relative to the runway. Since, in this case, the devices are 
nearby, small survey errors are expected. 

Baro Altimeter correlated errors (Group 12) are separated into 
bias, markov, scale factor and static defect components. Values for the 
bias and markov components are taken from Ref. 3, which treats a 
baro-inertial-transponder system involving overflight of a transponder. 

In that case (the CLASS filter evaluation study) it was found to be impor- 
tant to include the markov state to account for moderately rapid changes in 
local weather conditions. Values for the scale factor and static defect 
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components are taken from Ref. 14, which discussed them in the context of 
of the Shuttle landing problem. 


It is not possible to cite references for the values used to represent 
ILS bias errors (Group 13) and radar altimeter bias and scale factor error 
(Group 15)= These components were included simply to gain insight into the 
problem -« that is, to discover the effect of such error sources, should 
they exist. The ILS second-order markov errors (Group 14) are modeled 
using the same rms values as given for the first-order components (Eqs. 

(3. 1-44) and (3. 1-45)). The process noise elements driving these states 
are computed using 

Qc, =4 cr/r. 
n3. . i i 

1,1 


Q 


S. . . 

i,i+l 


= Q 


S 5 


= 4 cr. 


1+-1 , i 


%. , . , = 4 Vi 
1 + 1 , 1+1 


(3.1-46) 


The data given in this section, together with the equations given 
in the previous section, completely define the System A truth model. 


3. 2 SYSTEM B TRUTH MODEL 

The hardware elements of System B are the same as those of 
System A; namely, a KT-70 EMU, a computer and the equipment needed 
for six external measurements: range to two DME stations, baro altitude, 
radar altitude and ILS localizer and glide -slope measurements. The ground 
station locations and the specific measurement sequence are given in Section 2. 
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The computer navigation program includes a 6 -state Kalman filter called 
herein the MIT-type filter . The filter states and algorithms are defined 

in detail in Section 4 = 2 and Ref. 6. 


3. 2. 1 States and Error Sources 

The truth mode) states and other error sources used in evalu- 
ating System B are listed in Table 3. 2-1, which divides them into three 
major categories: 

• The six estimated states (position and velocity 
errors), platform misalignments and uncor- 
related measurement and process noises. 

• Non-estimated error source states related to 
the inertial system. 

• Non-estimated states related to the external 
aids. 


Error budget results corresponding to the first category above 
are generated using a 9 -dimensional truth model structure (Fg, Hg), 
whose first six states correspond exactly to the 6 -dimensional filter 
error model structure (F F , %). Results corresponding to the other 
two categories are generated using a higher dimensional model, which 
contains the basic 9 -dimensional structure plus other states represent- 
ing time-correlated error sources. These additional error sources are 
divided into 17 groups as indicated in Table 3=2-1= The group number- 
ing system follows that used in the System A evaluation, and is used in 
the Baseline Error Budget table in Section 5.2 and in the detailed results 

tabulated in Appendix C . 
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TABLE 3.2-1 

SYSTEM B TRUTH MODEL STATES AND ERROR SOURCES 


Error Source Names 

Number of 
States 

Number of 
Error Sources 

I. ESTIMATES STATES, MISALIGN- 
MENTS AND UNCORRELATED 
NOISES (Group 1) ' 

Position Errors 

3 

3 

Velocity Errors 

3 

3 

Platform Misalignments 

3 

3 

Uncorrelated Measurement Noise 
2 DMEs 


2 

Baro and Radar Altimeters 

- 

2 

ILS Localizer and Glide Slope 


2 

INS Quantization Noise 

- 

3 

II. NON- ESTIMATED, IMU-RELATED 
STATES 

Group 2. Accelerometer True 
Biases 

3 

3 

Group 3. Accelerometer Constant 
Scale Factor Errors 

3 

3 

Group 4. Accelerometer Misalign- 
ments 

6 

6 

Group 5. Accelerometer Non- 
linearities 

3 

3 

Group 6. Gravitational Deflection 
and Anomalies 

3 

3 

Group 7. Gyro True Bias Drifts 

3 

3 

Group 8. Gyro Mass Unbalances 

6 

6 

Group 9. Gyro Anisoelasticity 

3 

3 
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TABLE 3. 2-1 (Continued) 

SYSTEM B TRUTH MODEL STATES AND ERROR SOURCES 


Error Source Names 

Number of 
States 

Number of 
Error Sources 

IE. NON- ESTIMATED, EXT NAL-AID 
RELATED STATES 

Group 10. DME Bias Errors 

2 

2 

Group 11. DME Scale Factor 

2 

2 

Errors 

Group 12. Baro Altimeter Errors 
Bias 

1 

1 

Scale Factor 

1 

1 

First-Order Markov 

1 

1 

Static Defect 

1 

1 


Group 13. ILS Bias Errors 
Localizer 
Glide Slope 

Group 14. ILS Second-Order 
Markov Errors 
Localizer 
Glide Slope 

Group 15. Radar Altimeter 
Errors 

Scale Factor Error 
Bias Error 
First-Order Markov 

Group 16. DME Survey Errors 

Group 17. ILS Survey Errors 

Group 18. ILS First-Order 
Markov Errors 
Localizer 
Glide Slope 

Totals 
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3.2.2 T ruth Model Equations 

The System B truth model requires a 6 x 6 filter dynamics 
matrix, F^, , and a 70 x 70 system dynamics matrix, F g . The filter is 
mechanized to handle external data as a sequence of scalar measurements. 
As in System A there are six-distinct filter measurement matrices, 
through Hj, — each is a six dimensional row vector. The six correspond- 
ing system measurement matrices, H g through Hg^, are 70-dimensional 
row vectors. 


Figure 3. 2-1 presents the overall structure of the F g and H g 
matrices for System B. The upper left partition of F g is the 9 x 9 
matrix, F x v whose elements define the dynamic interaction between 
the Group 1 error states. The upper left 6x6 partition of F^ 1 is identi- 
cally equal to the filter matrix, F y . The horizontal row of sub-matrices, 

F through F ± g , define the effects of the IMU-r elated error sources 
(Groups 2 -9) on the velocity errors (states 4-6) and the platform misalign- 
ments (states 7-9). The sub-matrices along the main diagonal, 2 
through F lg 18 , define the dynamics of all correlated error sources. The 
six system measurement matrices (row vectors) are outlined in the lower 
half of Fig. 3. 2-1. In each case the first six elements are the same as 
those of the corresponding filter measurement matrix. Other non-zero sub- 
matrices, H. ., define the effects of the correlated error sources (Group j) 
in question of the measurement error (fiz.) in question. 

The sub-matrices, F^p F 15> 15 > F is, 18’ H 5, 15’ H 4, 18 ^ H 6, 18 
are defined in detail below. All of the others are exactly the same as in System 

A, and are defined in Section 3. 1. 2. 


I 


I 
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Group 1 — Position, Velocity and Alignment Errors 



where 



/j, = gravitational constant 
1=3x3 identity matrix 

r = position vector (in I frame; see Appendix A) 


0 





0 




0 


( 3 . 2 - 1 ) 


( 3 . 2 - 2 ) 


( 3 . 2 - 3 ) 
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specific force vector in the P frame 


T =3x3 matrix which transforms vector in the P (platform) 
frame to the I (navigation error analysis) frame; see 
Appendix A. 


Groups 2-14, 16, 17 — See Section 3. 1. 2 — Systems A and B 
are the same. 


Group 15 — Radar Altimeter , (3 states: bias, scale factor, 
markov). In the System A truth model the markov state is included in 
Group 1 since it is estimated by the filter. 


15,15 


ro 


0 


0 


0 


0 


0 -1/r i 


J3 x 3 


(3. 2-4) 



h 



(3.2-5) 
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Group 18 — ILS First-Order Markov Errors (2 states: 1 local- 
izer, 1 glide slope). In the System A truth model these states are in- 
cluded in group 1 since they are estimated by the filter. 


F 18, 18 


'-l/r 4 0 

0 -1/t 6 


(3. 2-6) 


H 4>18 -[l o] (3.2-7) 

H 6;18 - [o l] < 3 - 2 - 8 > 


A Matrix 

The 70 x 6 matrix relating the filter states to the truth model 
states has the form: 


A 


r 6 x 6 
0 


(3.2-9) 


The material of this section, augmented by the definitions of 
most of the F and H sub-matrices given in Section 3. 1. 2, completely de- 
fines the mathematical structure of the 70-state System B truth model. 
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3. 2. 3 Truth Model Data Base 

The numerical values assigned to the System B truth model are 
exactly the same as those given in Table 3. 1-2. for System A. The values 
and expressions, Eq. (3.1-43), (3.1-44) and (3.1-45), given for the radar 
altimeter, localizer and glide slope first-order markov errors, which 
are listed under Group 1 in the table, are properly considered part of 
Groups 15 and 18 in the System B truth model. With these minor book- 
keeping changes Table 3. 1-2 is the System B data base. 


3. 3 SYSTEM C TRUTH MODEL 

The hardware elements of System C include a KT-70 IMU, a 
computer and the equipment needed for three external measurements: 
a cycle count indicating the integral of doppler shift over measurement 
intervals, baro altitude and radar altitude. The ground doppler beacon 
location is given in Section 2. The computer navigation program in- 
cludes a 22 -state version of the JSC /TRW multi -phase filter. The filter 
states and algorithms are defined in detail in Section 4.3 and Ref. 10, 

3. 3. 1 States and Error Sources 

The truth model states and other error sources used in evaluat- 
ing System C are listed in Table 3.3-1, which divides them into three 
major categories: 

« The 22 estimated states and uncorrelated 
measurement and process noises. 

» Non -estimated states related tc the inertial system. 

® Non -estimated states related to the external aids. 
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TABLE 3.3-1 

SYSTEM C TRUTH MODEL STATES AND ERROR SOURCES 


Error Source Names 

Number of 
States 

Number of 
Error Sources 




Position Errors 

3 

O 

Velocity Errors 

3 

3 

Platform Misalignments 

3 

0 

Gyro Drifts (First-Order Markovs) 

3 

3 

Accelerometer Scale Factor Errors 

3 

3 

(First-Order Markovs) 



Accelerometer "Biases’' (First-Order 

3 

3 

Markovs) 

Correlated Measurement 

Errors (First-Order Markovs) 
Doppler Rate Bias Errors 

2 

2 

Baro and Radar Altimeters 

2 

6 

Uncorrelated Measurement Noise 


1 

Dopuier u”"le quantization 

** 

9 

Baro and Radar Altimeters 

" 


INS Quantization Noise 

~ 

3 

II. NON -ESTIMATED, EMU-RELATED 
STATES 



Group 2. Accelerometer True Biases 

3 

3 

Group 3. Accelerometer Constant Scale 

3 

3 

Factor Errors 



Group 4. Accelerometer Misalignments 

6 

6 

Group 5. Accelerometer Nonlinearities 

3 

3 

Group 6. Gravitational Deflections and 

3 

3 

Anomalies 



Group 7. Gyro True Bias Drifts 

3 

3 

Group 8. Gyro Mass Unbalances 

6 

6 

Group 9. Gyro Anisoelasticity 

3 

3 
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TABLE 3.3-1 (Continued) 

SYSTEM C TRUTH MODEL STATES AND ERROR SOURCES 



Error Sources Names 

Number of 
States 

Number of 
Error Sources 

III. NON-ESTIMATED, EXTERNAL-AID 
RELATED STATES 

Group 10. Doppler Errors 

Rapid Shift Errors 

1 

1 

Antenna Motion Compensation 

1 

1 

Error 

Group 12. Baro Altimeter Errors 
Bias 

1 

1 

Scale Factor 

1 

1 

First-Order Markov 

1 

1 

Static Defect 

1 

1 

Group 15. Radar Altimeter Errors 
Scale Factor Error 

1 

1 

Bias Error 

1 

1 

Group 16. Doppler Beacon Survey Errors 

3 

3 

Totals 

63 

69 


Error budget results corresponding to the first category above 
are generated using a truth model structure (F g , H g ), which corresponds 
almost exactly to the filter error model structure (F p , H p ). Results 
corresponding to the other two categories are generated using a higher- 
dimensional truth model, which contains the basic 22 -state structure 
plus other states representing time -correlated error sources not modeled 
explicitly in the filter design. These additional error sources are divided 
into 12 groups as indicated in Table 3.3-1. (The group numbering system 
is consistent with that used in the System A evaluation, and is used in the 
Baseline Error Budget table in Section 5.3. Groups 11, 13, 14 and 17 

have been deleted . ) 


3-36 



I 


1 


1 


I 


I 


I I 


THE ANALYTIC SCIENCES CORPORATION 


3.3.2 Truth Model Equations 

The System C truth model requires a 22 x 22 filter dynamics 
matrix, F^, and a 63 x 63 system dynamics matrix, F g . The filter 
mechanization treats external data as a sequence of scalar measurements, 
rather than as components of vector measurements. Up to three separate 
scalars may be incorporated at one update time. Therefore, there are 
three distinct filter measurement matrices, H Fl through Hp 3 -- each is a 
22 dimensional row vector. There are three corresponding system measure- 
ment matrices, Hg^ through Hgg — each of these is an 63-dimensional row 
vector. 


Figure 3.3-1 presents the overall structure of the F g and H g 
matrices for System C. The upper-left partition of F g is the 22 x 22 
matrix F^ whose elements define the dynamic interaction between the 
Group 1 error states. This sub-matrix of F g corresponds* to the filter 
matrix, Fp. The horizontal row of sub-matrices, F^ 2 through F^g, 
defines the effects of the non-estimated, IMU-related error sources 
(Groups 2-9) on the velocity errors (states 4-6) and the platform misalign- 
ments (states 7-9). The sub-matrices along the main diagonal, F^ 2 
through F 1 R define the dynamics of all non-estimated, correlated error 
sources. For a group of random constant error sources this sub-matrix 
is zero. The three system measurement matrices (row vectors) are out- 
lined in the lower half of Fig. 3-1. In each case the first 22 elements 
are the same as those of the corresponding filter measurement matrix. 
Other non-zero sub-matrices, H. define the effects of the correlated 
error sources (Group j) in question on the measurement error (o z^ m ques- 
tion. 

*F and F„ are identical except for the presence of a 3 x 3 coordinate 
1 ? 1 F 

transformation. See T^p in Eq. (3.1-1). 
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Estimated 

States 


Non-Estimated, 
IMU-Related States 


Non-Estimated 
External -Aid 
Related States 


/ V 

1 0 

F l,2 

V 

• • ♦ 

V 



a) State Dynamics Matrix 


One-Way Doppler: 


Laro Altimeter: 


\ ' [ 
V [ 


Radar Altimeter: 


Vo 0 


H 2, 12 ( 


b) Measurement Matrices 
Figure 3.3-1 System C Truth Model Structure 


The sub-matrices, F^j, F 1010 , f i 6 ,16> H l,10’ H l,! 

H 0 and H„ are defined in detail below. All of the others 
la 

are exactly the same as in System A, and are defined in Section 
3.1.2. 


Group 1 — Estimated States (See also Section 4„ 3) 



(3. 


where 


D = 13 x 13 diagonal matrix with D Tii = -1/r. 

G , and the other elements appearing in Eq. (3.3-1) 
are defined in Eqs. (3.1-3) and (3.1-4). 


(3 
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Groups 2 through 9 - See Section 3.1.2 — Systems A and C 
are the same. 

Group 10 - Doppler Errors (2 states: rapidly shifting rate bias 
error and antenna motion compensation error) . 


'-1/t 


0 


10, 10 


0 -1/t 


(3.3-3) 


H 1 1Q = [At 1 ] (3.3-4) 

(See discussion of doppler errors in Section 3. 3., 3). 


Group 11 - Notin the System C truth model. 


Group 12 - Baro -Altimeter Errors - See Section 3.1.2 — 
Systems A and C are the same, except that the measurement sub -matrix 
is ^2 instead of Hg 


Groups 13 and 14 . - Not in the System C truth model. 


Groups 15 - Radar Altimeter Errors - See Section 3.1.2 
Systems A and C are the same, except that the measurement sub-matrix 
is Hg jg instead of Hg 
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Group 16 - Doppler Beacon Survey Errors (3 states : up, east, 

north). 


F 16, 16 



(3.3-5) 




(3.3-6) 


where h F , h F and h F are the R frame components of a vector defined 
in the I frame by h^ h 2? hg, the first three components of the filter 
matrix given in Eq. (4.3-6) . (Survey errors are input in the R frame. ) 

A Matrix - The 63 x 22 matrix relating the filter states to the 
truth model states has the form: 


A = 


22 x 22 



(3.3-7) 


where 


69 18 22 


V 

^6x6 


0 

0 

0 

T p/I 

0 

T P/I 

0 

0 

0 

0 

T p/I 

0 

0 

0 

T p/I 

0 

0 

( 

X 


(3.3-8) 


3-41 



I 


THE ANALYTIC SCIENCES CORPORATION 


The role of the transformation Tp^j in Eq. (3.3-8) is discussed at 
Eq. (3.1-39). 

The material of this section, augmented by the definitions of 
many of the F and H sub-matrices given in Section 3.1.2, completely 
defines the mathematical structure of the 63 -state System C truth model. 


3.3.3 Truth Model Data Base 

The general discussion of Section 3.1.3, pertaining to the 
System A data base, applies just as well to System C and is not repeated 
here. Also, most of the specific error sources, such as gyro and acceler- 
ometer errors, IMU quantization, baro altimeter and radar altimeter errors, 
are modeled the same way again; so those discussions are not repeated either. 
Table 3.3-2 summarizes the truth model rms values and correlation times 
used in evaluating System C. As indicated above, most of the values and 
references are exactly the same as those previously listed in Table 3.1-2. 

The error sources unique to System C are all associated with the one-way 
doppler measurements; these are discussed below. 

Altogether there are eight error sources associated with dop- 
pler measurements in the System C truth model. Three of them cor- 
respond directly to error sources modeled in the filter; these three are 
covered under Group 1 in Table 3. 3-2. They are measurement noise 
(treated as uncorrelated) due to the doppler cycle quantization and the 
two correlated rate bias errors associated with shifts in the onboard 
and ground beacon frequency standards. The baseline rms values and 
correlation times for these processes are taken to be the same as those 
assumed bv the filter. The value given for the onboard rate bias error 
corresponds to a clock error of 1 part in 10 -- see Ref. 10. 
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TABLE 3.3-2 


SYSTEM C TRUTH MODEL DATA BASE 


Error Source 

Value 

Correlation 

Time 

(sec) 

Data 

Source 

I. ESTIMATED STATES AND UNCOR- 
ELATED NOISE (Group l) 




Initial Condition Errors 
Position 
Velocity 
Misalignment 


9x9 Covariance 
Matrix (see text) 


Ref. 7 

IMU Quantization Errors 


Qs . 4.48 X 10- 4 ^' 2 
11 

(i = 4, 5, 6 -- see text) 

— 

Ref. 8 

Doppler Cycle Quantization 

a ~ 0.0357 ft 

— 

Ref. 10 

Baro Altimeter Measurement 
White Noise 

a = 5 ft 

— 

Ref. 3 

Radar Altimeter Measurement 
White Noise 

or - 3.28 ft 


Ref. 11 

Doppler Rate Bias Errors 

Onboard Frequency Standard 
Ground Frequency Standard 

q = 0.0876 ft/sec 
a = 0.0044 ft/sec 

0*0 
o o 

CO ^ 

Ref. 10 

Radar Altimeter First-Order 
Markov Error 

see Eq. (3.1-43) 

4.0 

Ref. 11 

II. NON-ESTIMATED, IMU- 

RELATED 




STATES 





Group 2. Accelerometer Biases 

o- = 60 ug 

00 

Ref. 9 

Group 3. Accelerometer Scale 
Factor Errors 

a ~ 34 ppm 

03 

Ref. 9 

Group 4. Accelerometer Misalignments 

rr = 40 sec 

03 

Ref. 9 

Group 5. Accelerometer Nonlinearities 

of = 3,5 ug/g 2 

00 

Ref. 9 
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TABLE 3.3-2 (Continued) 
SYSTEM C TRUTH MODEL DATA BASE 


Error Source 

Value 

Correlation 

Time 

(sec) 

Data 

Source 

II. (Continued) 




Group 6. Gravity 

-4 . 9 



Anomaly 

a = 13.3x10 re/sec 

492, 000/Vrel ) 


Deflection about East 

o = 5.1 sec? 

144, 000/Vrel [ 

Ref. 12 

Deflection about North 

a = 6.5 sec 

98, 500/Vrel ) 


Group 7. Gyro Bias Drifts 

<y = 0.01 deg/hr 

co 

Ref. 9 

Group 8. Gyro Mass Unbalances 

a = 0.015 deg/hr/g 

CD 

Ref. 9 

Group 9. Gyro Anisoelasticity 

a = 0.005 deg/hr/g^ 

CO 

Ref. 9 

IK. NON-ESTIMATED, EXTERNAL-AID 




RELATED STATES 




Group 10. Doppler Errors 




Fast Varying Frequency Standard 

a = 0.01 ft/sec 

30) 


Error 



See text 

Antenna Motion Compensation Error 

p£ = 0.5 ft, see Eq. (3.3-9) 

At) 


Group 12. Bar o Altimeter Errors 




Bias 

or = 100 ft 

00 ) 


Markov 

a = 20 ft 

1001 

Kei . o 

Scale Factor 

a = 0.03 „ 

oo 1 


Static Defect 

a = 1.52 xlO" 4 ft/V^ el 

“ / 

Kef. 14 

Group 15. Radar Altimeter Errors 




Bias 

o = 1 ft 

CO ) 


Scale Factor 

n = 0.025 

co} 

See text 

Group 16. Doppler Beacon Survey Errors 

cr = 1 ft 

00 

Ref. 15 
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Two other doppler measurement error sources are covered 
under Group 10 in Table 3. 3-2. The first is a rapidly -varying rate bias 
error assigned a correlation time of 30 seconds. No reference can be 
cited for this error source. It was included to gain insight into the prob- 
lem -- that is, to discover the effect of such an error source should it 
exist. (We know, for example, that a crystal frequency source might be 
subject to g-sensitive shifts. On the other hand, the rubidium device 
suggested in Ref. 10 is not expected to be subject to such errors.) The 
other source under Group 10 represents an error in compensating for 
antenna motion due to vehicle angular motion over the doppler measure- 
ment intervals. It is clear that compensation (via software) is required 
to account for the fact that the onboard doppler antenna is not located at 
the IMU, and will not experience the same relative motion with respect 
to the ground antenna whenever there is angular motion during a measure- 
ment interval. The inputs to the compensation calculation must include 
a pre -stored moment arm value, gimbal angle changes over the interval 
and a description of the geometrical relationship between the relative 
range vector (to the ground beacon) and the IMU axes. For error analysis 
purposes we use the following expression for the forcing function of this 
doppler measurement error source. 


Qg = (p£) 2 A0 2 /At 


(3. 3-9) 


where 


p = a percentage error in compensation 
H = IMU-antenna moment arm 
At = measurement interval 
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A0* = A 0 ^ + A02 + A 0 3 (3.3-10) 

A 0 A0 Afi = the three gimbal angle changes 
1 ^ during the interval, At. 

The value, p£ = 0. 5 ft, given in Table 3. 3-2 corresponds, for 
example, to a compensation error of 5 percent with a 10 ft moment arm. 
With the correlation time chosen as At, the effect of this model is to pro- 
duce an rms error of p^AS^ for the first time interval after the beginning 
of an angular motion. If the motion continues over successive intervals, 
the successive errors will be correlated to some extent. Aftur the motion 
ceases the error will die out quickly. 

Finally, three survey errors, representing inaccurate knowl- 
edge of the ground beacon location relative to the runway, are included 
in Group 16. Since the devices are nearby, small survey errors are 
expected. 


The data given in this section, together with that given in 
Section 3. 1, completely define the System C truth model. 


3. 4 SYSTEM D TRUTH MODEL 

The hardware elements of System D include a KT-70 IMU, a 
computer and the equipment needed for five external measurements: 
range and azimuth to a single VOR/DME station, baro altitude, radar 
altitude and ILS localizer measurements. The ground station locations 
and the specific measurement sequence are given in Section 2. The 
computer navigation program includes a 23 -state version of the JSC/ 
TRW multi -phase Kalman filter. The filter states and algorithms are 
defined in detail in Section 4. 4 and Ref. 4. 
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Initial condition errors are treated somewhat differently in 
the System D evaluation than they were in the preceeding three cases. 

The difference arises from the assumption, for System D, that drag 
acceleration is used to update the navigation state prior to the "initial” 
point at an altitude of 130,000 feet. This has no effect on the structure 
of the truth model, but does alter the details of the data base, as dis- 
cussed in Section 3. 4. 3. 

3. 4. 1 States and Erroi Sources 

The truth model states and other error sources used in evaluat- 
ing System D are listed in Table 3. 4-1, which divides them into three 
major categories: 

® The 23 estimated states and uncorrelated 
measurement and process noises. 

® Non-estimated states related to the inertial 
system. 

® Non-estimated states related to the external 
aids. 

Error budget results corresponding to the first category above 
are generated using a truth model structure (F g , H g ), which corresponds 
almost exactly to the filter error model structure (Fp, Hp). Results 
corresponding to the other two categories are generated using a higher- 
dimensional truth model, which contains the basic 23 -state structure 
plus other states representing time -correlated error sources not modeled 
explicitly in the filter design. These additional error sources are divided 
into 17 groups as indicated in Table 3. 4-1. (The same group numbers are 
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TABLE 3. 4-1 


SYSTEM D TRUTH MODEL STATES AND ERROR SOURCES 


Error Source Names 


Number of 
States 


Number of 
Error Sources 


L ESTIMATED STATES AND UNCORREL- 


ATED NOISES (Group 1) 


Position Errors 

Velocity Errors 

Platform Misalignments 

Gyro Drifts (First-Order Markovs) 

Accelerometer Scale Factor Errors 
(First-Order Markovs) 

Accelerometer 'Biases’ (First-Order 
Markovs) 

Correlated Measurement 
Errors (First-Order Markovs) 

DME 

VOR 

Baro and Radar Altimeters 
ILS Localizer 

Uncorrelated Measurement Noise 
DM r 
VOR 

Baro and Radar Altimeters 
ILS Localizer 

INS Quantization Noise 

H. NON-ESTIMATED, IMU-RELATED 


STATES 

Group 2, Accelerometer True Biases ’ 

Group 3. Accelerometer Constant Scale 
Factor Errors 

Group 4. Accelerometer Misalignments 

Group 5. Accelerometer Nonlinearities 

Group 6. Gravitational Deflections and 
Anomalies 

Group 7. Gyro True Bias Drifts 
Group 8. Gyro Mass Unbalances 
Group 9. Gyro Anisoelasticity 
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TABLE 3. 4-1 (Continued) 

SYSTEM D TRUTH MODEL STATES AND ERROR SOURCES 


Error Sources Names 

Number of 
States 

Number of 
Error Source;? 

III. NO’ T - ESTIMATED, EXTERNAL- AID 



RELATED STATES 



Group 10. DME Bias Error 

1 

1 

Group 11. DME Scale Factor Error 

1 

1 

Group 12. Baro Altimeter Errors 


1 

Bias 

1 


Scale Factor 

1 

1 

First-Order Markov 

1 

1 

Static Defect 

1 

1 

Group 13. ILS Bias Error 



Localizer 

1 

JL 

Group 14. ILS Second-Order Markov 



Error 



Localizer 

2 

i 

Group 15. Radar Altimeter Errors 



Scale Factor Error 

1 

i 

Bias Error 

1 

i 

Group 16. VOR/DME Survey Errors 

6 

6 

Group 17. ILS Survey Errors 

3 

3 

Group 19. VOR Bias Error 

1 

1 

Totals 

74 

81 


used in the Baseline Error Budget table in Section 5.4 and in the detailed 
results tabulated in Appendix C . There is no group 18 in the System D 
truth model. ) 
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3. 4. 2 Truth Model Equations 


The System D truth model requires a 23 x 23 filter dynamics 
matrix, and a 74 x 74 system dynamics matrix, F g . The filter 
mechanization treats external data as a sequence of scalar measurements, 
rather than as components of vector measurements. Up to five separate 
scalars may be incorporated at one update time. Therefore, there are 
five distinct filter measurement matrices, through ■ - each is a 
23 -dimensional row vector. There are five corresponding system measure- 
ment matrices, H c through H a -- each of these is a 74-dimensional row 
b l b 5 

vector. 


Figure 3. 4-1 presents the overall structure of the F g and H g 
matrices for System D. The upper -left partition of F g is the 23 x 23 
matrix whose elements define the dynamic interaction between the 

Group 1 error states. This sub-matrix of F g corresponds * to the filter 
matrix, F^. The horizontal row of sub-matrices, F^ g through F^ g, 
define the effects of the non -estimated, IMU-related error sources 
^Groups 2-9) on the velocity errors (states 4-6) and the platform mis- 
alignments (states 7-9). The sub-matrices along the main diagonal, 

Fg 2 through F^ ^g, define the dynamics of all non-estimated, cor- 
related error sources. For a group of random -constant error sources 
this sub-matrix is zero. The six system measurement matrices (row 
vectors) are outlined in the lower half of Fig. 3. 4-1. In each case the 
first 23 elements are the same as those of the corresponding filter mea- 
surement matrix. Other non-zero sub-matrices, H. define the effects 
of the correlated error sources (Group j)in question on the measurement 
error (6s^) in question. 


* 

F^ ^ and F^ are identical except for the presence of a 3 x 3 coordinate 
transformation. See T^ p in Eq. (3. 1-1). 
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The sub -matrix ^ is identical to that given by Eq. (3. 1-1) 
for System A except that the 24th row and column are omitted — since 
there is no glide slope measurement in System D. The sub-matrices 

F 10,10’ F ll,ll’ F !3, 13’ F 14, 14’ F 16,16’ F 17,17’ F 19, 19’ H l, 10’ 

H l,ll? H l, 16’ H 2, 16’ H 2, 19’ H 4, 13’ H 4, 14’ H 4, 17’ H 5, 16’ andH 5,19 
are defined in detail below. All of the others are exactly the same as in 

System A, and are defined in Section 3. 1-2. 


Groups 2-9, 12, 15 — See Section 3. 1-2 -- Systems A and D 
are the same. 


Group 10 — DME Bias Error {1 station, 1 state) 

(3. 4-1) 
(3. 4-2) 


F 10, 10 0 


H l,10 - 1 


Group 11 — DME Scale Factor Error (1 station, 1 state) 

(3. 4-3) 
(3. 4-4) 

where is the distance from the vehicle to the DME station. 


F ll,ll = 

H i,n = p i 


Group 13 — ILS Localizer Bias (1 state) 

(3.4-5) 
(3.4-6) 


F 13, 13 ~ 0 


H 4, 13 = 1 
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Group 14 - ILS Localizer Second-Order Markov Error (2 states) 


14,14 


“ 1//t lc l " 1//t lcl 


0 


-1/t 


LCL 


(3. 4-7) 


J 2 x 2 


H 


4,14 


[1 


0 ] 


(3. 4-8) 


Group 16 - DME and VOR Survey Errors (2 antennae, 3 com- 
ponents each) 


F 16, 16 W 6 x 6 


(3. 4-9) 


H 


1.16 


■Pi -P-t “Pi , 

i R/p 1 i C/p 1 i D/p 1 


0 0 0 


(3. 4-10) 


H 


2,16 


0 0 0 -u 0 -a, -u 9 

2 R/o C/p„ D/p. 


S 1 ' S 1 


(3.4-11) 


where 


p = relative position vector: VOR/DME 
station to vehicle; see Fig. 3.4-2. 


p ,p, , o.. - components of in R frame 

1 R X C X D 
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u ,u 9 , u„ = components of unit vector U 2 > P er P en ^ cu ^ ar 

2 R 2 C 2 D to and lying in horizontal plane; see 
Fig. 3.4-2. 

p = projection of p 1 onto horizontal plane. 

S 1 




Figure 3. 4-2 VOR/DME - Shuttle Geometry 


R-9623 


Group 17 — ILS Localizer Survey Errors (3 components) 


F 17,1'7 C°] 3 x 3 


(3.4-12) 


H 


4,17 


-u -u -u 0 0 0 

e R/p S C/p "D /P 

4 4 4, 


(3. 4-13) 
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where 


u ,u , u = components of cross -runway unit vector; 
e R e C e D see Fig. 3.1-3. 

p = projection onto horizontal plane of P 4 the 
s 4 relative position vector from localizer 
antenna to the vehicle; see Fig. 3. 1-2. 


Group 19 - VOR Bias (1 state) 


F 19, 19 


= 0 



(3. 4-14) 
(3.4-15) 


A Matrix 


The 74 x 23 matrix relating filter states to the truth model 
states has the form: 


A 


A ' 
23 


x 23 


(3.4-16) 


where A ' has the same structure given in Eq. (3. 1-39) for System A, 
except that the lower -right -hand corner is a 5 x 5 identity matrix rather 
than 6 x 6 . 


The material of this section, augmented by the definitions 
of common F and H sub -matrices given in Section 3. 1. 2, completely 
defines the mathematical structure of the 74-state System D truth 
model. 
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3. 4. 3 Truth Model Data Base 

The general discussion of Section 3. 1. 3, pertaining to the 
System A data base, applies just as well to System D and is not repeated 
here. Also, most of the specific error sources, including IMU errors, 
baro altimeter, radar altimeter, DME and localizer errors are modeled 
the same way again. Table 3.1-1 represents, therefore, most of the 
System D truth model data base. There are two ways in which the Sys- 
tem D data base is different: in its treatment of initial condition errors 
and in the inclusion of VOR errors; these are discussed below. 

As mentioned earlier it is assumed, for the System D evalu- 
ation, that drag-updating is employed prior to end -of -blackout to improve 
the navigation state estimate. This assumption has significant effects 
on both the initial filter covariance and the initial truth model errors. 
(The former is discussed in Section 4. 4. ) A recent study* at NASA/ 

JSC, involving 30 monte carlo runs through the black-out portion of 
a typical entry profile, produced a 6x6 covariance matrix r epresenting 
the position and velocity error statistics at end -of -blackout -- assuming 
the use of a drag -update scheme as proposed by Lear in Ref. 16. This 
covariance matrix represents the combined effects of IMU initial errors 
(at 400,000 ft, taken from Ref. 7), measurement noise and uncertainties 
in vehicle aerodynamic coefficients; it is one of the inputs used to gen- 
erate the System D error budget given m Section 5. 4. The same NASA 
study generated a measure of the effect of a non-standard atmosphere 
profile by inserting such a profile and holding it constant through, 'it the 
30 monte carlo runs. The result was a non-zero mean vector of position 

_ 

Summary results of the study were communicated to TASC by Messers „ 
P. Pixley and P. Michell of NASA. 
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and velocity errors. This vector is another input to the System D error 
budget calculations., The initial errors at 400, 000 ft, which were used 
in the study cited, correspond to the one-rev mission analyzed in Ref. 7 
with a KT-70 IMU and pre-launch alignment based on accelerometer 
leveling and optics azimuth alignment. In order to be consistent with 
the above we have treated the accelerometer errors (Groups 2-5) and 
gyro errors (Groups 7**9) by inserting the contributions to initial mis - 
alignments at 130,000 ft of individual error sources, along with the 
error sources themselves. These contributions are taken from the 
"1 a Error Tables” (Table D-IH-b) of Ref. 7. One additional error source 
is the initial pre -launch azimuth misalignment -- this did not appear in 
the System A, B or C evaluations because pre -launch gyrocompassing 
was assumed. 

A drawback of the procedure outlined above is that certain 
initial correlations, such as between position errors and misalignments 
at 130, 000 ft, are not properly accented for. (In fact .ere is no way 
to deduce them based on the information available to us. ) Fortunately, 
however, it happens that the missing correlations cannot have major 
effects in the case studied. This is shown by the result that the final 
contributions .it touchdown of all initial position and velocity errors are 
small (Table 5. 4-1). Thus, whether or not these quantities are highly 
correlated with other error sources is of little consequence. (Note that 
this result may not be the case for every system studied. ) 

Table 3.4-2 summarizes the features of the System D data 
base which are different from the System A data base given in 
Table 3.1-1. Altogether there are six error sources associated with VQR 
measurements in the System D truth model. Two of them correspond 
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TABLE 3.4-2 

SYSTEM D TRUTH MODEL DATA BASE 
(WHERE DIFFERENT FROM SYSTEM A) 


Error Source 

Value 

Correlation 

Time 

(sec) 

Date 

Source 

Group 1 

Initial Condition Errors 

Position 1 
Velocity f 

Azimuth Misalignment 

16 x 6 Covariance Matrix! 
| Mean Error Vector ( 

60 sec 

— 

see text 
Ref. 7 

VOE Measurement White 
Noise 

8. 5 mrad 

— 

Ref. 5 

VOR First-Order Markov 
Error 

18.67 rarad 

800 

Ref. 5 

Group 16 

VOR Survey Errors 

1 ft 


Ref. 15 

Group 19 
VOR Bias Error 

5. 0 mrad 

— 

Ref. 9 

Groups 1, 13, 14, 17 
Glide Slope Errors Omitted 

— 

— 



directly to error sources modeled in the filter. These are an exponentially 
correlated error and an uncorrelated measurement noise; they appear 
under Group 1 in Table 3. 4-1. Three error sources not modeled by the 
filter are the three VOR survey error components listed under Group 16. 
These are treated as distinct from the DME survey errors, since dif- 
ferent antennae are involved. The rms magnitude is taken to be one foot 
per component, as assumed for all survey errors considered in this report. 
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Finally, there is a state representing a VOR biag error. This is not 
modeled by the filter and is listed as Group 19 in the table. 

The data and equations given in this section, augmented by that 
given in Section 3. 1, completely define the System D truth model. 

Summary - The four system truth models have been pre- 
sented using a consistent format. The truth model data bases for 
Systems A and B are identical (except for minor bookkeeping differences) 
since these two systems have identical hardware elements. System C 
has a number of added features, all related to the one-way doppler mea- 
surements; the elements related to the DME and ILS measurements are 
deleted. The System D truth model is very similar to that of Systems A 
and B, with VOR elements added, glide slope elements deleted and initial 
condition assumptions altered. 


3-59 



I 

I 


1 


t 


t 


THE ANALYTIC SCIENCES CORPORATION 


4. FILTER COVARIANCE AND GAIN CALCULATIONS 


Two of the preliminary steps required prior to performing 
error budget calculations are the detailed definition of the filter which 
forms part of the system and the preparation of a filter covariance pro- 
gram. This program contains the covariance update and propagate algo- 
rithms specified by the filter designer and generates a sequence of gain 
vectors corresponding to a particular trajectory and measurement sched- 
ule. This sequence is saved on tape and repeatedly used by the truth 
model covariance program in generating the svstem error budget. Pre- 
sented below are the specific algorithms used to represent the filter gain 
calculations corresponding to the four candidate systems evaluated in 
this study. Systems A, C and D employ different versions of the 
JSC/TRW multi -phase filter. System B employs a low -dimensional square- 
root filter patterned after those studied by the MIT Draper Laboratory, 


4. 1 SYSTEM A FILTER 

Figure 4. 1-1 is a macro flow chart indicating the overall orga- 
nization of the filter covariance programs used for both System A and 
System B. The principal inputs are the starting time, T , and the initial 

filter covariance matrix, P„ . 

0 

The principal output is the file of gain vectors, through Kg 
for the first 303 update times (606 seconds) and through Kg for the 
remaining 164 update times (144 seconds with a 2 second interval and 
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46 seconds with a half second interval). At the time of the first localizer 
measurement (T = T lcj ) the current P p matrix can be punched on cards,, 
This gives us the capability to restart the problem at using these 

cards to represent the new P„ . 

0 

The large square blocks in Fig. 4. 1-1 contain the detailed 
covariance update and propagate algorithms. These can vary greatly 
from one filter design to another, depending on filter state size, state 
definitions, whether or not a square root formulation is used, reinitiali- 
zations, process noise treatment, etc. 

The JSC/TRW System A filter has been represented as a 24- 
state filter with the following state definitions. 


State Numbers 

Variables 

1-3 

Position Errors 

4-6 

Velocity Errors 

7-9 

Platform Misalignments 

10-12 

Gyro Drifts 

13-15 

Accelerometer Scale Factor Errors 

16-18 

Accelerometer Bias Errors 

19 

DMEl Correlated Error 

20 

DME2 Correlated Error 

21 

Baro Altimeter Correlated Error 

22 

ILS Localizer Correlated Error 

23 

Radar Altimeter Correlated Error 

24 

ILS Glide Slope Correlated Error 


The upper-left 18 x 18 submatrix of P has been supplied to 

*0 

TASC by JSC; the numbers are not repeated here. The last six rows 
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and columns are all zeros except for the six diagonal elements. These 
are the squares of the initial rms values given or indicated in the "cr c ” 
column in Table 4. 1-1. The transition matrix elements are programmed 
in accordance with the description given in Ref. 4 (see pages 46-53). 

TABLE 4. 1-1 


FILTER STATE MEASUREMENT ERROR STATISTICS; SYSTEM A* 




Correlated Error (Filter State) 

Measurement 


V. jusoromont 



Noise 


CT c 

7 


1 . 

DM El 

295 ft 

400 sec 

48 ft 

2. 

DME2 

295 ft 

400 see 

48 ft 

3. 

Baro Altimeter 

Eq. (4.1-10) 

40 sec 

Eq. (4.1-11) 

4. 

ILS Localizer 

Eq. (4.1-13; 

0. 55 sec 

0 

5. 

Radar Altimeter 

Eq. (4.1-12) 

4 sec 

3. 28 ft 

L_ 1 

6. 

ILS Glide Slope 

Eq. (4.1-14) 

0. 55 sec 


The last six diagonal elements of the filter transition matrir 
have the form: 


1 - At/(2 r.) 

0 F.. = 1 + At/(2 r.) ;i = 19,...,24 (4.1-1) 

n l 

Data from Ref. 5 and personal communication with W. L«ar of TRW. 

* ^ 

Except where 2 r. < At; see Eq. (4. 89) Ref. 4. 
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where At is the propagation step size. The corresponding discrete 
process noise elements have the form 




(4.1-2) 


where the r. are the correlation times given in Table 4. 1-1. Each filter 
covariance propagation employs the following standard formula: 



(4. 1-3) 


At each update time a sequence of covariance updates are cal- 
culated, corresponding to a sequence of scalar measurements. In each 
case a 24 x 1 gain matrix is calculated using 

K f = p; h£/(H f P- h£ u 70 + % 2 ) (4. 1-4) 

where is a 1 x 24 measurement matrix, is the covariance just 
prior to the updating, u^^ is a scalar measurement underweighting 
factor , * and is the assumed measurement noise variance (see the 
last column of Table 4. 1-1). The covariance matrix is then updated 
in accordance with 

= (I-K f H f )P f (4.1-5) 


* 

The scalar u^q is normally unity, but has the value 1. 2 during the 
early portion of the trajectory (see Refs. 4 and 17). 


CL '2- 
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The 1x24 matrices can be viewed as 24-dimensional row 
vectors, h. The non -zero elements of h for each of the six scalar mea- 
surements are given below: 


DME1 

DME2 


hj = p x /p 

h l = p x^ p 


h 2 =Py/p 

h 2 = p y ^ p 


h 3 = p z /p 

hg = Pgf / P 


h 19 =1 

h 20 = 1 

(4.1-6) 

where p is the magnitude of the vector displacement of the vehicle from 

the DME station in question, and p x> 

p and p are its components in 
y z 

the I frame (see Appendix A). 



Baro Altimeter 

Radar Altimeter 


h i - x/r 

h^ = x/r 


h 2 =y/r 

h 2 = y/r 


h 3 = z / r 

h 3 = z/r 


*21 = 1 

h 23 = 1 

(4. 1-7) 


In the above r is the magnitude of the current position vector from the 
center of the earth and x, y and z are its components. 


ILS Localizer 
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( 4 . 1 - 8 ) 


In the above p s ^ is the magnitude of the projection of the displacement 
vector p 4 onto the earth's surface, as shown in Fig. 3.1-3. The quan- 
tities u , u , and u are components of the unit vector u g shown in 

e x e y e z 

the figure. 


ILS Glide Slope 

h i ■ V° e 

h 2 = U p/ D 6 

h 3 = U p z /p 6 

h 24 = 1 (4. 1-9) 


hi the above p 0 is the magnitude of the displacement vector shown in 

Fig. 3. 1-3, and u , u and u are components of the unit vector u . 
Px Py Pz 

The rms values of the correlated errors associated with the 
last four external-aid devices are programmed as a function of position 
and (in one case) relative velocity. The baro altimeter measurement 
noise 1-sigma value is also given as a function of altitude. These func- 
tions are indicated below. 

Baro Altimeter , (units in feet and ft/sec) 

ct 2 = 17,222 + 3. 99 x 10" 14 Alt 4 + 3. 22 x 10" 8 V 4 0l 
0 

ct 2 = 9.84 + (0.00025 Alt) 2 
m 


(4. 1-10) 
(4.1-11) 
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Radar Altimeter, (units in feet) 




r. 


0. 05 Alt 
0. 02 Alt 
2.0 




(Alt > 500) 
(500 > Alt >100) 
(100 > Alt) 


ILS Localizer . (g in milliradians, p in feet) 


cr 

c 


5.6 




(p> 27,300) 
1.4 + 4.2(^^°) (27,300>P>3,500) 


1.4 


v 


(3,500 2: p) 


(4. 1-12) 


(4. 1-13) 


ILS Glide Slope, (cr„ in milliradians, p in feet) 
■ ' 11 " * ' *" 


CT 

c 


r 


- \ 


(p > 27,300) 

+ (27,300 >p> 3,500) 

(3,500> P) 


2.5 
1 


1.4 




(4. 1-14) 


The filter covariance program based on the 24-state filter out- 
lined above has been exercised over the trajectory and measurement sched- 
ule described in Section 2. The time history of the filter -indicated per- 
formance is given in Table 4. 1-2. RMS values of position and velocity 
errors in R frame components (runway coordinates; see Appendix A) are 
given at key times — just before and just after the first measurement of 
each type, and at touchdown. At the bottom of the table the filter-indicated 
rms estimation errors, at touchdown, for the other 18 filter states are 
also given. The results given in Table 4. 1-2 represent what optimal per - 
formance would be if this 24-state filter model were a completely accurate 
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TABLE 4.1-2 

SYSTEM A FILTER-INDICATED PERFORMANCE 


i i 

u 


i ; 

u 

n 


Event 

(see Table 2.1-1) 

Time 

(sec) 

RMS Position Errors 
(R frame (ft)) 

RMS Velocity Errors 
(R frame (ft/sec)) 

Vertical 

Downrange 

Crossrange 

Vertical 

Downrange 

Crossrange 


1686' 

49,400 

03,100 

25,000 

160 

93.6 

134 


1686 

15,300 

25,900 

20,200 

40.8 

57,2 

64.1 

t ra 

1806" 

2,920 

1,170 

874 

6, 71 

5. 79 

7.69 


1800* 

2,080 

1,030 

817 

5. 12 

5. 73 

6. 94 

T 

2292" 

473 

208 

509 

1.69 

1.98 

3.98 

LCL 

2292 

470 

203 

416 

1.65 

1.74 

3.03 

T 

2430' 

655 

231 

53.1 

1.77 

1.26 

0. 92 

RA 

2430* 

122 

213 

50.4 

0. 58 

1.14 

0.89 

T 

2449. 5' 

43.4 

218 

16.9 

0.49 

1.15 

0.62 

GLS 

2449. 5* 

24.8 

206 

16.3 

0.47 

1.12 

0.62 

T 

a TD 

2482 

1.9 

30.6 

5.5 

0.16 

0.90 

0.34 


OTHER INDIC ATED ESTIMATION ERRORS AT TOUCHDOWN, T. 


Misalignments (rad) 

Gyro Drift (rad/sec) 
Accelerometer Scale Factors ~ 
Accelerometer Hlasos (ft/sec 4 ) 
DME Errors (ft) 

Baro Altimeter Error (ft) 
Localizer Error trad) 

Radar Altimeter Error (ft) 
Glice Slope Error (rad) 


TD 

— 0.287x10-3, 0.079xl0“ 3 , ai23xl0‘ 3 ' 

— 0.279x10*6, 0.2B3xl0- 6 , 0.278xl0 -8 ’ 

— 4.94 xlO-5, 4.55 xlO‘5, 4.90 xl0-S| 

— 5.67 xlO-4, 5.60 xlO' 4 , 5.66 XlO* 4 ' 

— 41.4, 29.2 

— 2.5 

— 1.54 xl0“ 3 

— 1.8 

— 1.4 xlO’ 3 


representation of real-world error dynamics. The results given in 
Section 5.1 are a projection of the performance of this 24-state sub- 
optimal filter in an 82 -state model of the real world. 


* 

The rms misalignment estimation errors are in runway coordinates 
(Vertical, Downrange, Crossrange); the gyro and accelerometer errors 
are given in I frame (X, Y, Z) coordinates — see Appendix A„ 
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4. 2 SYSTEM B FILTER 

An MIT -type square root filter used in System B has been 
represented as a 6-state filter with the following state definitions. 

State Numbers Variables 

1-3 Position Errors 

4-6 Velocity Errors 

The square root matrix W is defined by 

WW T = P^ 

F 

For this study the following initial diagonal matrix has been used: 


a 

P 


where 



o = 20,000 ft 
P 

a = 20 ft/sec 
v 


(4. 2-1) 
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The filter error dynamics matrix is 


process noise matrix is 


Q, 


0 

I 

G 

0 

r 


lent matrix defined 

0 

0 

0 

Q 


( 4 . 2 - 2 ) 


J6 x 6 


( 4 . 2 - 3 ) 


J6 x 6 


where 


Q = 


0 0 


o q. 


v 


0 0 q 


( 4 . 2 - 4 ) 


q v = 0. 01 (ft/sec) /sec 


(See the dicussion of the choice of q y at the end of this section. ) 

The propagation of W between measurements is based on the 
following equation (from Ref. 6). 
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W = F F W + i QyW’ 1 


(4. 2-5) 


where is a ’'diagonalized" form of W. The "average G" numerical 
integration procedure (Ref. 18) is mechanized as follows: 


W(t + At) = < I + At 


0 

I 

At 2 

G r (t) 

0 

G 

avg 

0 

+ 2 

0 

G r (t+ At) 


W(t) 


J) 


where 


G_ 


+ 


At 

0 

At . 

2 Q 

2 

0 

Q 




G p (t) + G r (t+At) 


avg 


(4. 2-6) 


(4. 2-7) 


To update W at the i^ scalar measurement incorporation the 
following equations are used (see Ref. 6). 


z. = W T H f T (i = 1, — ,6) 

i 


(4. 2-8) 


D. = z z+R-p (a scalar) 
1 r . 

1 


K. = Wz./D. 
i l i 


(4. 2-9) 
(4. 2-10) 
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w 




w - 


K. z.' 
1 1 




(4.2-11) 


The K. are the 6-dim ensional gain vectors making up the gain file shown 
in Fig. 4. 1-1. This file is used to produce the System B error budget re- 
sults given in Section 5. 2. Expressions and values for the measurement 
matrices (6-dimensional row vectors) and measurement error vari- 

ances R_, are given below. The H_ have the form 
F i F i 


H 


F. 

l 


h l h 2 h 3 ° ° ° 


(4. 2-12) 


The three non-zero elements, hp hg and hg, are the same as the corres- 
ponding System A elements, given in Section 4. 1. For the DME measure- 
ments (i = 1 and 2): 


DME 1 

DME 2 


h l = p x /p 

h l = p x /p 


h 2 = p /° 

h 2 = 0y /p 

(4. 2-13) 

h 3 = P J P 

h 3 = °z /o 



where p is the magnitude of the vector displacement of the vehicle from 
the DME station in question, and p , p , and p z are its I frame 
components. 
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For altimeter measurements ( i = 3 and 5): 


Baro Altimeter 

Radar Altimeter 

h^ = x/r 

\ = x/r 

h 2 = y/r 

hg = y/r 

h 3 = z/r 

h 3 = z/r 


(4. 2-14) 


where r is the magnitude of the current position vector, from the center of 
the earth, and x, y and z are its I frame components. 

For ILS measurements (i = 4 and 6): 


Localizer Glide Slope 



The unit vector u and p S4 are defined in Fig. 3.1-3. The unit vector 
u and p„ are defined in Fig. 3.1-4. 

-p -6 

The six filter measurement error variances R F (i = 1, . . . , 6), 

used in the System B filter covariance calculations were computed in the 

following way: In each case the values used in the System A model for 

both correlated errors, tx , and uncorrelated errors, cr were combined 

c 

to yield 
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2 



(4. 2-16) 


The values or expressions indicated in Table 4. 1-1, for System A, were 
used in Eq. (4. 2-16) for System B. 

MIT personnel suggested the following range of possible values 
for the velocity process noise: 

q v = 0. 01 to 0. 10 (ft/sec ) 2 /sec 

The minimum and maximum values in this range were tried in the 
System B filter covariance program. The smaller value produced 
a significantly better filter -indicated vertical velocity performance at 
touchdown, so the gain values corresponding to that run were used to 
generate the System B error b idget. The time history of the filter- 
indicated performance is given in Table 4. 2-1. (The final values 
corresponding to the larger process noise are also given, at the bottom 
of the table. ) These results present what optimal performance would be 
if this 6-state filter model were a completely accurate representation of 
real-world error dynamics. The results given in Section 5. 2 are a pro- 
jection of the performance of this 6-state sub-optimal filter in a 70-stp ^e 
model of the real world. 

Note on Subsequent MIT Studies — Since the completion of the 
System B evaluation, MIT personnel have conducted further studies cor- 
responding to similar measurement sequences. These have included 
trail-and-error procedures to achieve better choices of process noise 
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TABLE 4.2-1 

SYSTEM B FILTER-INDICATED PERFORMANCE 


Event 

(see Table2.1-1) 

Time 

(sec) 

RMS Position Errors RMS Velocity Errors 

(R frame (ft)) (R frame (ft/sec)) 


Vertical 

Downn age 

Crossrange 

Vertical 

Down range 

Crossrange 

t dme 

n 

20,000 

19,800 

20,000 

11,500 

20,000 
10, 500 

im 

■■ 

— 

t ba 

e 

1 2 E I 

2,900 

1,100 

1,790 

894 

mt 


14.9 
10. 9, 

t lcl 

m 

42.5 

41.1 

70.7 

67.2 

184 

174 


wm 

HI 

t ra 


35.0 

32.6 

50.3 

49.2 

50.5 

48.4 


0.66 
0. 64 

m 

t gls 

2449. 5' 
2449. 5 

11.3 

8.4 

38.3 

37.2 

AES 


M 


t td 

(q v = 0. 01) 

2482 

1.4 

21,3 

4.1 

0. 10 

0.46 


t td 

K = 0 ' l) 

2482 

1.9 

29.0 

m 

0.43 

1.30 

0. 64 


anrl measurement noise variances, as well as the use of additional filter 
states such as level misalignment and DME bias error estimates. One 
of their conclusions is that better performance is possible with a six- 
state filter than with the six-state version described above. It is felt 
that the general comparison and summary discussion given in Section 6 
of this report remains valid. In addition, our suggestions concerning the 
use of additional states seem to be borne out by the results of the MIT 
experiments along these lines. 


4.3 SYSTEM C FILTER 

The JSC /TRW System C filter has been represented as a 
22 -state filter with the following state definitions. 
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Variables 

Position Errors 
Velocity Errors 
Platform Misalignments 
Gyro Drifts 

Accelerometer Scale Factor Errors 
Accelerometer Bias Errors 
Onboard Doppler Rate Bias Errors 
Beacon Doppler Rate Bias Error 
Baro Altimeter Correlated Error 
Radar Altimeter Correlated Error 

The upper -left 18 x 18 submatrix of P Fq is identical to that 
used in the System A evaluation. The last four rows and columns are all 
zeros except for the four diagonal elements. These are the squares of 
the initial rms values given or indicated in the 'Vr c M column in Table 4.3-1. 

TABLE 4. 3-1 

FILTER STATE MEASUREMENT ERROR STATISTICS: SYSTEM C 


State Numbers 

1-3 
4-6 
7-9 
10-12 
13-15 
16-18 . 

19 

\ 20 
21 
22 


Measurement 


1. One-Way Doppler 

Onboard Rate Bias 
Beacon Rate Bias 

2 . Baro Altimeter 

3 . Radar Altimeter 


Correlated Error 
(Filter State) 



0.0876 ft/sec 
0.0044 ft/sec 

Eq. (4.1-10) 

Eq. (4.1-14) 


800 sec 
400 sec 

40 sec 

4 sec 


Measurement 

Noise 


0.0357 ft 

Eq. (4.1-11) 
3.28 ft 
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The transition matrix elements are programmed in accordance with the 
description given in Ref. 4 (see pages 46 and 53) and discussed in Section 4. 1. 

The filter covariance propagation and update calculations are 
those given in Eqs. (4. 1-3), (4. 1-4) and (4. 1-5) except for one additional 
term used to account for measurement nonlinearities. The denominator 
of the gain equation (4. 1-4) is increased by an amount AR as follows: 

K F - Vf /(h f p f h f “70 + V "• AR > < 4 - 3 - a) 

The increase is calculated using the second partial derivatives of a non- 
linear measurement equation which is expanded in a Taylor series about 
the estimated state. The second partial derivatives with respect to position 
and velocity make up a 6 x 6 matrix, S p . Then, a 6 x 6 matrix product 
is calculated 

M p = Sp Pp (4.3-2) 

where P p is the 6x6 upper -left partition of P p . M p is used to compute 
a double -sum term, DSUM, and a quadruple -sum term, QSUM. 

DSUM = Trace (Mp) (4.3-3) 


where 


QSUM 




c..M p (i, ]) M p (j, i) 


(4. 3. 4) 



1: i = j 
2: i tM 
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Then, 

AR = QSUM + ^ DSUM) 2 (4.3-5) 

It is shown in Ref. 10 that i DSUM is the expected value of the measurement 
bias error associated with the second-order terms of the Taylor series. 
The factor u ?0 in Eq. (4.3-1) is the measurement underweighting factor 
discussed in Section 4.1. 

The 1 x 22 matrices can be viewed as 22-dimensional row 
F 

vectors, h. The non-zero elements of h for each of the three scalar 
measurements are given below. 


One-Way Doppler 





At 

2 




h 


3 



P z X At/^z; 

p 3 ‘ 2 v° 



h 4 = AtP x /P 
h 5 = AtP y /P 
h 0 = A t P Z /P 
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Where the p k , p, d z , P x> etc. are I frame components of the displace- 
ment vector, o, from beacon to shuttle, and its derivatives £ and £. 
The other parameters appearing above are defined as follows: 

o = JE ■ E 

X = P * P 


7? =p*0 


At = 

update interval 

(4.3-7) 

Raro Altimeter 

Fnitar Altimeter 


mjL 

II 

X 

h^ = x/r 


JL 

h 2 = y/r 

4 = y/r 


h 3 = z/r 

h 3 = z/r 


h 21 = 1 

h 22 ’ 1 

(4.3-8) 


In the above r is the magnitude of the current position, from the center of 
the earth, and x, y and z are its I frame components. 

The filter covariance program based on the 22 -state filter out- 
lined above has been exercised over the trajectory and measurement 
schedule described in Section 2. The time history of the filter -indicated 
performance is given in Table 4. 3-2. RMS values of position and velocity 
errors in R frame components are given just before and after the first 
measurement of each type and at a few other times, including touchdown. 
The results given in Table 4.3-2 represent what optimal performance. 
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TABLE 4.3-2 

SYSTEM C FILTER-INDICATED PERFORMANCE 


Event 

(see Table 2.1-1) 

Time 

(sec) 

RMS Position Errors 
(R frame (It)) 

RMS Velocity Errors 
(R frame (It/sec)) 

Vertical 

Downrange 

Crossrange 

Vertical 

Downrange 

Crossrange 

t owd 

1686" 

1686 + 

49,400 

34.300 

93.100 

54.100 

25,000 

24.900 

160 

llo 

93.6 

49.6 

134 

ELS 

t ba 

1806' 

lBoe 4 - 

10,300 

4,910 

15,400 

9,610 

19,100 

11,200 

34.2 

17.0 

16.0 

10.6 

9.77 

9.02 


2202 

140 

63.4 

88.2 

0.87 

0.31 

0.85 

t ra 

2430' 

2430 + 

35.8 

33.7 

65.8 

65.7 

77.6 

77.1. 

0.32 

0.31 

0.11 

0.11 

0.84 

0.83 

t td 

2432 

2.1 

0.7 

3.6 

0.22 

0.07 

0.33 

OTHER INDH 

Misaligns 
Gyro Dril 
Aeceleron 
Acceleron 
Onboard T 
Ground Ra 
Baro Altli 
Radar Alt 

CATED ES" 

ents (rad) 
(rad/sec) 
leter Scale 
neter Biasi 
.ate Bias (I 
itc Bias (It 
neter Errc 
Imeter Err 

TIMATION ERRORS AT TOUCHDOWN, T TD 

- 0.266X10' 3 , 0.066*10' 3 , 0.056xl0- 3 ) 

— 0.256 x 10*3, 0.281x10-6, 0.276x10-° l. 

Factors - 4,91 xlO’ 3 , 4.48 xlO'5 4.84 x 10'5 \ 

>s (ft/sec 2 ) - 5.67 x 10-4, 5,60 x 10* 4 , 5.64 xlO" 4 ' 
t/sec) - 0.051 

/sec) — 0.007 

r (ft) - 9.53 

or (It) — 1.89 


would be if this 22-state filter model were a completely accurate repre- 
sentation of real-world error dynamics. The results given in Section 5. 3 
are a projection of the performance of this 22 state sub-optimal filter in 
a 63-state model of the real world 0 


4. 4 SYSTEM D FILTER 

The JSC/TRW System D filter has been represented as a 
23 -state filter with the following state definitions. 

The rms misalignment estimation errors are in runway coordinates 
(Vertical, Downrange, Crossrange); the gyro and accelerometer errors 
are given in I frame (X, Y, Z) coordinates — see Appendix A. 
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State Numbers 

Variables 

1-3 

Position Errors 

4-6 

Velocity Errors 

7-9 

Platform Misalignments 

10-12 

Gyro Drifts 

13-15 

Accelerometer Scale Factor Errors 

16-18 

Accelerometer Bias Errors 

19 

DME Correlated Error 

20 

VOR Correlated Error 

21 

Baro Altimeter Correlated Error 

22 

ILS Localizer Correlated Error 

23 

Radar Altimeter Correlated Error 


The upper left 18 x 18 submatrix of Ppg has been supplied by NASA 
(a different set from that used in the System A and C evaluations). The last 
five rows and columns are all zeros except for the five diagonal elements. 
These are the initial rms values given or indicated in the "cr c " column of 
Table 4.4-1. 


TABLE 4.4-1 

FILTER STATE MEASUREMENT ERROR STATISTICS: SYSTEM D 


Measurement 

Correlated Error (Filter State) 

Measurement 

Noise 

cr 

m 

a 

c 

T 

1. DME 

295 ft 

400 sec 

48 ft 

2. VOR 

18. 67 mrad 

800 sec 

8. 5 mrad 

3. Baro Altimeter 

Eq. (4.1-10) 

40 sec 

Eq. (4.1-11) 

4. ILS Localizer 

Eq. (4.1-13) 

0. 55 sec 

0 

5. Radar Altimeter 

Eq. (4.1-12) 

4 sec 

3.28 ft 
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The transition matrix elements are programmed in accordance 
with the description given in Ref. 4 and discussed in Section 4. 1. The 
filter covariance propagation and update calculations are those given in 
Eqs. (4. 1-3), (4. 1-4) and (4. 1-5) 


The 1 x 23 matrices can be viewed as 23 -dimensional row 
r 

vectors, h. The non -zero elements of h for the first two scalar measure- 
ments are given below: 



DME 

VOR 


h l 

= V p i 

h l = u 2 x/ 

% 

h 2 

= °y/ p i 

h * = u v 

,P 1 

h 3 

= 9 %/°i 

h3 = \/ 

/Ps l 

h 19 

= l 

h 20“ 1 



where is the magnitude of the vector displacement p ^ of the vehicle 

from the V OR/D ME station, p , p and p are is I frame components, 

x y z 

p s is the magnitude of the projection of p- onto the earth’s surface, and 
U 2 X , ug and ug z are the I frame components of the unit vector Ug shown 
in Fig. 3.4-2. 


The non-zero elements of the h vectors associated with baro 
altimeter , radar altimeter, and ILS localizer measurements are the same 
as in the System A filter; see Eqs. (4. 1-7) and (4. 1-8). 


The filter covariance program based on the 23 -state filter out- 
lined above has been exercised over the trajectory and measurement sched- 
ule described in Section 2. The time history of the filter -indicated per- 
formance is given in Table 4. 4-2. RMS values of position and velocity 
errors in R frame components (runway coordinates; see Appendix A) are 
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TABLE 4.4-2 

SYSTEM D FILTER-INDICATED PERFORMANCE 



Time 

RMS Position Errors 
(R frame (ft)) 

RMS Velocity Errors 
(R frame (ft/sec)) 



Vertical 

Downrange 

Crossrange 

Vertical 

Downrange 

Crossrange 

T DME 

m 

4,600 

4,580 

27. 200 

14.200 

22,600 

12,800 

22.8 

17.6 

36.9 

35.5 

16.6 

15.9 

>-3 

to 

> 

1806" 

wmm 

■PH 

wmm 



12.4 

180G + 



m 

S9- 

BmHI 

11.8 

Tt nr 

2292" 

780 

290 

580 

2,82 

2.52 

3.55 

LiCLt 

2292 + 

773 

285 

421 

2.78 

2.38 

2.81 


ESI 


284 

12.3 

KB 


m 

t ra 

■ ' 


270 

11.6 

MEM 

KB 1 

H M 

t td 

2482 

2.5 

54.2 

5.7 

gg 

1.33 

0.30 

t td 

(old trajectory) 

2482 

2.3 

43.4 

5.6 

0.31 

1.29 

0.29 


OTHER INDICATED ESTIMATION ERRORS AT TOUCHDOWN, T TD 


Misalignments (rad) 

Gyro Drift (rad/sec) 
Accelerometer Scale Factors 
Accelerometer Biases (ft/ sec®} 
DME Error (ft) 

VOR Error (raid) 

Baro Altimeter Error (It) 
Localizer Error (rad) 

Radar Altimeter Error (ft) 


ixlO" 3 , 0.122xl0" 3 , 0.211 xlO -3 
1x10"®, 0.546x10"®, 0,526x10"® 
x 10"®, 9.21 xlO" 5 , 9.91 x 10"® 


0.545 
0.543 

9.97 - - , - — . ... , — , 

1.59 xlO" 3 , 1.56 xlO" 3 , 1.59 xlO" 3 
2.48 
2.38 
8.61 
0.54 
2.04 


-3 


xlO 
xlO" 3 


given at key times — just before and just after the first measurement of each 
type, and at touchdown. At the bottom of the table the filter -indicated rms es- 
timation errors, at touchdown, for the other 17 filter states are also given. 

The added row of position and velocity errors at touchdown, 
labeled ”T T ^(old trajectory)”, corresponds to an additional run through 
the System D filter covariance program with the trajectory tape which had 
been used in the System A, B and C evaluations. This extra run was 
motivated by the observation that the System D filter -indicated touchdown 

performance (with the new trajectory) is very similar to the System A 

— 

The rms misalignment estimation errors are in runway coordinates 
(Vertical, Downrange, Crossrange); the gyro and accelerometer errors 
are given in I frame (X, Y, Z) coordinates — see Appendix A. 
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filter -indicated performance (Table 4.1-2) in every component except 
vertical velocity — 0,50 ft/sec compared with the System A value of 
0.16 ft/sec. The added run, in which the terminal rms vertical velocity 
error is 0.31 ft/sec, shows that a major part of the difference between 
the System A and System D results is caused by the steeper approach 
path for System D. This point is discussed further in Section 6.1. 

The results given in Table 4.4-2 represent what optimal per- 
formance would be if this 23-state filter model were a completely accu- 
rate representation of real-world error dynamics. The results given in 
Section 5.4 are a projection of the performance of this 23 -state sub- 
optimal filter in a 74-state model of the real world. 
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5. 


RESULTS 


This section presents detailed results for the four candidate 
systems . Overall performance curves are given for each system, show- 
ing time histories of position and velocity errors due to all sources com- 
bined. In each case a "Baseline Error Budget" table is given, showing 
the contributions of individual error sources, or small groups of error 
sources, at touchdown. Smaller tables of "Alternative Contributions" 
are provided, allowing the reader to see the effects of different truth 
model assumptions. Sensitivity curves illustrating the effect of variations 
in the rms value of a given major error source (or group of error sources) 
on overall system performance are also given , Section 6 provides a 
general su mmar y and comparison of the four sets of results, including 
discussions of how particular groups of error sources contribute dif- 
ferently to different systems and of filter states which might be safely 
deleted or usefully added. 


5.1 SYSTEM A EVALUATION 

The results presented below correspond to the trajectory seg- 
ment, lasting 796 seconds, pictured in Fig. 2.1-1. The measurement 
schedule is described in Section 2.1; it employs a 2 second interval be- 
tween updates for 750 seconds, followed by a 0.5 second update inter /al 
for the final 46 seconds. A file of gain values, generated using the filter 
covariance program outlined in Section 4.1, has been used repeatedly to 
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compute detailed error contributions corresponding to this trajectory and 
measurement schedule. 


5.1.1 Overall System Performance 

Figures 5.1-1 and 5.1-2 present overall performance curves 
for System A, showing position error components and velocity error com- 
ponents, respectively. The curves plotted represent rms errors due to 
the combined effects of all error source groups in the baseline error bud- 
get. They were generated by root-sum-squaring individual contributions 
at one minu te intervals (tabulated in Appendix C). At T^^g, when the first 
DME measurements are made, and at other key times, such as T BA , 

T _ T and T~ t a , the root -sum-square calculation was performed 
LC Xj R-A CjXjfe » « I 

both before and after update. Thus, the large jumps in certain component 

errors, which occur at these times, are accurately shown. Otherwise, 
the curves are faired -in over the one-minute intervals between the cal- 
culated points. 

At the time of the first DME measurements, T DME , some of 
the component errors increase — notably, the crossrange position and 
velocity errors. This reflects the fact that the initial filter covariance 
matrix does not accurately represent the assumed real world covariance, 
resulting in non-optimal performance. Clearly, if the filter were 
''tuned" for the particular case considered, it could do better at the first 
update time. Navigation errors then decrease rapidly over the several 
minutes following the first update. Note, especially, the large decrease 
in vertical position and velocity errors at the time of the first baro- 
altimeter update. The exception is the increase in crossrange velocity 
error during the first minute. Examination of individual contributions 
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(as tabulated in Appendix C) shows that the major cause of this increase 
is the group of accelerometer scale factor errors (Group 3). This group 
is also the major cause of the increase at the first update. The effect is 
associated more with initial correlations caused by the accelerometer 
scale factor errors during boost than with the accelerometer errors 
occurring in the time interval plotted. 

An interesting feature of the velocity performance curves of 
Fig. 5. 1-2 is the sharp decrease in the downrange velocity error at the 
time of the first localizer measurement, T LCL - Clearly, the crossrange 
and downrange position and velocity errors must be highly correlated at 
this time,, We note in passing that a filter scheme which automatically re- 
initialized to a diagonal filter covariance matrix at this point in time would 
lose the benefit illustrated. 

A potentially serious problem for System A is illustrated in 
Figs. 5.1-1 and 5.1-2 by the large vertical errors just prior to the time 
of the first radar altimeter measurements, T RA - This occurs less than 
one mi nute prior to touchdown, and the projected rms vertical position 
error of more than 2000 feet is probably unacceptable -- despite the fact 
that the performance is quite good after the radar altimeter comes into 
use. The large rms vertical position error prior to T RA is caused by the 
large vertical velocity error integrated over time. The major contributors 
to vertical velocity errors are the gyro bias and mass unbalance drifts 
(Groups 7 and 8). Evidently, the filter "believes” it has accurate esti- 
mates of platform alignment and vertical velocity information, and is 
making heavy use of mis-resolved accelerometer outputs. The problem 
illustrated here is not inherent to the type of conventional navaids under 
investigation (since it does not occur in the System B or D studies). It is 
a filter design problem and is, thus, subject to fixing by software modification. 
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Section 6 contains further discussions of the implications of 
the performance curves presented in this section. 


5,1,2 Detailed Error Budget 

System error covariances have been calculated for each group 
of error sources over the entire 796 second trajectory (t = 1686 sec to 
2482 see). The total amount of data produced is therefore quite large, 
since rms values of all significant quantities, before and after each 
measurement, have been generated. In order to summarize important 
results in one table of manageable size, a "snapshot" view of eonditionr 
at the nominal touchdown time (2482 sec) is presented in Tables 5. 1-1 
and 5. 1-2. Table 5. 1-1 is the System A Baseline Error Budget showing 
rms estimation errors in position, velocity and platform alignment com- 
ponents at touchdown. Each value is the rms contribution oi the error 
source or sources indicated in the left hand column; it comes from a 
computer run in which the elements of Pg(0), Qg or Rg, corresponding to 
those errors alone, are set at truth-model values • — and all others are 
zero. The overall System A performance, or root -sum- square values, 
are given at the bottom of the table and compared, first with the filter's 
own estimate of performance based on its simplified 24- state model of 
system error dynamics, and second with the Shuttle landing navigation 
specification given in Ref, 1, Table 5,1-2 lists the System A Alternative 
Contributions showing how different assumptions or additional error 
sources would produce other contributions to position, velocity and 
alignment errors at touchdown. 

Particular numbers have been circled in the position and 
velocity error columns of Table 5. 1-1 to focus attention on major 
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TABLE 5ol-l 

SYSTEM A BASELINE ERROR BUDGET 


Error Source 


ESTIMATED STATES AND NOISES 
Group 1. IMU Quantization Noise 
Measurement Noises 
DME's 

Radar Altimeter 
Radar Altimeter Markov 

BIU-RELATED STATES 
Group 2. Accelerometer Biases 

Group 3. Accelerometer Scale 
Factor Errors 

Group 4. Accelerometer Misalignments 
Group 5. Accelerometer Nontlnearltles 

Group 6. Gravity Anomalies and 
Deflections 

Group 7. Gyro Bias Drifts 
Group u. Gyro Mass Unbalances 
Group 9. Gyro Anlsoelasticlty 

EXTERNAL-AID RELATED STATES 


Value 


One Pulse - 1 cm/scc 

24 It 
3.28 ft 
Eq. (3.1-43) 

60 pg 

34 ppm 

40 sec 
3.5pg/g 2 
see Table 3.1-2 

0.01 deg/hr 
0.015 dog/hr/g 
0.005 deg/hr/g 2 


Total Projected Performance (root-sum-square) 


Filter-Indicated Performance 
Specification (Ref. 6) 


RMS Navigation Errors at Touchdown 


Position (ft) 


Vertical Dcwnrangc Crosarangc 


0.1 4.2 

CajD 4.o 


CgTp CUD 

0.3 

0.2 3.0 


0.2 

rprr 


0.7 


0.1 


7.8J 

3.6 

0.1 


o.e 

0.3 

0.0 

0.0 


Group 10, DME Bias Errors 

295 It 

0.2 

1.7 

0.2 

0.01 


0.02 

0.010 

0.004 

0.004 

G oup 11. DME Scale Factor Errors 
Group 12. Baro Altimeter Errors 

100 ppm 

— 

0.8 

0.1 


0.04 

0.01 

0.005 

0.002 

0.007 

Bias and 
Scale Factor 

100 ft \ 
0.03 J 

— 

0.1 

— 

— 

0.01 

• — 

— 

0.001 

0.002 

Markov and 

20 rt\ 










Measurement NoIsb 

5 ill „ 










Static Defect 

1.52 k 10*4 ft/(ft/sec) z 

— 

0.1 

— 

— 

0.01 

— 

— 

— 

0.001 

GGroup 13. ILS Bias Errors 











Localizer 

0.5 mrad 


1.2 

<XP 

— 

0.07 

C0T2O 

0.050 

0.025 

0.014 

Glide Slope 

0.5 mrad 

CoTD 

CUD 

0.1 

0.13 

0.07 

0.01 

0.003 

0.C01 

0.004 

Group 14. ILS 2 Hd Order Markovs 











Localizer 

Eq. (3.1-44) 

0.1 

2.4 

CITO 

— 

0.09 

CTScD 

0.049 

0.038 

0.017 

Glide Slope 

Eq. (3.1-45) 

CoTQ 


0.1 

0.01 

COD 

— 

0.007 

0.001 

0.016 

Group 15. (Net Included In Baseline 











Error Budget) 











Group 16. DME Survey Errors 
Group 17. ILS Survey Errors 

1 ft 

“““ 

0.6 

■““ “ 


0.01 


0.001 

— 

0.001 

Localizer 

1 ft 

— - 

— 

1.0 

— 

— 

— 

0.001 

— 

— 

Glide Slope 

1 ft 

0.2 

CYD 

— 

0.03 

0.11 

— 

0.005 

— 

0.011 


1.9 


20.7 


8.0 


1.9 


30.6 


5.5 


3.0 


80.0 


4.7 


Velocity (ft/soc) 


Vertical [ Dircnmriic [crojarargrj 



0.22 


1.07 


0.49 


0.16 


0.90 


0.34 


0.20 


3.00 


2.00 


Misalignment (mrad) 
Vertical j Docrnransc |CfQ3arancc 


0.048 0.032 0.034 


0.008 

0.004 

0.000 

0.541 

1.608 


0.002 


0.075 

0.114 


0.739 

0.2GB 

0.005 


U.043 

0.012 

0.001 


1.660 0.232 


0.267 0.079 


0.006 

0.00B 

0.012 

0.078 

0.013 


0.178 0.172 0.201 


0.025 0.021 0.021 


0.018 

0.03G 

0.001 


0.226 


0.123 


Nat Specified 


A " — " entry refers to a negligibly small contribution/ 
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TABLE 5.1-2 

SYSTEM A ALTERNATIVE CONTRIBUTIONS 

! RMS Navi uni ion Errors at Touchdown 


Error Source 


Group 1. Initial 9*9 '"^variance 

Matrix (Position, Velocity 
and Alignment) 

Initial 6*6 Covariance 
Matrix (Position and Velocity 
only) 

ILS First-Order Markovs 
Localizer 
Glide Slope 

Group 15. Radar Altimeter Errors 
Bias 

Scale Factor 


See Section 3. 1. 3 


Eq. (3.1-44) 
Eq. (3.1-45) 


Position (ft) 


Velocity (ft/sec) 


Misalignment (mrad) 


Vertical 

j Down ran go 

CroB3range 

Vertical 

Dawn range 

Crossrange 

Vertical 

Der m ranee 

CroFSrarce 

0.4 

19.9 

2.2 

0.06 

0.96 

0.30 

1.920 

0.126 

0.199 

0.4 

13.0 

2.1 

0.04 

0.63 

0.30 

1.827 

0.122 

0.121 


0.082 0.064 

0.007 0.001 


0.005 — 

0.005 — 
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contributors. The rule followed in deciding which numbers to circle was 
to include, in any one cblumn, the largest value, second largest value, 
etc. , until 95% of the RSS total value was reached. The following 
general, statements may be drawn from an examination of the System A 
baseline error budget: 


© The projected downrange performance, in both 
position and velocity, is well within specs, and 
the filter-indicated performance is within 20% 
of projected actual performance. The largest 
contributors to do wnrange errors are gyro bias 
drifts (Group 7). Other significant error sources 
are accelerometer errors (Groups 2, 3 and 4), 
radar altimeter markov error (Group 1), glide 
slope errors (Groups 13,14 and 17), IMU quanti- 
zation noise (Group 1), and DME bias errors 
(Group 10). 

« The projected vertical performance is within spec 
;.n position, but out-of -spec in velocity. The 
filter -indicated performance is the same as the 
projected actual performance in position, but 
almost 30% off in velocity — erring on the opti- 
mistic side. The largest contributors to vertical 
velocity errors are glide slope bias error (Group 13) 
and the accelerometer mi salignments (Group 4). 
Other significant contributors to vertical errors are 
radar altimeter errors (Group 1), glide slope 
markov errors (Group 14), gyro biases and mass 
unbalances (Groups 7 and 8), accelerometer bias 
and scale factor errors (Groups 2 and 3) and IMU 
quantization (Group 1). 

c The projected crossrange performance is within 
spec in velocity, but out-of-spec in position. The 
filter -indicated performance , in both position and 
velocity, is optimistic — predicting roughly two 
thirds of the projected actual rms errors. The 
largest contributors to crossrange errors are 
localizer errors (Groups 13 and 14) . Other signi- 
ficant contributors to crossrange' velocity errors 
are accelerometer scale factor errors (Group 3). 


5-8 


THE ANALYTIC SCIENCES CORPORATION 


* The projected errors in estimating platform misalign- 
ments are approximately 0.2 milliradians, for mis- 
atignments about level axes, and nearly 2.0 millira- 
dians for azimuth misalignment. The latter is signifi- 
cantly larger than the projected rms actual misalign- 
ment which is 0.620 milliradians — -the rss value at 
the bottom of the System B error budget iable. (.Mis- 
alignments are not estimated by the System B filter). 

Considering all components of position and velocity errors at 
touchdown, the list of major contributors includes gy ro errors, accelero- 
meter errors, and error s associated with all of the extern al devices ex- 
cept the baro-altimeter . More detailed discussion of the important error 
mechanisms and possible ways to reduce some of the contributions is 
given in Section 6 . 

Detailed tabulations of the error contribution time histories 
are given in Appendix C. For every row in Tables 5. 1-1 and 5. 1-2 there 
is a page in Appendix C, which is a reproduction of a computer print-out 
page summarizing important results from one error budget run. 


5.1.3 Sensitivity Curves: System A 

This section contains several curves illustrating the sensitivity 
of System A performance to variations in error source statistics. These 
"fixed-filter" sensitivity calculations answer the question: "What is the 
effect of an unknown variation in the rms value or values or an error 
source or group of error sources?" These calculations can be made 
easily, given the type of error budget information summarized in Table 
5. 1-1, because the appropriate error covariance equations are linear. 

All of the example curves given in this section correspond to 
major contributors to the system performance. Similar sensitivity curves 
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may be constructed for any error source group for which error budget 
data exists in Table 5 0 1-1. Sensitivity curves corresponding to error 
sources which produce minor contributions when their nominal values are 
assumed are quite flat and of little interest. 

To illustrate the means by which the data points for the sensi- 
tivity curves were calculated an example will be given. The sensitivity 
of downrange position at touchdown to gyro bias drift is shown in Fig. 5.1 -3a. 
The baseline data point, 28.7 ft, is the total system error, which includes 
the effect of a 0. 01 deg/hr bias drift about each platform axis. The con- 
tribution of these bias drifts is shown in Table 5.1-1 to be 17. 8 ft. To 
compute the effect of a 0.03 deg/hr gyro drift the 0.01 deg/hr bias drift 
contribution is removed from the total system error and replaced with an 
error which is three times as large. Thus: 


D 



-(17.8) 2 + (53.4) 2 



= 57.9 ft 


This result is indicated in Fig. 5.1-3 as a boxed -in point. Also shown on 
the curve is the specification value for downrange position error (80 ft) 
and the corresponding gyro drift which would cause it (0.043 deg/hr) can 
be found from the curve. The dashed line through the origin is the asym- 
tote approached by the total system error curve as the gyro bias drift con- 
tribution becomes large and dominates all others. The remaining points 
on the curve are obtained as illustrated above. 

Figure 5. 1-4 shows the sensitivity of vertical position and 
velocity errors to variations in the accelerometer misalignment angles. 

It is evident that the vertical velocity error won’t meet the spec (0. 2 ft/ sec) 
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time of touchdown 


8 o SPECIFICATION 11MIT A f. 

A/ 

// 

// 

6 ° - 

As~' VALUE 
AS COMPUTED 
Ay IN THE 

Ay EXAMPLE 

40 - BASELINE A^ / 

VALUE / 

V/ s 
JA / 

/ 

20 - f 

/\ ■ X x I 

o 0.01 0.02 0.03 0.04 0.05 

RMS GYRO BIAS DRIFT (deg/hr) 


(a) Position 


TIME OF TOUCHDOWN 


SPECIFICATION LIMIT 


BASELINE 
VALUE A/ 

\ // 

// 


0.01 0.02 0.03 0.04 0,05 


RMS GYRO BIAS DRIFT (deg /hr) 

(b) Velocity- 


Figure 5,1-3 Sensitivity of Downrange Position and 
Velocity "’rrors at Touchdown to 
Gyro Bias '‘rift: System A 




RMS VERTICAL VELOCITY ERROR {ft/sec) RMS VERTICAL POSITION ERROR (ft) 
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Figure 5.1-4 Sensitivity of Vertical Position and 
Velocity Errors at Touchdown to 
Accelerometer Misalignment: System A 
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even if the accelerometer misalignment is very small. That is, the other 
system error sources cause the system to just barely meet spec even 
with no accelerometer misalignment. 

The data of Fig, 5,1-5 shows the sensitivity of vertical posi- 
tion error to radar altimeter markov error. The radar altimeter markov 
error standard deviation is a function of altitude. Thus, the ratio of the 
assumed function to the baseline function is used as the independent vari- 
able in this figure. 



RELATIVE RMS RADAR ALTIMETER MARKOV ERROR 


Figure 5.1-5 Sensitivity of Vertical Position Error at Touchdown 
to Radar Altimeter Markov Error: System A 

Figures 5.1-6, 5.1-7 and 5,1-8 plot the sensitivity of cross - 
range errors to their three largest contributors. Figure 5.1-6 is a plot 
of crossrange position error as a function of localizer markov error. As 
in Fig. 5.1-5, the independent variable of Fig. 5.1-6 is the ratio of the 
assumed to the baseline function used to generate the data of Table 5.1-1. 
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Figure 5=1 


R-9676 



-6 Sensitivity of Crossrange Position Erroi at Touchdown 
to Localizer Markov Error : System A 


R-9675 



Figure 5=1-7 


Sensitivity of Crossrange Position Error 
at Touchdown to Localizer Bias Error: 
System A 





RMS CROSSRANGE POSITION ERROR (ft] 
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(a) Position 



Figure 5.1-8 Sensitivity of Crossrange Position and Velocity 
Errors at Touchdown to Accelerometer Scale 
Factor Errors: System A 
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This is a result of the fact that the localizer markov error standard de- 
viation is a function of the slant range to the ILS localizer antenna. Note 
that the spec is exceeded for all values of this error because of the other 
error source contributions. The crossrange position error is also exceeded 
for all values of localizer bias error as shown in Fig. 5.1-7. The final 
figure in this section, Fig. S.l^B, shows the sensitivity of crossrange 
position and velocity to accelerometer scale factor error. Again, there 
is no value of accelerometer scale factor error for which the crossrange 
position error specification is met. 


5.2 SYSTEM B EVALUATION 

The results presented below correspond to the trajectory and 
measurement schedule described in Section 2 -- exactly the same as that 
used in the System A evaluation. A file of gain values, generated using the 
6-state, square root filter covariance program described in Section 5. 2, 
has been used repeatedly to compute detailed error contributions for this 
trajectory and measurement schedule. 


5.2.1 Overall System Performance 

Figures 5.2-1 and 5.2-2 present overall performance curves 
for System B, showing position error components and velocity error com- 
ponents, respectively. The curves plotted represent rms errors due to 
the combined effects of all error source groups in the baseline error bud- 
get. They were generated by root-sum-squaring individual contributions 
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at one minute intervals. (This data is tabulated in Appendix C = ) At . 

when the first DME measurements are made, and at other key times, such 

a a rp T m an( j the root-sum -square calculation was per- 

as a bA’ LCL’ RA GL3 

formed both before and after update. Thus, the large jumps in certain 
component errors, which occur at these times, are accurately shown. 
Otherwise, the curves are faired-in over the one-minute intervals between 
the calculated points. 

At the time of the first DME measurements, T DME , all of the 
position error components decrease. (Since a diagonal initial filter co- 
variance matrix was used, no ’’wrong way” corrections can occur. ) 
Following the first update the position errors generally increase over the 
next one or two minutes. All are reduced considerably at the time of the 
first baro-altimeter measurement, T BA - Two of the velocity component 
errors increase significantly over the first three minutes. Examination 
of the individual contributions (as tabulated in Appendix C ) shows that the 
major cause of these increases in position and velocity errors are the 
unmodeled DME bias and scale factor errors. The ’’humps” in the 
velocity component error curves near t = 8 minutes and t = 10 minutes 
are also largely due to the unmodeled DME errors, especially scale fac- 
tor errors. 

All three position error components for System B are reasonably 
small and well-behaved after the baro-altimeter comes into use, at T BA - 
While the errors are generally larger than the corresponding System A 
errors, the problem of the large vertical errors prior to T RA (exhibited 
by System A — see Section 3=1-1) is avoided= 

Section 6 contains further discussions of the implications of the 
performance curves presented in this section. 
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5.2.2 Detailed Error Budget 

System error covariances have been calculated for each group 
of error sources over the entire 796 second trajectory. Table 5. 2-1 is 
the System B Baseline Error Budget showing rms estimation errors in 
position, velocity and platform alignment components at touchdown. Each 
value is the rms contribution of the error source or sources indicat ed in 
the left hand column; it comes from a computer run in wh ich the elements 
of those errors alone are set at truth-model values -- and all othe rs are 
zero. The overall System B performance, or root -sum- square values, 
are given at the bottom of the table and compared, first with the filter’s 
own estimate of performance based on its simplified 6-state model of 
system error dynamics, and second with the Shuttle landing navigation 
specification given in Ref. 1. Table 5.2-2 lists the System B Alternative 
Contributions showing how different assumptions or additional error 
sources would produce other contributions to position, velocity and align- 
ment errors at touchdown. 

Particular numbers have been circled in the position and velocity 
error columns of Table 5.2-1 to focus attention on major contributors. The 
rule followed in deciding which numbers to circle was to include, in any 
one column, the largest value, second largest value, etc., until 95% of 
the RSS total value was reached. The following general statements may 
be drawn from an examination of the System B baseline error budget. 

• The projected downrange performance , in both posi- 
tion and velocity, is slightly out-of-spec, and the 
filter -indicated performance is over -optimistic by a 
wide margin. The largest contributors to downrange 
errors are DME bias errors (Group 10). Other signi- 
ficant error sources are gyro bias drifts (Group 7 ), 
baro-altimeter bias and scale factor errors (Group 12), 
radar altimeter markov error (Group 15) and glide 
slope bias error (Group 13). 
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TABLE 5.2-1 

SYSTEM B BASELINE ERROR BUDGET 




RMS Navigation Errors at Touchdown 

Error Source 

Value 

Position (ft) 

Velocity (ft/sec) 

Misalignment (mrad) 



Vertical 

| Do* orange 

| Cro« range 

Vertical 

Dow orange 

Crossrangc 

Vertical j Dffvrorango j Crao* range 

1. UNCORRELATED NOISE 











Group 1. IMU Quantization Noise 

One Pulse - cm/sec 

0.2 

2.7 

0.4 

0.05 

0.14 

0.08 

0 

0 

0 

Measurement Noises 











DME's 

24 ft 

— 

0.7 

— 

— 

0.02 

— 

0 

0 

0 

Radar Altimeter 

3.28 ft 

cud 

6.9 

0.1 

COD 

0.10 

0.01 

0 

0 

0 

n. IMU -RE LATED STATES 











Group 2. Accelerometer Biases 

60 pg 

0.1 

1.6 

0.2 

0.02 

0.05 

0.02 

0.027 

0.060 

0.054 

Group 3. Accelerometer Scale 

34 ppm 

— 

0.7 

0.1 

0.01 

0.03 

0.01 

— 

— 

— 

Factor Errors 











Group 4. Accelerometer Misalignments 

40 sec 

0.2 

4.1 

0.5 

0.04 

0.10 

0.08 

0.088 

0.192 

0.173 

Group 5. Accelerometer Nonlincarities 

3.5 pg/g2 

— 

0.1 

— 

— 

0.01 

— 

0 

0 

0 

Group 6. Gravity Anomalies and 

see Table 3. 1-2 

0.1 

1.2 

0.2 

0.01 

0.05 

0.02 

0 

0 

0 

Collections 











Group 7. Gyro Dias Drifts 

0.01 deg/hr 

0.1 

dED 

CHD 

0.04 

CUD 

CHID 

0.604 

0.372 

0.637 

Group 8. Gyro Mass Unbalances 

O.OI5deg/hr/g 

— 

1.0 

0.2 

— 

0.10 

0.02 

0.105 

0.032 

0.057 

Group 9. Gyro Anisoefnsticity 

O.OOSdeg/hr/g 2 

— 

— 

— 

— 

— 

— 

0.001 

0.001 

0.001 

El. EXTERNAL-AID RELATED STATES 


C£D 









Group 10. DME Bias Errors 

295 ft 

1.2 

COO 

COD 

0.10 

0 

0 

0 

Group 11* DME Scale Factor Errors 

100 ppm 

CoTQ 

6.0 

CHD 

COD 

0.46 

CoTesD 

0 

0 

0 

Group 12. Baro Altimeter Errors 


C~3G.(L? 







Bias and 
Scale Factor 

0.03/ 

CoX> 

1.0 

0.01 

CjToO 

0.07 

0 

0 

0 

Markov and 

20 ft) 


2.7 



0.07 





Measurement Noiso 

5 ft J 

1 1 


' '* 

"" " 

0 

0 

0 

Static Defect 

1.52xl0' 4 ft/(ft/sec)2 

0.1 

7.4 

0.1 

— 

0.20 

0.01 

0 

0 

0 

Group 13. ILS Dias Errors 











Localizer 

0.5 mrad 

— 

0.B 

C HD 

— 

0.01 

0.18 

0 

0 

0 

Glide Slope 

0.5 mrad 

<£E> 

oro 

— 

0.03 

0.16 

0.01 

0 

0 

0 

Group 14. ILS 2 nd Order Markovs 











Localizer 

Eq. (3.1-44) 

0*1 

1.2 

Cs.O 

0.01 

O.Oz 

CeUD 

0 

0 

0 

Glide Slope 

i Eq, (3.1-45) 

0.2 

11.9 

0.2 

0.01 

0.02 

U.01 

0 

0 

0 

Group 15. Radar Altimeter Markov 

Eq. (3.1-43) 

COD 

CHID 

0.2 

COD 

0.26 

0.01 

0 

0 

0 

Group 16. DME Survey Errors 
Group 17. ILS Survey Errors 

1 ft 

■■ ■ 

0.2 

■* 

- 

0.01 


0 

0 

0 

Localizer 

1 ft 

— 

0.1 

1.0 

— 

— 

— 

0 

0 

0 

Glide Slope 

1 ft 

— 

14.4 

0.1 

0.02 

0.13 

— 

0 

0 

0 

Total Projected Performance [root-sum-squa 

re) 

2.2 

BS.6 

9.1 

0.34 

3.60 

0.89 

0.620 

''.424 

0.664 

Filter-Indicated Performance 

1.4 

21.3 

4.1 

0.19' 

0.48 

0.26 

Not Estimated 

Specification (Ref. 6} 

3.0 

80.0 

4.7 

0.20 

3.00 

2.00 

Not Specified 


A " — " entry refers to a negligibly email contribution. 
A "O'' entry represents an Identically zero contribution. 
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• The projected vertical performance is within spec 
in position, but out -of -spec in velocity „ The filter 
indicated performance is somewhat optimistic — 

35% off in position, 45% off in velocity. The largest 
contributors to vertical errors are the radar alti- 
meter markov error (Group 15) and DME scale fac- 
tor errors (Group 11). Other significant error 
sources are DME bias errors (Grump 10), radar 
altimeter measurement noise (Group 1), baro- 
altimeter bias and scale factor errors (Group 17) 
and glide slope bias (Group 13). 

e The projected crossrange performance is within 
spec in velocity, but out-of-spec in position. The 
filter -indicated performance is over-optimistic by 
a wide margin. The largest contributors to cross - 
range errors are localizer errors (Groups 13 and 14). 
Other significant contributions are those of DME 
scale factor errors (Group 11) and gyro bias drifts 
(Group 7). 

• The projected platform misalignments (which are 
not estimated by the System B filter) are in the 
range from 0.4 to 0.7 milliradians . 


Considering all components of position and velocity errors at 
touchdown, the list of major contributors includes gyro errors, accel- 
erometer errors, and errors associated with all of the external devices 
including, in this case, the baro altimeter. More detailed discussion of 
the important error mechanisms and possible ways to reduce some of the 
contributions is given in Section 6 = 


Detailed tabulations of the error contribution time histories 
are given in Appendix C. For every row in Tables 5. 2-1 and 5. 2-2 there 
is a page in Appendix C, which is a reproduction of a computer print- 
out page summarizing important results from one error budget run. 
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5.2,3 Sensitivity Curves: System B 

This section contains several curves illustrating the sensitivity 
of System B performance to variations in error source statistics. These 
curves were produced using the error budget data of Table 5. 2-1. The me- 
thod of generating the data for these curves is discussed in Section 5.1.2. 

Similar sensitivity curves may be constructed for any error 
source group for which error budget data is available in Table 5.2-1. All 
of the error curves in this section correspond to major contributors to 
system performance. Sensitivity curves corresponding to error sources 
which produce minor contributions are quite flat. 

Figure 5.2-3 shows the sensitivity of downrange position and 
velocity errors to DME bias errors. For System B the downrange errors 
do not meet the spec for either position or velocity unless the DME bias 
is less than the baseline value. This inability to meet spec results from 
the fact that the filter does not calibrate the DME 1 bias error, as the 
Shuttle passes over it, as is done by the System A filter. A further dis- 
cussion of this problem may be found in Section 6. 

The downrange position and velocity error sensitivities to gyro 
bias drift are shown in Fig. 5.2-4. Note that the sensitivity of System B 
to gyro bias drift is greater than that of System A. ' is is also discussed 
in Section 6 . 

The position and velocity error sensitivity to radar altimeter 
markov error is shown in Fig. 5,2-5. The independent variable in this 
curve is the ratio of the assumed standard deviation of the markov error 
to the baseline value. This is done since the baseline standard deviation 
is a function of altitude as given by Eq. (3.1-43). 
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RAAS DOWNRANGE VELOCITY ERROR (ff/sec) RMS DOWNRANGE POSITION ERROR (ft) 


Figure 


R-968S 




.2-3 Sensitivity of Downrange Position and Velocity 
Errors at Touchdown toDME Bias Errors: 
System B 





RMS DOWNRANGE POSITION ERROR (ft) 


THE ANALYTIC SCIENCES CORPORATION 


Figure 5=2 



(a) Position 



-4 Sensitivity of Downrange Position and Velocity Errors 
at Touchdown to Gyro Bias Drift: System B 




) 
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RELATIVE RMS RADAR ALTIMETER MARKOV ERROR 


(a) Position 



RELATIVE RMS RADAR ALTIMETER MARKOV ERROR 


(b) Velocity 

Figure 5.2-5 Sensitivity of V ertical Position and V elocity 

Errors at Touchdown to Radar Altimeter 
Markov Error: System B 
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Figure 5.2-6 shows the sensitivity of vertical velocity to DME 
scale factor error. The vertical velocity does not meet the spec for any 
value of DME scale factor error „ 

Figure 5.2-7 shows crossrange position sensitivity as a func- 
tion of localizer markov error. The independent variable is again a ratio 
of assumed standard deviation to the baseline given by Eq. (3.1-44). It 
should be noted that the contributions to System A and to System B of this 
error source are similar. This is observed by noting the slopes of the 
straight lines in Figs. 5.1-6 and 5.2-7. 

The final figure of this section, Fig. 5.2-8, indicates the sen- 
sitivity of crossrange position errors to gyro bias drift rate. This sen- 
sitivity of System B to gyro bias error results from the fact that the 
System B filter does not try to estimate platform misalignment. 



RMS DME SCALE FACTOR ERROR (ppm) 


Figure 5.2-6 Sensitivity of Vertical V elocity Error 

at Touchdown to DME Scale Factor 
Error: System B 


5-27 




RELATIVE RMS LOCALIZER MARKOV ERROR 


Figure 5.2-7 Sensitivity of Crossrange Position Error 
at Touchdown to Localizer Markov 
Error: System B 



RMS GYRO BIAS DRIFT |de 0 /hr] 

Figure 5.2-8 Sensitivity of Crossrange Position Error 
at Touchdown to Gyro Bias Drift: 

System B 
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5.3 SYSTEM C EVALUATION 

The results presented below correspond to the trajectory and 
measurement schedule described in Section 2, employing a 2 second 
update interval for the first 720 seconds and a 0. 5 second interval for the 
final 76 seconds. A file of gain values, generated using the 22-state fil- 
ter described in Section 4.3, has been used repeatedly to compute detailed 
error contributions for this trajectory and measurement schedule. 


5.3.1 Overall System Performance 

Figures 5.3-1 and 5.3-2 present overall performance curves 
for System C, showing position error components and velocity error com- 
ponents, respectively. The curves plotted represent rms errors due to 
the combined effects of all error source groups in the baseline error bud- 
get. They were generated by root-sum-squaring individual contributions 
at one minute intervals. (These are tabulated in Appendix C. ) At Tq^jv 
when the first doppler measurement is made, and at T rA , when the first 
baro altimeter measurement is made, the root-sum-square calculation 
was performed both before and after update. Thus, the large jumps in 
certain component errors, which occur at these two times, are accurately 
shown. Otherwise, the curves are faired -in over the one -minute intervals 

between the calculated points. 

At the time of the first one-way doppler measurement, TqwD’ 
some of the position and velocity component errors increase. This re- 
flects the fact that the initial filter covariance matrix does not accurately 
represent the assumed real-world covariance — just as in the System A 
evaluation (although the same components do not increase in both cases). 
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Errors then decrease rapidly over the three or four minutes following 
the first update. Note, especially, the large decreases in vertical posi- 
tion and velocity errors at T BA , the time of the first baro-altimeter mea- 
surement. Figures 5. 3-1 and 5.3-2 also show pronounced "humps” in 
certain component errors near t = 5 minutes and t = 8 minutes. The 
former corresponds to the interval, near Mach 1, during which baro- 
altimeter measurements are inhibited. The latter corresponds to a rather 
sharp turning maneuver. The largest single contributor during this turn- 
ing maneuver is the error in antenna motion compensation. However, 
other contributions increase significantly at this time also. 

Figure 5.3-3 re-plots the position error components over the 
final two minutes with an expanded scale. The most striking feature of 
these curves is the sharp reduction in the downrange error just after the 
t = 13 minutes point. This is the time at which the vehicle overflies the 
doppler beacon. The filter evidently is able to remove the remaining 
"constant of integration 1 :’ error (see discussion by Lear in Ref. 10) when 
this overflight occurs. It is likely that the reduction would occur at an 
earlier time if the beacon were moved farther from the end of the run- 
way. 

Section 6 contains further discussion of the implications of the 
performance curves presented in this section. 


5.3.2 Detailed Error Budget 

The System C baseline error budget is given in Table 5.3-1, 
showing rms estimation errors in position, velocity and platform align- 
ment components at touchdown. Each value is the rms contribution of_ 
the error source or sources indicated in the left hand column, it comes_ 
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Figure 5.3-3 System C Final Approach Performance 


from a computer run in which the elements of those errors alone are set 
at truth-model values — and all others are zero. The overall System C 
performance, or root -sum-square values, are given at the bottom of the 
table and compared, first with the filter’s own estimate of performance 
based on its simplified 22 -state model of system error dynamics, and 
second, with the Shuttle landing navigation specification given in Ref. 1 . 
Particular numbers have been circled in the position and velocity error 
colu mn s of Table 5.3-1 to focus attention on major contributors. The 
rule followed in deciding which numbers to circle was to include, in any 
one column, the largest value, second largest value, etc., until 95% of 
the RSS total value was reached. The following general statements may 
be drawn from an examination of the System C baseline error budget: 
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TABLE 5,3-1 


SYSTEM C BASELINE ERROR BUDGET 
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® The projected downrange performance, in both 
position and velocity, is far more accurate than 
called for in the preliminary specifications. 

The filter-indicated performance is even better. 
Significant contributors are doppler rate bias 
error (Group 1), accelerometer biases (Group 2), 
antenna motion compensation errors (Group 10), 
doppler survey errors (Group 16) and IMU 
quantization noise (Group 1). 

® The projected vertical performance i s within 
spec in position, but out-of-spec in velocity. 

The filter-indicated performance is within 20% 
of the projected actual performance in position, 
but almost 50% off in velocity -- erring on the 
optimistic side. The largest contributors to 
vertical velocity errors are the accelerometer 
misalignments (Group 4) — contributing 
0.26 ft/sec compared to a total of 0.39 ft/sec, 
against a spec of 0.20 ft/sec. Other significant 
contributor, s to vertical errors are radar alti- 
meter errors (Group 1), doppler onboard rate 
bias error (Group 1), IMU quantization noise 
(Group 1), accelerometer biases (Group 2), 
gravity anomalies (Group 6) and antenna motion 
compensation errors (Group 10). 

e The projected crossrange performance is within 
spec in velocity but out -of -spec in position. The 
filter -indicated performance, in both position 
and velocity, is optimistic — predicting approxi- 
mately one third of the projected actual rms 
errors. The largest contributors to crossrange 
errors are accelerometer scale factor errors 
(Group 3) and antenna motion compensation errors 
(Group 10). Other significant contributors are 
IMU quantization noise (Group 1), accelerometer 
misalignments (Group 4) and gyro bias drifts 
(Group 7). 
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® The projected errors in estimating platform mis- 
alignments are approximately 0.2 milliradians 
for misalignments about level axes, and 2.3 milli- 
radians for azimuth misalignment. The filter- 
indicated performance is considerably better for 
all three axes. The projected rms azimuth 
estimation error is larger than the actual rms 
misalignment (0.62 mrad) — shown as the rss 
vertical misalignment at the bottom of the 
System B error budget. 

Considering all components of position and velocity errors at 
touchdown, the list of major contributors includes gyro errors, accel- 
erometer errors, and errors associated with all of the external devices 
except the baro altimeter. The overall performance projected for Sys- 
tem C is similar to that of System A. Significant improvement is probably 
possible in crossrange position and vertical velocity performance, as dis- 
cussed in Section 6. 

Table 5=3-2 lists the System C Alternative Contributions show- 
ing how different assumptions or additional error sources would produce 
other contributions to position, velocity and alignment errors at touch- 
down. Detailed tabulations of the error contribution time histories are 
given in Appendix C„ For every row in Tables 5.3-1 and 5.3-2 there is 
a page in Appendix C, which is a reproduction of a computer printout page 
summarizing important results from one error budget run. 


5.3=3 Sensitivity Curves: System C 

This section contains several curves illustrating the sensitivity 
of System C performance to variations in major error source statistics. 
These curves were produced using the error budget data of Table 5.3-1. 
The method of generating the data for these curves is discussed in 
Section 5.1.3. Similar sensitivity curves may be constructed for any error 
source group for whi. h error budget data is available in Table 5.3-1. 
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TABLE 5,3-2 


SYSTEM C ALTERNATIVE CONTRIBUTIONS 





RMS Navigation Errors at Touchdown 


Error Source 

Value 

Position (ft) 

Velocity (ft/sec) 

Misalignment (mrad) 




Vertical 

Down range 

Crossrangc 

Vertical 

Downranga J Croaarajige 

Vertical 

Down range 

Cross range 

Group 1. 

Initial 3x9 Covariance 
Matrix {Position, Velocity 
and Alignment) 

See Section 3. 1.3 

0.3 

0.3 

0.9 

0.05 

0.01 0.05 

2.082 

0.030 

0.028 


Initial 6x6 Covariance 
Matrix (Position and Velocity 
only) 


0.2 

0.3 

6.2 

0,03 

0.01 0.5B 

1.933 

0.018 

0.024 

Group 15. 

Radar Altimeter Errors 
Bias 

Scale Factor 

Eq. (3.1-44) 
Eq. (3.1-45) 

1.0 

0.0 

0.1 

0.7 

1.4 

0.01 

0.12 

— 0.01 

0.02 0.01 

0.001 

0.001 

0.002 

0.005 
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Figure 5.3-4 shows the sensitivity of vertical position and 
velocity errors to accelerometer misalignments. The vertical velocity 
error does not meet the spec for any value of accelerometer misalignment, 
even though this group is the largest single contributor. The combined 
effect of all other groups is represented by the intercept of the curve with 
the ordinate; it is 0.23 ffc/sec, exceeding the spec value of 0.20. The 
straight line, representing the contribution of this group alone, shows that 
the rms misalignment would have to be less than 30 sec for the effect of 
this single error source to be less than the spec value. 

Figures 5.3-5 and 5.3-6 illustrate the effects of radar altimeter 
markov error and doppler rate bias error, respectively, on vertical posi- 
tion and velocity errors. The baseline radar altimeter rms error is 
taken to be a function of altitude, as given in Eq. (3.1-43). Thus, the 
ratio of the assumed function to the baseline function is used as the inde- 
pendent variable in Fig. 5.3-5. This is the largest contributor to vertical 
position error — a ratio of 1.3 causes the total curve to exceed the spec 
value. The effects of radar altimeter error and doppler rate bias error 
on vertical velocity error are approximately the same, as indicated by the 
lower halfs of the two figures. 

Figures 5.3-7, 5.3-8, and 5.3-9 show the sensitivity of cross- 
range position and velocity errors to three major contributors: antenna 
motion compensation error, accelerometer scale factor errors and gyro 
bias drifts. In all three cases the velocity error sensitivities are small 
relative to the spec value. The crossrange position sensitivities, how- 
ever, are large and the spec cannot be met by reducing any one of the 
three major contributors alone. 
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(a) Position 



(b) Velocity 


Figure 5.3-4 Sensitivity of Vertical Position and 

Velocity Errors at Touchdown to 
Accelerometer Misalignment : System C 


5-38 





THE ANALYTIC SCIENCES CORPORATION 


R-9697 



RELATIVE RMS RADAR ALTIMETER MARKOV ERROR 


(a) Position 



RELATIVE RMS RADAR ALTIMETER MARKOV ERROR 


(b) Velocity 

Figure 5.3-5 Sensitivity of Vertical Position and Velocity 
Error at Touchdown to Radar Altimeter 
Markov Error: System C 
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(a) Position 



(b) Velocity 

Figure 5=3-6 Sensitivity of Vertical Position and Velocity 
Errors at Touchdown to Onboard Doppler 
Rate Bias Error: System C 
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Figure 5 


R-9692 



(a) Position 



(b) Velocity 

3-7 Sensitivity of Crossrange Position and Velocity 
Errors at Touchdown to Antenna Motion 
Comoensation Error: System C 
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(a) Position 



(b) Velocity 

Figure 5=3-8 Sensitivity of Crossrange Position and Velocity 
Error at Touchdown to Accelerometer Scale 
Factor Errors: System C 
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TIME OF TOUCHDOWN 


BASELINE 

VALUE 




SPECIFICATION LIMIT 


RMS GYRO BIAS DRIFT {deg/hr) 

(a) Position 


TIME OF TOUCHDOWN 


SPECIFICATION LIMIT 


BASELINE 

VALUE 


Figure 5.3-9 


RMS GYRO BIAS DRIFT (deg/hr) 

(b) Velocity 

Sensitivity of Crossrange Position and Velocity 
Errors at Touchdown to Gyro Bias Drift: 
System C 
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5.4 SYSTEM D EVALUATION 

The results presented below correspond to the trajectory seg- 
ment, lasting 796 seconds, pictured in Fig. 2.1-1. The final five minute 
portion of the trajectory has the modified altitude time history described 
in Section 2.3 — ending in a steeper final approach during the last minute, 
compared to that used in the other three system evaluations. The mea- 
surement schedule is described in Section 2.2; it employs a 2 second 
interval between updates for 750 seconds, followed by a 0. 5 second update 
interval for the final 46 seconds. A file of gain values, generated using 
the filter covariance program outlined in Section 4.4, has been used re- 
peatedly to compute detailed error contributions corresponding to this 
trajectory and measurement schedule. 


5.4.1 Overall System Performance 

Figures 5.4-1 and 5.4-2 present overall performance curves 
for System D, showing position error components and velocity error com- 
ponents, respectively. The curves plotted represent rms errors due to 
the combined effects of all error source groups in the baseline error bud- 
get. They were generated by root-sum-squaring individual contributions 
at one minute intervals. (These are tabulated in Appendix C. ) At T DME , 
when the first DME and VOR measurements are made, and at T^^ , 
when the first baro altimeter measurement is made, the root-sum-square 
calculation was performed both before and after update. Thus, the large 
jumps in certain component errors, which occur at these two times, are 
accurately shown. Otherwise, the curves are faired -in over the one- 
minute intervals between the calculated points . 
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Figure 5.4-2 


System D Overall Performance: Velocity 
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At the time of the first DME and VOR measurements there are 
no large increases in any of the position or velocity component errors . 
This result contrasts with the other three candidate system evaluations, 
in which drag-updating was not used prior to the end of blackout. The 
vertical position and velocity errors do increase significantly in the next 
two minutes, prior to the first baro-altimeter measurement. Examina- 
tion of the individual contributions (tabulated in Appendix C ) shows that 
the major cause of these increases is the set of initial condition errors 
caused by non-standard atmosphere during the drag-update interval. The 
increase in crossrange velocity error over the first minute is caused by 
the VOR markov error. 

After the first baro-altimeter measurement all component 
errors appear to be of reasonable size and well behaved. (We are not 
speaking here of errors at touchdown, which cannot be judged from plots 
scaled as in Figs. 5.4-1 and 5.4-2. Final error performance of System D 
is treated in the following section. ) 

Section 6 contains further discussion of the implications of the 
performance curves presented in this section. 


5.4.2 Detailed Error Budget 

The System D baseline error budget is given in Table 5.4-1, 
showing rms estimation errors in position, velocity and platform align- 
ment components at touchdown. Each value is the rms contribution of 
the error source or sources indicated in the left hand column; it comes 
from a computer run in which the elements of those errors alone are set 
at truth-model values — and all others are zero. The overall System D 
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TABLE 5.4-1 

SYSTEM D BASELINE ERROR BUDGET 


Error Source 


J. ESTIMATED STATES AND NOISES 
Croup l. - Initial 6*6 Covariance 

klzlrlx l Position and Velocity) 

Initial Position and Velocity 
Errors Due 1o Non-Standard 
Almospherc during Drag 
Update 

Initial Prelaunch Aalmuth 
Misalignment (Optical) 

IMU Quantitation Nolae 
Measurement Nolseo 
DME 
VG.' 

Radar Altimeter 
First Order Markov Processes 
VOR 

Radar Altimeter 

[n. IMU-RELATED STATES 

Croup 2. Accelerometer Diases 
Croup 3. Accelerometer Scale 
Factor Errors 

Group 4. Accelerometer Misalignments 
Groups. Accelerometer NcnUnearitles 

Group fi. Gravity Anomalies and 
Deflections 

Croup 7. Gyro Dias Drill Ratss 
Group B. Gyro Mass Unbalances 
Group 9. Gyro Anlsoo last! city 

jm. EXTERNAL-AID RELATED STATES 
Group 10. DME Dias Error 
Group 11. DME Seals Factor Error 
Group 12. Raro-AlUmeter Errors 

Dias and 
Scale Factor 
Markov Process and 
Measurement Noise 

Static Defect 

Group 13. ILS Loeatizer Dias 
Group 14. 1LS Localizer Second-Order 
Markov Proccns 

Group 15. fNot Included in Baseline 
Error Budget) 

Group lfl. VOR/DME Survey Errors 
DME 
VOR 

Group 17. ILS Localizer Survey Errors 
Group 18. {Not In Sy Mem D Truth Model) 
Group 10, VOR Diaa Error 



RMS Navigation Errors at ToucWown 

Value 

Position 10) 1 

Velocity (ft /seel ! 

Misalignment (mrad) 


vmtoi 

Pmuu[t |Crfitirupj 

Ytrtlul 

DomLrJifcf* 

Cfolaittip 

v*rtit*J 


Crvtrup 

See Section 3.4.3 

0.1 

<3? 

o.s 

0.05 


0.04 

0.455 

0.010 

0.029 

See Section 3.4.3 

0.1 

0.7 

1.1 

0.04 

0.02 

<sI3> 

0.BB5 

0.049 

0.014 

00 sec 

_ 

— 

0.1 

— 

— 

0.01 

0.133 

0.005 

0.001 

One Pulse » 1 ern/sec 

0.3 

<33? 

0.5 

cos 

(pTB> 

0.08 

0.051 

0.037 

0.031 

24 ft 


<03 

0.2 

0.01 

0.07 

0.00 

0.000 

o.on 

0.009 


0.1 

<33? 

0.3 

0.03 

d£> 

<i>35> 

0.013 

0.021 

0.008 

3.28 ft 

<o? 

2.5 

— 

COD 

0.0B 

— 

0.007 

0.003 

0.003 

18.07 mrad 

0.1 

<G3> 

0.5 

0.05 

<od 

(6j3> 

0.027 

0.027 

0.021 

Eq.<3.1-43) 

<33? 

<S> 

0.1 

<£B> 

(STB) 

0.03 

0.012 

0.007 

0.005 

80 pg 

0.4 

2.3 

0,1 

COD 

0.05 

0.01 

0.088 

o.osi 

0.058 

34 ppm 

0.2 

1.2 

— 

0.10 

0.01 

— 

0.034 

0.017 

0.0D9 

40 tec 

<533 

<33? 

0.1 

<0> 

<5TT& 

0.02 

0.115 

0.172 

0.180 

3.5 hb/c 2 

<s? 

1.8 

0.9 

CB? 

0,08 

CS> 

0.020 

0.015 

0.003 

see Table 3.1-2 

0.3 

2.4 

0.2 

(£733 

<535> 

0.02 

0.042 

0.023 

0.019 

0.01 deg/hr 

. 

2.7 

0.3 

0.02 

<52? 

0.04 

0.233 

0.01B 

0.015 

0,015 den/tir/g 

— 

1.1 

0.1 

— 

0.05 

0.01 

0.072 

0.005 

0.012 

C.OOSdcE/W* 1 

— 

0.2 

— 

— 

0.01 

0.01 

0.028 

0.011 

0.017 

295 ft 

0.2 

<Sq3> 

0.2 

0.10 

0.04 

<ojD 

O.OQ*. 

0.022 

0.007 

100 ppm 

- 

1.8 

0.1 

0.01 

<OJ> 

— ‘ 

0.000 

0.001 

0.010 

100 ft \ 


0.1 


0.01 


__ 

0.013 

0.003 



0.03 ft 1 










1 20 ft I 







__ 

— 

— 

5 ft J 










1.52.10 -4 rt/(n/..r! z 

. — 

— 

— 

— 

— 

— 

0.004 

0.001 

■ * 

0.5 mrad 

— 

<o? 

<o> 

— 

0.03 

ds> 

0.101 

0.032 

0.002 

Eq. (3.1-44) 

— 

4.3 

<53? 

0.01 

0.05 

dS§> 

0.107 

0.051 

0.000 


1 ft 
1 ft 
I ft 

5 mrad 


Total Projected Performance Imt* -sum-square) 


0.8 

— 0.8 

— - 0.7 

— 0.2 


0.1 

1.1 


Fitler-lndtratcri Performance 


Speclfirallnn (Ret. 11 


— 0.01 

— 0.01 

— 0.01 

— 0.01 


0.03 

0.03 


0.001 

0.003 

0.003 

0.002 


0.005 

0.005 


0.001 

0.001 

0.001 


— 0.001 


1.153 0.211 


0.122 0.211 


Not Specified 


A " — '* entry rrfrrs to a negligibly small contribution. 
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performance, or root-sum -square values, are given at the bottom of the 
table and compared, first with the filter’s own estimate of performance 
based on its simplified 23-state model of system error dynamics, and 
second, with the Shuttle landing navigation specification given in Ref. 1. 
Particular numbers have been circled in the posits jn and velocity error 
colu mn s of Table 5.4-1 to focus attention on major contributors. The 
rule followed in deciding which numbers to circle was to include, in any 
one column, the largest value, second largest value, etc., until 95% of 
the RSS total value was reached. The following general statements may 
be drawn from an examination of the System D baseline error budget: 

® The projected downrange performance , in both 
position and velocity, is considerably more 
accurate than called for in the preliminary 
specifications; and, in fact, is more accurate 
than filter -indicated performance. Significant 
contributors are DME errors (Groups 1, 10 and 
11), VOR errors (Group 1), initial condition 
errors (Group 1), IMU quantization noise 
(Group 1), radar altimeter errors (Group 1), 
accelerometer misalignments (Group 4), gravity 
anomalies (Group 6) and localizer bias 
(Group 13). 

» The projected vertical performance is within spec 
in position, but out-of-spec in velocity. The 
filter -indicated performance is quite accurate, 
being slightly pessimistic in vertical position and 
slightly optimistic in vertical velocity. The 
largest contributors to vertical velocity errors 
are accelerometer misalignments, nonlinearities 
and biases (Groups 4, 5 and 2) and radar altimeter 
markov error and measureme nt noise (Group l). 

Other significant contributors to vertical errors 
are IMU quantization (Group 1) and gravity 
anomalies (Group 6). 
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® The projected crossrange performance is within 
spec in velocity, but out-of-spec in position. 

The filter -indicated performance, in both position 
and velocity, is somewhat optimistic — predicting 
approximately two -thirds of the projected actual 
rms errors. The largest contributors to cross - 
range errors are localizer bias and markov 
errors (Groups 13 and 14). Other significant 
contributors to crossrange velocity error are 
initial condition errors (Group l),VOR errors 
(Group 1), accelerometer nonlinearities (Group 5) 
and D ME bias (Group 10). 

e The projected errors in estimating platform 
misalignments are approximately 0.2 milli- 
radians for misalignments about level axes, 
and 1 . 2 milliradians for azimuth misalign- 
ment. The filter -indicated performance is 
accurate in tilt-about-crossrange and some- 
what optimistic in both azimuth and tilt-about- 
downrange — predicting approximately one- 
half the projected actual performance for these 
two components. 

Considering all components of position and velocity errors at 
touchdown, the list of major contributors includes gyro errors, accel- 
erometer errors, and errors associated with all of the external devices 
except the baro altimeter. Evidently, there are strong correlations 
between various components of error at touchdown. This accounts, for 
example, for VOR errors contributing significantly to downrange errors 
and DME bias contributing significantly to crossrange velocity error. 

The overall performance projected for System D is similar to that of 
System A, except for the larger rms vertical velocity error. The major 
reason for the latter is the steeper final approach trajectory used in the 
System D evaluation. This result and prospects for improving System D 
performance are discussed in Section 6. 
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Table 5.4-2 lists the System D Alternative Contributi ons show- 
ing how different assumptions or additional error sources wou*d produce 
other contributions to position, velocity and alignment errors at touch- 
down. Detailed tabulations of the error contribution time histories are 
given in Appendix C. For every row in Tables 5.4-1 and 5.4-2 there is 
a page in Appendix C, which is a reproduction of a computer printout page 
summarizing important results from one error budget run. 


5.4.3 Sensitivity Curves: System D 

This section contains several curves illustrating the sensitivity 
of System D performance to variations in error source statistics. These 
curves were produced using the error budget data of Table 5.4-1. The 
method of generating the data for these curves is discussed in Section 5. 1. 3. 
Similar sensitivity curves may be constructed for any error source group 
for which error budget data exists in Table 5.4-1. 

Figures 5.4-3 and 5.4-4 show the sensitivity of downrange posi- 
tion error to variations in the two largest contributors, DME bias error 
and VOR markov error, respectively. In both cases the overall perform- 
ance stays well within spec for substantial increases in the assumed error 
source rms values. 

Figures 5.4-5 and 5.4-6 show the sensitivity of vertical posi- 
tion a nd velocity errors to variations in accelerometer misalignment and 
radar altimeter markov error, respectively. The baseline performance 
is within spec in position, but out-of-spec in velocity. Accelerometer 
misalignments are the major contributor to vertical velocity error at 
touchdown, However, even for zero accelerometer misalignment the 
velocity error remains out-of-spec due to the combined effect of the other 


5-50 



THE ANALYTIC SCIENCES CORPORATION 




5-51 




THE ANALYTIC SCIENCES CORPORATION 


R- 9703 



Figure 5.4-3 Sensitivity ofDownrange Position 

Error at Touchdown to DME Bias 
Error: System D 


R-9698 



Figure 5.4-4 Sensitivity ofDownrange Position 
Error at Touchdown to VOR 
Markov Error : System D 
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R-9699 



(a) Position 



(b) Velocity 


Sensitivity of Vertical Position and Velocity 
at Touchdown to Accelerometer Misalign- 
ment: System D 


Figure 5 , 4-5 
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RELATIVE RMS RADAR ALTIMETER MARKOV ERROR 


(b) V elocity 

Figure 5 4-6 Sensitivity of Vertical Position and Velocity 
Errors at Touchdown to Radar Altimeter 
Markov Error : System D 
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error sources. An rms misalignment accuracy of 20 sec or better is 
required to make this single contribution less than 0.2 ft/sec, the spec 
value a 


Figures 5 = 4-7 and 5.4-8 show the sensitivity of crossrange 
position and velocity errors to localizer markov and bias errors, 
respectively o The baseline performance is well within spec in velocity, 
but out -of -spec in position. In order for the combined effect of these two 
major contributors to be within spec both would have to be reduced to 
approximately one-half their baseline values. 

Summary — Projected overall position and velocity component 
errors have been plotted versus time for the four candidate systems. 
Detailed error budget tables and sensitivity curves related to contribu- 
tions at touchdown have also been presented . All four systems are shown 
to be out-of-spec in crossrange position error and vertical velocity error 
at touchdown. System B is also slightly out-of-spec in downrange posi- 
tion and velocity errors at touchdown. The following section summarizes, 
compares and interprets these results. 
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(a) Position 



(b) Velocity 


Figure 5,4-7 Sensitivity of Crossrange Position and 
Velocity Errors at Touchdown to 
Localizer Markov Error : System D 
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RMS LOCALIZER BIAS ERROR (mrod) 


(a) Position 



RMS LOCALIZER BIAS ERROR (mrad) 


(b) Velocity 

Figure 5.4-8 Sensitivity of Crossrange Position and Velocity 
Errors at Touchdown to Localizer Bias 
Error: System D 
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6 C SUMMARY COMPARISON AMD DISCUSSION 


In this section the results of the four candidate system evalua- 
tions are summarized, compared, interpreted and discussed. Overall 
projected performance at touchdown is compared in a single table, which 
also lists the dominant error sources for each system. The general 
implications of these touchdown -point comparisons as well as compari- 
sons of performance at earlier times of interest are reviewed. Related 
groups of error sources are then taken up one at a time, and their dif- 
ferent contributions to system errors are contrasted and explained. 
Finally, the choices of filter states are reviewed and possibilities for 
improving performance through software modifications are discussed. 


6 1 OVERALL PERFORMANCE COMPARISON 

Table 6,1-1 presents a summary of the projected rms touch- 
down errors for the four systems. The table is divided vertically into 
four major sections corresponding to Systems A, B, C and D, respec- 
tively, The top two rows in each division are the overall system per- 
formance projections taken from the bottom of the Baseline Error Budget 
tables given earlier. Numbers greater than the preliminary specifica- 
tion values given in Table 1.2-1 represent out-of-spec performance and 
are boxed-in. The rows labeled "filter performance ratio” relate the 
filter -indicated performance, based on the filter's covariance calcula- 
tion, to the projected system performance, based on the higher -dimensional 
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TABLE 6.1-1 

COMPARISON OF PROJECTED SYSTEM ACCURACY AT TOUCHDOWN 


j Vertical 

Downrango 

Crossrange 

SYSTKM A 




RMS Errors, eg 

Position (ft) 

1.3 

28.7 

OlpJ 

Velocity (ft/sec) 

[ 0. 22 | 

1.07 

0.49 

Filter Performance 
Ratio, a j,’/ {7 s 

Position 

1,00 

1.04 

0.C9 

Velocity 

0.73 

0.84 

0.09 

Dojninant Error Sources 

Glide Slope Errors 

Accelerometer 

Misalignments 

Radar Altimeter 
Errors 

Gyro Bias Drifts 

Localizer Errors 

SYSTEM B 




RMS Errors, ug 

Position (ft) 

2.2 

|85. 6 | 

1 9.U 

Velocity (ft/sec) 

|0.34i 

1 3.601 

0.89 

Filter Performance 
Ratio, crp/og 

Position | 0.64 

0.25 

0.45 

Velocity 

0.95 

0.13 

0.29 

Dominant Error Sources 

Radar Altimeter 
Errors 

DME Errors 

DME Errors 

Localizer Errors 

SYSTEM C 




RMS Errors, eg 

Position (ft) 

2.6 

3.0 

110-0 J 

Velocity (ft/sec) 

|0. 39 1 

0.21 

1.01 

Filter Performance 
Ratio, Op/o s 

Position 

.0.81 

0.20 

0.35 

Velocity 

0.54 

0.29 

0,33 

Dominant Error Sources 

Accelerometer 

Misalignments 

Radar Altimeter 
Errors 

Doppler Rate 
Bias Errors 

Accelerometer 

Biases 

Accelerometer Scale 
Factor Errors 

Antenna Motion 
Compensation Errors 

SYSTEM D 




RMS Errors, Og 

Position (ft) 

2. 3 

23.5 

|8. 4 | 

Velocity (ft/sec] 

10-021 

0.57 

0.44 

Filter Performance 
Ratio, <Jp/og 

Position 

1.09 

2.30 

0.68 

Velocity 

0.81 

2.34 

0.68 

Dominant Error Sources 

Accelerometer 

Misalignments 

Radar Altimeter 
Errors 

Accelerometer 
Nonlinear ities 

DME Errors 
VOR Errors 

Localizer Errors 


■1 „ J denotes out-of-spec performance. 

* ■ 
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truth model covariance analysis . Values less than unity indicate opti- 
mism on the part of the filter. The other entries list the dominant 
(single largest or two or three largest) error sources — those making 
major contributions to vertical, downrange or crossrange errors at 
touchdown. 

The most striking feature of the results summarized in 
Table 6.1-1 is that all four systems are out-of-spec in crossrange posi- 
tion and vertical velocity. While System B is also shown to be out-of- 
spec in downrange position and velocity, this is not of great concern, as 
discussed below. First, we discuss the two components which have con- 
sistently emerged as troublesome. 

The projected crossrange position errors for the three con- 
ventional navaid systems (A, B and D) are approximately equal and are 
all dominated by the same error sources -- the ILS localizer errors. 

This is true despite the presence of independent crossrange data, coming 
from the second DME station in Systems A and B, and from the VOR data 
in System D. These independent data sources are not nearly as accurate 
as the localizer data and, therefore, do not really help in reducing cross- 
range errors at touchdown. A related observation is that Systems A and 
D perform only marginally better than System B, even though the former 
include a state corresponding to correlated localizer errors. The state 
estimate is not accurate because of the lack of accurate, independent data. 
In summary, the crossrange performance of these conventional navaid sys- 
tems is only as good as the localizer data; * sophisticated software tech- 
niques are not likely to change this result. 

The projected crossrange position error for System C, the one- 
way doppler system, is 10 to 20 percent higher than that of the other three 

The use of some independent source of accurate data., such as precision 
ranging systems, would of course be effective in reducing the error. 
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systems,, However, when the major contributors (accelerometer scale 
factor errors and antenna motion compensation errors) are considered, 
we find that there are genuine possibilities for improvement through 
careful software design; these are discussed later in this section and in 
6. 2 and 6. 3. In other words, the candidate system with the largest pro- 
jected crossrange error is the one most likely to be able to achieve the 
spec, given sufficient filter design improvement effort- We note, also, 
that a one-way doppler system based on two or more ground beacons, 
rather than the single station of System C, is likely to be capable of much 
improved crossrange performance- On the other hand, a potential short- 
coming of the one-way doppler scheme, especially a single -station ver- 
sion, is that it may not work well with a , 'straight-in , ’ approach- That is, 
the spiral approach path used in this evaluation may produce important 
correlations which account for the system’s relatively good crossrange 
navigation performance. 

Although the projected vertical velocity errors exceed the spec 
for all four candidate systems, System A comes close — exceeding the 
prelimr ary specification value by ten percent- The System A filter 
generally strikes a good balance between glide slope data, IMU data and 
radar altimeter data, as shown by the fact that elements of all three are 
major (therefore somewhat equal) contributors to vertical errors. The sen- 
sitivity curve of Fig. 5- l-4(b) shows that a tighter IMU hardware specifica- 
tion, reducing the allowable rms accelerometer misalignment, could reduce 
the projected overall vertical velocity error at touchdown essentially to the 
spec level — with no other changes- It may also be possible to gain im- 
proved vertical performance by filter design modifications, as suggested 
in Section 6-3„ Since System B consists of the same hardware elements 
as System A, it is of course possible to reduce its errors to the level 
shown for System A through the addition of filter states and proper choice 
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of filter parameters, A careful comparison of the two error budgets 
indicates that only a few additional states are needed, as discussed in 
Section 6.3. Conversely, quite a large number could be safely omitted 
from the System A filter. 

Before discussing the vertical performance of Systems C and 
D, we make the following general observation. Radar altimeter errors 
aio a dominant source of vertical errors for r.ll four candidate systems. 

Since the radar altimeter is a direct source of vertical data, it is clear 
that any given level of vertical navigation accuracy is possible if this 
data alone can be made sufficiently accurate arid if the filter is designed 
to place sufficient weight on it. 

Vertical velocity err ore are also a problem for System C, the 
single-station one-way doppler system. Accelerometer misalignments 
and other IMU error sources contribute to vertical velocity errors to a 
somewhat greater degree than for System A. Comments already made 
about improving vertical velocity performance also apply to System C. 

That is, a tighter IMU accelerometer misalignment spec, better radar 
altimeter performance and filter design m .>tiifications (discussed further 
in Section 6.3) can all be significant. 

Comparison of the projected touchdown performance of 
System A and System D, shown in Table 5.5-1, shows that the only im- 
portant difference is in their vertical velocity accuracies. The System D 
rms vertical velocity error is nearly three times that of System A, There 
are two reasons for this difference: the lack of glide slope measurements 
in SystemD and the steeper final approach path used in the SystemD evalua- 
tion. The steeper path affects vertical velocity accuracy in several ways. 

The radar altimeter measurements begin later in the flight, there are less of 
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them and they are on the average less accurate (since the assumed rms 
error is a function of altitude, as given by Eq. (3.1-43)). Also, since 
the vertical velocity is changing more rapidly, the filter is less able to 
estimate it accurately from the sequence of vertical position measure- 
ments. An additional filter covariance calculation was made to gain in- 
sight into this problem (see Section 4.4). The results indicate that the 
steeper approach path is probably the major cause of increased vertical 
velocity error for System D . This is perhaps the most negative finding 
of this study -- calling into question the ability of the new Orbiter con- 
figuration to be landed safely and automatically with conventional navaids. 
To do so requires that the radar altimeter accuracy be considerably 
better than assumed in our baseline studies. 

The projected downrange performance for the three systems 
using high -dimensional filters (A, C and D) are well wiSin specs. 

System B, which includes a 6 -state filter, is slightly out-of-spec in both 
downrange position and velocity accuracy. Since DME errors are the 
do minant contributors to downrange errors for System B, it seems clear 
that the addition of two DME bias states, as included in System A, will 
provide the necessary improvement. (This point is amplified in 
Section 6.2. ) MIT/DL personnel have reported on simulations of a 
6 -state filter which satisfies the downrange specification. This is ac- 
complished by proper selection of assumed measurement error variances 
to achieve more nearly optimal mixing of glide slope, DME and radar 
altimeter data. They have also simulated further improvements through 
the use of additional states, such as DME bias. In summary, meeting 
the preliminary specifications for downrange performance does not appear 
difficult for the hardware combinations considered in this study. 
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The remainder of this section is devoted to a summary of 
points related to earlier portions o f the trajectory. The drag -updating 
procedure, assumed in the System D evaluation, takes places during the 
blackout portion of the trajectory and has a substantial effect on the 
"initial conditions" of the post-blackout portion. An obvious benefit to 
System D, as compared to Systems A and C, is evident in the error time 
histories of Fig. 5.4-1, as compared to those of Figs. 5.1-1 and 5.3-1. 

The initial errors are smaller and do not experience large increases — 
due to incorrect correlation assumptions — at the time of the first ex- 
ternal measurements. A more subtle benefit involves the accelerometer 
scale factor errors (discussed in Section 6. 2) which lead to incorrect azimuth 
alignment estimates in Systems A and C. This effect is absent in System D 
due to the changed initial correlation matrix. The inclusion of drag- 
updating, particularly in the case of System C(one-v/ay doppler) should 
lead to significant improvement in crossrange error at touchdown. 

A significant advantage of the VOR/DME measurements used 
in System D, compared to the two DME stations used in Systems A and B, 
involves the avoidance of filter initialization problems. These are 
numerical problems which can arise when large navigation corrections 
are made. Situations of poor measurement geometry, such as when all 
of the accurate measurements are associated with a given "direction" in 
state space, are particularly vulnerable. Nonlinear effects can also play 
a role in these difficulties — for example, when an accurate range mea- 
surement is made in the presence of large uncertainties in perpendicular 
directions. Filter initialization problems have not occurred in generating 
the four baseline error budgets reported herein, using linear covariance 
methods. Others, however, have reported such difficulties in nonlinear 
simulation work. It is clear that the "angle -range" data from a single 
VOR/DME station always provide "good geometry" while the "range-range" 
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data from two DME stations provide relatively poor geometry in the early 
stages of the trajectory and for the station locations employed in this 
study o The position error time histories for System D plotted in 
Fig. 5.4-1 reflect this geometrical advantage in the ,T midcourse" portion 
between t = 5 minutes and t = 10 minutes. The errors are consistently 
low in this region as compared to those of the other three systems pre- 
sented in Figs. 5.1-1, 5.2-1 and 5.3-1. 

As discussed in Section 5.1.1, System A has a problem involv- 
ing large vertical errors just prior tc the time of the first radar alti- 
meter measurements. The problem is associated with mis -estimates of 
platform tilts about level axes. System B apparently avoids the need for 
accurate platform date, of this type by putting more faith in the baro 
altimeter in this region. System D minimizes the problem by actually 
having better tilt estimates in this region. This is apparently another 
advantage of the good geometry of the VOR/DME measurements. 

The following sections take up individual groups of error 
sources, one at a time, and discusses and compares their effects on the 
four candidate systems. 


6.2 COMPARATIVE INDIVIDUAL CONTRIBUTIONS 
6,2.1 IMU Errors 

Gyro bias drifts a^e major contributors to downrange errors 
for Systems A and B and significant, but relatively minor, contributors for 
System D. The System B contributions are somewhat larger than those 
of System A (21.4 ft and 1.11 ft/sec compared to 17.8 ft and 0.85 ft/sec), 
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even though the System A filter weights IMU data more heavily during the 
final approach o Examination of the detailed results and comparison with 
the contributions of initial errors (given in the Alternative Error Con- 
tributions, Tables 5.1-2 and 5.2-2) shows that these effects are not 
primarily associated with the drifting that occurs during approach and 
landing, but are related to the ’’initial” platform tilt (at 130, 000 ft) caused 
by drifting during the major portion of the one-rev mission. The better 
performance for System A is due to its estimation of platform tilt. The 
actual rms tilt-about-crossrange is 127 sec (System B Error Budget; 

0.637 mrad), while the rms estimation error is 4 sec (System A Error 
Budget; 0.018 mrad). On the other hand, System A's estimates of the 
drift rates themselves are poor and do not help. (They don't hurt either 
because the filter "knows” that these estimates are poor, as shown in 
Table 4.1-2. ) Gyro bias drifts are also significant contributors to cross - 
range errors for System C (3.8 ft) and B(2.0 ft). The filters of Systems A, 

C and D do a reasonably good job of estimating tilt-about-downrange, which 
helps to hold these contributions down. For these three systems the azi- 
muth misalignment is a greater factor in producing crossrange error at 
touchdown, and System C (one -way doppler) is affected the most. Other 
significant, but not serious, effects of gyro oias drifts involve vertical 
velocity errors in Systems A and C. Gyro mass unbalances ( g-sensitive 
terms) also produce significant, but not serious, vertical errors for 

System A. 

Accelerometer misalignments are dominant contributors to 
vertical velocity errors for Systems A, C and D (0.10, 0.26 and 0.36 ft/sec, 
respectively). The key to understanding these effects is the realization 
that there are six independent misalignment components, two for each in- 
strument, resulting in a nonorthogonal triad of sensitive (input) axes. 
Therefore, although the orientation of the (initially) level instrument 
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sensitive axes may be well known, the vertical accelerometer sensitive 
axis direction will not be. This results in an incorrect resolution of the 
horizontal deceleration existing prior to touchdown, producing vertical 
velocity and position error growth. Between measurements these effects 
are about the same for System A and System B„ But, the overall effect 
on System A is greater because it places, more faith in the compensated 
IMU outputs. (The same is true for Systems C and D. ) This faith is 
based on the 18 states of the filter which relate directly to the IMU. 
Sophisticated as this 18-state model is, it does not take nonorthogonality 
into account. Possible approaches to improving the filter in this regard 
are discussed in Section 6,3. Other significant, but not serious, effects 
of accelerometer misalignments involve downrange errors for Systems A 
and D, and crossrange errors for P ;tem C . 

Accelerometer scale factor errors contribute to crossrange 
errors in a dominant way for System C and significantly for System A 
(5.0 and 1.8 ft respectively). Examination of the detailed results shows 
that this effect is due to a mis-estimation of azimuth misalignment . The 
azimuth estimation error due to this one error source group is 1.688 
milliradians for System C, compared with a total actual azimuth mis- 
alignment of 0.620 milliradians due to all sources, as shown in the 
System B error budget. Apparently Systems A and C would perform better 
if they ignored their own azimuth updates. System D, on the other hand, 
avoids this problem because of the altered initial correlations due to drag- 
updating, ""he scale factor errors also produce significant, but not serious, 
contributions to vertical and downrange e rrors for System A. 

Accelerometer nonlinearities produce negligible effects for the 
first three systems, but surprisingly large effects for System D. The 
vertical velocity error contribution (0.24 ft/sec) is among the dominant 
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ones for this system . This result is apparently related to the steep 
approach path and more abrupt flare maneuver associated with the new 

configuration. 

Accelerometer biases are major contributors to vertical 
velocity error for System D (0.21 ft/sec) and significant, but not serious, 
contributors for Systems A and C, in both vertical and downrange^ direc- 
tions. The size of the contributions compared with those of other signi- 
ficant error sources indicate that these systems generally strike a good 
balance between inertial errors and external measurement errors. At- 
tempts to estimate the accelerometer biases in Systems A, C and D are 
not effective. (In a one-g flight condition the level biases are not dis- 
tinguishable from tilts and level axes. It is only maneuvering accelera- 
tion which provides any hope of separately observing the accelerometer 
errors. ) 

IMU Quantization produces significant effects for Systems A, 

C and D; most notably the crossrange position error (3.4 ft) for System C. 
These are the systems with high-dimensional filters. They place more 
faith in IMU data than does System B, and are therefore more affected by 
the quantization. This faith is justified because of their better knowledge 
of platform tilts about level axes. 

The interpretation of the accelerometer scale factor effects 
given above and further examination of the detailed error budget tables, 
leads to a set of imp ’lied requirements for initial azimuth alignment accu- 
racy. In order to hold the IMU contribution to rms crossrange error at 
touchdown to approximately 2 feet or less, the initial rms azimuth mis- 
alignment should be approximatelv 2 milliradians for conventional navaid 
systems and approximately 1 miiliradian for a single station one-way 
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doppler system. These requirements relate to conditions at the start of 
the reentry phase and apply regardless of the method used to achieve 
alignment prior to reentry. 


6.2.2 Conventional Navaids 

DME bias errors are the dominant source of downrange errors 
for System B (59.2 ft and 3.13 ft/sec) and significant sources for Systems 
A and D. This result clearly stems from the fact that DME biases are 
not estimated in the System B filter while the other filters estimate them 
quite effectively. The detailed results in Appendix C for System A, 

Group 10 show that the rms estimation errors for the DME filter states 
are 292 and 294 feet for Stations 1 and 2, respectively, compared to actual 
rms bias errors of 285 feet. This implies that the filter effectively cali- 
brates the DME biases with errors of seven and nine feet. It is ob- 
viously useful to calibrate Station 1 (the station directly under the 
approach path) with this degree of accuracy, since this station provides 
the only good measure of downrange position in the final 20 or 30 seconds 
before touchdown. It is interesting to note that the DME noise • ontr ibu - 
tions for Systems A and D are greater than for System B. This reflects 
the fact that, having calibrated the correlated portion of DME errors, the 
higher -state -size filters can (and should) put more faith in these mea- 
surements -- striking a good balance between the effects of the residual 
DME error an^ other error sources. The DME biases also contribute in 
a major way to vertical errors for System B, and in a sign. 1 i Leant, but not 
serious way to c rossrange velocity error for System D. The latter effect 
is an illustration of the correlation previously noted between downrange 
and crossrange errors in the System D (VOR/DME system) evaluation. 


6-12 



L 


i 

i 


I 


THE ANALYTIC SCIENCES CORPORATION 


DME scale factor errors are major contributors to vertical 
errors for System B (0.3 ft and 0.19 ft/sec) but relatively minor con- 
tributors for Systems A and D. This difference is a result of the lack 
of DME error states in the System B filter and the dynamic interaction 
between downrange and vertical errors in the final approach pnase. 

The VOR system is employed in System D only. The VOR 
maikov error i s a significant contributor to crossrange error for that 
system and is listed as a dominant contributor to downrange errors. 

The latter contributions (11.0 ft and 0.32 ft/sec) are not serious; they 
are listed as dominant sources unly because the overall system down- 
range performance is very good — there are no large contributions com- 
pared to the spec values. V OR measurement noise also produces signi- 
ficant, but not serious, downrange and crossrange errors for System D. 

The ILS glide slope bias is a dominant contributor to vertical 
errors for System A and a sip 7 ’ *icant but much ? nailer contributor for 
System B. The glide slope markov error also produces larger contri- 
butions for System A than for System B. This is actually a fault of the 
System B filter which should give more weight to glide slope measure- 
ments, striking a better balance between glide slope and radar altimeter 
errors. The glide slope errors also produce significant, but not serious, 
contributions to downrange errors for Systems A and B. 

Localizer error contributions, whether modeled as second- 
order markov processes or bias errors, are remarkably similar for the 
three conventional navaid systems; even though the System A and D fil- 
ters include a localizer error state and the System B filter does not. 

The c rossrange position contributions due to the second-order markov 
error are 5.9, 6.9 and 6.2 ft for Systems A, B and D, respectively. This 
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result reflects the fact that localizer error is essentially unobservable 
in the postulated situation — there is no independent source of accurate 
crossrange measurement, and the filter is forced to treat localizer mea- 
surements as unbiased. These contributions can only be reduced by im- 
proving the ILS system itself — not by filter design changes. 


6.2.3 One-Way Doppler Errors 


The doppler beacon is included in System C only. Doppler rate 
bias error is listed as a dominant contributor to downrange errors (1.8 ft 
and 0.15 ft/sec) and a significant contributor to vertical velocity error 
(0.09 ft/sec). The former are not important since they are very small 
compared to the downrange specification values. The postulated rapidly 
varying rate bias error (30 second correlation time — Group 10) is not 
a significant contributor. It is of interest to note that if the postulated 
rms value were scaled up to match that of the slowly varying component 
(400 second correlation time — Group 1), the system error contributions 
would be approximately the same for both sources. Thus, it is not at all 
critical for the correlation time selected for this error state by the filter 
designer to be an accurate representation of the real-world error charac- 
teristic . 


The antenna motion compensation error is a dominant contri- 
butor to crossrange errors (5.6 ft and 0.65 ft/sec) and a significant con- 
tributor to vertical and downrange errors. The effect on crossrange 
errcrs at touchdown is similar to that of the accelerometer scale factor 
errors discussed previously. During banking maneuvers, when the mea- 
surement errors are significant, the filter mis -interprets the observations 
and develops a pl/itform azimuth alignment estimation error. (Note the 
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1.171 milliradian entry in the System C error budget.) This persists 
into the final approach phase and leads to significant crossrange errors 
through a mis -resolution of the deceleration prior to landing. 


6. 2.4 Altimeter Errors 

The baro-altimeter errors produce negligible touchdown errors 
for System A, C and D and significant, but not serious, downrange and 
vertical errors for System B. It is worth mentioning again (see 
Section 6.1) that System B out-performs System A, in the region just 
prior to the start of the radar altimeter measurements, because it makes 
more use of the baro-altimeter measurements. 

The radar altimeter errors are dominant contributors to 
vertical errors at touchdown for all four systems. This is an anticipated 
result since the radar altimeter is a direct source of accurate vertical 
data in the final approach phase. Significant contributions to downrange 
errors are also produced for Systems A, B and D, the conventional 
navaid systems. As mentioned previously the larger vertical velocity 
contributions for System D appear to stem from the steeper final approach 
path associated with the new aerodynamic configuration. 


6.2.5 Survey Errors 

The assumed rms value of all survey errors used in generating 
the baseline error budgets was 1.0 ft (for each of three components). In 
most cases this produced an easily predicted result such as the 1.0 to 
1.1 ft crossrange errors due to localizer survey error for Systems A, B 
and D, or the 1.0 ft downrange error due to doppler beacon survey error 
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in System C. The only significant touchdown error contributions due to 
survey errors are the downrange position errors (7 . 1 and 14 . 4 ft) due to 
glide slope survey errors for Systems A and B. Examination of the de- 
tailed computations associated with these results shows that most of these 
contributions stem from the vertical component of survey error. The 
implied requirement is that the glide slope antenna should be surveyed - 
in carefully enough to insure that the accuracy assumed herein is met, 
particularly with respect to the vertical displacement from the runway 

level. 


6.3 FILTER DESIGN 

This section contains comments and suggestions pertaining 
to the design of the filters which form a part of each of the candidate 
systems. Many of these comments have been mentioned or implied 
in preceeding sections j some are mentioned below for the first time. First, 
it should be noted that this investigation has not been a filter design study 
per se. Such a study would involve successive calculations with specific 
filter states added or deleted, and would be aimed at reaching a snecific 
recommended sub-optimal filter design for a given hardware combination. 
However the detailed results of this study, along with their interpre- 
tation, do provide many useful insights and inputs to such a design effort. 

This is especially true with respect to the comparative results for Systems A 
and B, which have identical conventional navaids (in the hardware sense) 
but very different filter designs. In fact, the System B filter was de- 
liberately chosen to have only six states (the minimum conceivable num- 
ber for this application) in order to make this comparison as useful as 
possible — to indicate what additional states are most needed to achieve 
"performance comparable to that of the 24-state filter of System A. 
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Based on the comparative individual error contributions dis- 
cussed above some or all of the following filter states could usefully _ bg_ 
added to System B: 

® 2 DME Error States (definitely one for Station 1, 

probably one for Station 2 also) 

® 1 Radar Altimeter Error State 

« 2 Platform Level Tilt States 

The Station 1 DME error state is clearly useful in bringing downrange 
errors within spec. It is not so obvious, from the results presented, 
that the Station 2 state is needed, but it may be important to have both 
states included in order to achieve the kind of calibration accuracy in- 
dicated in the System A results for either one. There may also be good 
reason to include both when consideration is given to failures and re- 
dundancy management questions. The addition of a radar altimeter error 
state probably will reduce the radar altimeter contribution to vertical 
velocity error for System B. On the other hand, if the suggestion re- 
garding a vertical acceleration error state given later in this section is 
followed, this extra state may not be needed. The addition of the level 
tilt states will clearly reduce significant gyro drift contributions to down- 
range and crossrange errors. The resulting improvement in overall 
performance may or may not justify these additions. Perhaps the best 
argument for including them is that the gyro drift rate specifications, for 
this kind of one-rev pure -inertial-navigation (prior to approach) scenario, 
could be relaxed. 

Based on the comparative individual contributions discussed 
above there are two places where major System A contributions arc signi- 
ficahtlv larger than the corresponding System B contributions. These are: 
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» The contribution of accelerometer misalignments 
to vertical velocity error (0.10 ft/sec — note the 
even larger contributions for System C and D). 

® The contribution of accelerometer scale factor 
errors to crossrange position error (1.8 ft — 
note the even larger contribution for System C). 

A brute-force approach to the reduction of the accelerometer misalign- 
ment effect would be to add three additional filter states representing the 
nonorthogonality of the triad of accelerometer sensitive axes. A simpler 
approach would involve the use of an artificially high value of velocity 
state noise in the vertical direction* This would effectively result in less 
weight being given to the IMU vertical accelerometer outputs, and more 
weight given to the glide slope and radar altimeter measurements. A 
third approach, which might be the most fruitful of all, is given later in 
this section. The contribution of accelerometer scale factor errors to 
crossrange position error is almost entirely due to the mis -estimate of 
azimuth misalignment. Systems A and C would perform better if they did 
not att empt to estimate the azimuth component of misalignment. Alter - 
natively, there may be value in maintaining the azimuth estimation as part 
of a failure detection scheme, for example, but in omitting the IMU com- 
pensation associated with it as long as it remains within some reasonable 
level. It should be recalled, also, that when drag -updating was included 
in the System D evaluation, no appreciable azimuth misalignment or cross- 
range error was caused by the accelerometer scale factor errors. Thus, 
the entire matter is largely a function of initial correlations as determined 
by events prior to the end of blackout. 

Based on a detailed examination of the error budget calculations 
for the first two systems it appears that System A landing performance 
would not suffer if the following states were omitted from its filter: 
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o 3 gyro drifts 

• 3 accelerometer scale factor errors 

* 3 accelerometer biases 

» 1 ILS localizer error 

® 1 ILS glide slope error 

These states do not appear to help or hurt the System A performance in 
any serious way — see the discussions in Sections 6.2.1 and 6.2.2. (It 
is recognized that the gyro and accelerometer states may be useful in other 
mission phases. It is also possible that they can be useful during approach 
and landing as part of a failure detection scheme. ) The above statements 
apply equally well for System D. That is, the nine IMU sensor error 
states and the ILS localizer error state could be omitted without appre- 
ciable effect on touchdown performance. 

It is possible that the nine IMU states can also be omitted from 
the System C (one-way doppler) filter without loss of accuracy, but the 
calculations presented in this report do not provide sufficient evidence to 
justify such a claim. It is clear that for a single -station one-way doppler 
system there is no advantage in including two separate doppler rate bias 
error states, since they are not separately observable. (For a multiple- 
station doppler system, however, separate states for each station plus an 
added state for the onboard clock error do make sense, since ali are 
theoretically observable.) The large contribution of antenna motion com- 
pensation error to crossrange position error for System C can be reduced 
by one or both of two approaches. First, the error can be reduced at the 
source simply by designing a highly accurate compensation scheme. 

Second, the contributions could probably be reduced by underweighting the 
doppler measurements whenever unusual angular motions occur. 


6-19 


THE ANALYTIC SCIENCES CORPORATION 


A possible method for reducing vertical velocity errors of all 
the systems considered is to incorporate a vertical ac celeration error 
state o This suggestion does not follow from the results presented in this 
report,, A vertical acceleration state was used in the NASA/Ames 
INS/PRS aircraft landing test program (Ref. 21). It was also included in 
the Shuttle landing design study discussed in Ref. 14 — also based on a 
set of PRS (Precision Ranging System) transponders aiding an inertial sys- 
tem. The purpose of this state is to compensate for the effects of vertical 
accelerometer errors and gravity anomaly errors. It is expected that this 
state would be useful, also, in compensating for the vertical acceler om- 
eter misalignment (nonorthogonality) because of the fairly steady decelera- 
tion (approximately 0. 2g) during the terminal approach. With a state cor- 
relation time of about 2U seconds there should be a good opportunity to ob- 
serve this combined error source while the accurate radar altimeter and 
glide slope measurements are being filtered. Systems A, C or D could 
create such a state simply by reducing the correlation time associated with 
the vertical accelerometer bias error state. System B would, of course, 
require an additional filter state . 

As stated in Section 6 . 1 the radar altimeter accuracy will deter- 
mine the ultimate vertical navigation performance achievable by the kinds 
of systems evaluated in this study. The sensitivity curves show, 
however, that none of the four systems studied can meet the preliminary 
vertical velocity spec (0.2 ft/sec), no matter how accurate the radar 
altimeter is if the filter designs are not modified. (The curves plotted 
are the ’’fixed filter ’’ type of sensitivity curves.) Thus, if better accu- 
racy can be obtained in radar altimeter measurements, it is necessary 
to use lower measurement variances, etc. in the filter in order for it to 
place sufficient weight on this good data. In this instance the critical 
factor is not so much a matter of which states are used, but the 



THE ANALYTIC SCIENCES CORPORATION 


parameter values associated with certain state or measurement sta- 
tistics. 

A similar comment can be made concerning the parameters 
associated with the baro altimeter state used in Systems A, C and D. 

The navigation error time histories for Systems A and B reveal an 
advantage for System B in the trajectory portion just prior to the start 
of the radar altimeter measurements. As discussed in Section 6. 1, 
System B apparently makes better use of baro-altimeter data in this 
region. System A would probably benefit from using a more realistic 
expression and values for its correlated altimeter error state. For 
example, an expression combining a bias term and a scale factor term 
(error proportional to altitude) should result in a more proper weighting 
of data than that of Eq. (4.1-10) if reasonable parameters, such as those 
used in the truth model, are used. 

In summary, it appears that software modifications alone cannot 
cause the systems studied to satisfy all of the preliminary specifications 
tabulated in Table 1.2-1. None of the candidate systems can meet the 
vertical velocity spec unless the radar altimeter accuracy is superior to 
that assumed in our truth model. Given sufficient radar altimeter accu- 
racy, the filt er design details discussed above will be very important in 
achieving performance within spec „ Similarly, the conventional navaid 
s ystems are inherently limited in achievable crossrange position accu- 
racy by the correlated errors of the ILS localizer, but careful atten- 
tion to filter design details is necessary to keep other contributions to 
crossrange error small by comparison. In the case of the single -station 
one-way doppler syste m it appears that software improvements alone are 
poten tially capable of reducing rms crossrange error at touchdown to 
within the specification value. 
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1 7.1 SUMMARY OF FINDINGS 

{ Four candidate landing navigation systems for the Space Shuttle 

Grhiter have been evaluated in detail. Three of them are conventional 
1- navaid systems using combinations of DME, 'V OR, ILS and baro and radar 

altimeters aiding a gimbaled inertial system through a Kalman filter. One 
J; . employs a single one-way doppler beacon and baro and radar altimeters, 

als o aiding a gimbaled inertial system through a Kalman filter . Compre- 
•1 h hensive truth models representing all potentially significant sources of 

error have been formulated and used to generate detailed error budgets 
J and sensitivity curves. Table 6 . 1-1 summarizes overall performance at 

| i . touchdown for the four candidate systems. 

The major findings, based on interpretation of the detailed 
results, are as follows: 

e Vertical velocity and crossrange position errors 

at touchdown present the only difficult challenges 

to the design of a system of sufficient terminal 
accuracy,, The projected errors for all four can- 
ditate systems are out-of-spec in these two 
components. 

W 

• The vertical velocity problem has been made much 
more difficult by the steeper approach path of the 
new aerodynamic configuration (assumed in the 
fourth candidate system evaluation). The ultimate 
vertical performance, for the kinds of systems 
studied herein, depends on the radar altimeter 

y. accuracy achievable. RMS measurement errors 

of less than 1. 0 ft are required in the terminal 
phase. 
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» For the conventional navaid systems the crossrange 
performance at touchdown is limited by the ILS 
localizer errors. RMS errors of less than 
1.0 milliradian are required — less than 0,5 milli- 
radian if the error is a true bias or slowly varying 
correlated error. 

* It appears that the single -station one-way doppler 
system could be made to satisfy the crossrange spec 
by software modifications alone. However, given a 
"straight-in" approach this system may not work 
well at all. On the other hand, a multi-beacon 
doppler system should work well with any trajectory 
and provide better accuracy as well. 

* There are significant advantages in employing a 
single VOR/DMEstation rather than two DME 
stations. These advantages relate to the filter 
initialization phase and the "midcourse" approach 
phase, not primarily to touchdown accuracy. 

There are significant advantages in employing the 
drag -update procedure prior to the end of blackout 
(assumed in the fourth system evaluation). This 
procedure results not only in smaller initial errors, 
but prevents erroneous azimuth estimation errors 
related to initial correlations existing prior to the 
drag -updating. -■ 

A better location for the one-way doppler beacon 
would be 6,000 to 10,000 feet from the nominal 
touchdown point, rather, than the 3280 feet assumed 
in the evaluation described in this report. This 
would cause the sudden large reduction in down- 
range error to occur at an earlier time. 

The initial azimuth misalignment at the start of 
the entry phase should have an rms value of 
2.0 milliradians or less for the conventional 
navaid systems and 1.0 milliradians or less for 
a single -station one-way doppler system — in 
order to hold the XMU contribution to terminal 
crossrange error to 2.0 feet or less. 
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A reduction in the rms accelerometer misalign- 
ment error to something less than the 40 sec 
assumed herein would be helpful in reducing 
terminal vertical velocity, errors for all the 
candidate systems studied, ; , - . 

1 A number of comments and suggestions regarding the choice of 

filter states and filter design parameters are given in Section 6.3. Com- 
parison of the three conventional navaid systems studied indicates that a 
[ good sub-optimal filter for the approach and landing phase would include 

approximately 12 states. A similar number is probably appropriate for 
| the one-way doppler system, although additional states related to accel- 

1 erometer errors may be necessary. In general, software m odifications. 

| alone cannot cause the systems studied to satisfy all of the preliminary 

specifications tabulated in Table 1 ,2-1. But, if sufficient sensor accu- 
| racy, especially in the radar altimeter measurements, is available, care- 

ful attention to the filter design details discussed will be vital in achieving 
I the full potential of the overall system hardware and in avoiding various 

pitfalls uncovered in the course of this study, 

7.2 RECOMMENDED FUTTJRE STUDIES 

It is recommended that similar detailed error budgets be 
generated for additional candidate systems of interest. These can include, 
for example, multi-station one-way doppler systems, microwave landing 
J. systems and navigation satellite based systems. The studies should be 

extended backwards in time to include the blackout portion of the entry 
trajectory; this is necessary to obtain a proper accounting of 

correlations between all truth model states during drag-updating. Similar 

studies related to bther Space Shuttle mission phases such as boost and 
rendezvous are also desirable. 
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The generation of detailed error budgets presented in this report 
provide a useful starting point for a filter design improvement effort. 

Given such an error budget based on a comprehensive truth model of about 
80 states, one can define a "reduced state" approximate truth model of per- 
haps 30 states which account for most of the total system error. This 
reduced state model can then be the basis for a set of experimental cal- 
culations in which filter states, filter design parameters, ground antenna 
locations, trajectory, etc. are varied. Once a "good" design is approached, 
the designer can verify its performance using the full 80 staj;e truth model. 

It is recommended that such a filter improvement effort be conducted for 
whatever system or systems emerge as leading candidates. 
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